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Résumé. — A l’aide d’une décharge à fréquence radio passant dans un courant d’oxygène à la 
pression atm sphérique, il a été possible, en utilisant un spectrographe à grande dispersion, d'étendre 
l’analyse de vibration et de rotation des bandes du système de Schumann-Runge (B 32 — X 3%) 
de O, jusqu’à 5 700 À. Nous avoas été également en mesure de mettre en évidence une prédisso- 
ciation à v' = 3 et v’ — 4. Le niveau v’ = 4 est fortement prédissocié et il ne semble pas y avoir 
de dissociation appréciable à D’ = 5. 

Dans ces conditions, nous observons dans la même décharge quelques nouvelles bandes atmo- 
sphériques. Les séquences Av = + 1 et Av — 0 montrent chacune 5 bandes. Nous suggérons 
que l’excitation des bandes atmosphériques est due à l’interaction entre les états Æ 3X et b1Z; 
elle serait favorisée par la pression élevée du gaz. 


Abstract. — In a high frequency discharge through flowing oxygen at atmospheric pressure, 
it was possible, with high dispersion spectrographs, to extend the vibrational and rotational 
analysis of the Schumann-Runge bands (B 52 — X #E system) up to 5 700 À. We have been 
able also to show a predissociation at v” — 3 and 4. The level v’ — #4 is strongly predissociated. 
There seems to be no predissociation at 0’ = 5. 

In the same discharge, we have observed some new atmospheric bands. Each of the Av = +1 

. and Av — 0 sequences have five bands. We suggest that the excitation of the atmospheric 
bands is due to an interaction between the X °Z and b ?X states enhanced at high pressure. 


Introduetion. — On sait que les raies aurorales 
interdites de l'oxygène [07]6 300 — 6 364 peuvent 
être excitées par dissociation du niveau supérieur 
du système de bandes de Schumann-Runge de O, 
(Cabannes et Garrigue, 1936; Herman et Weniger, 
1950a ; Chamberlain, 1958). En raison de son 
importance pour la physique de la haute atmos- 
phère, nous avons repris l’étude de ce système avec 
une dispersion plus grande. Nous avons utilisé 
un générateur de haute fréquence (45 Mz) et un 
courant d'oxygène à la pression atmosphérique. 

L'aspect du spectre obtenu dépend du débit 
d'oxygène et du voltage employés. Pour un fort 
débit et un fort voltage, on obtient un spectre 
atomique assez intense comprenant les raies de O7 
et les raies intenses de OZ], superposées aux bandes 
de Schumann-Runge. Au fur et à mesure que l’on 
diminue le voltage, l’intensité des raies atomiques 
permises diminue et l’intensité relative des raies 
interdites et des bandes de Schumann-Runge aug- 

(*) Ce travail a été subventionné en partie par Air Force 
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mente (fig. 1). Cette série d'expériences ayant mon= 
tré les conditions favorables pour l’émission des 
bandes de Schumann-Runge correspondaient à 
une faible excitation, on a fait varier le débit d’oxy- 
gène. L’intensité des raies interdites varie peu ; par 
contre, on constate que, sans débit d'oxygène, 
c’est-à-dire lorsque la décharge a lieu dans l’air, 
on observe les séquences Av = 2,3 et 4 du premier 
système positif de l’azote. 

Lorsqu'on emploie un faible débit, seules sub- 
sistent les bandes les plus intenses de ce système 
et les bandes atmosphériques commencent à appa- 
raitre. Lorsqu'on augmente encore le débit, les 
bandes de N, disparaissent et les bandes atmosphé- 
riques deviennent intenses. Lorsqu'on augmente 
l'intensité du courant, il apparaît un fond continu 
d'émission. 

Nous étudierons tout d’abord les bandes de 
Schumann-Runge. 


Système de Schumann-Runge. — Ce système 
montre une structure compliquée s'étendant du 
rouge à l’ultra-violet. Les bandes disparaissent 
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région visible par Lochte-Holtgreven et Dieke 
(1929) jusqu’à 4 350 À. Nous avons pu les photo- 
graphier à l’aide d’un spectrographe à 3 prismes 


au-dessous de 2 100 À en raison de l'absorption par 
des couches extérieures plus froides de l’étincelle. 
La structure des bandes avait été étudiée dans la 
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Fic. 1. — Ces 3 spectres ont été obtenus en diminuant la chute de tension dans l’arc de a vers bet c. 
Où voit que les raies interdites [OZ] et les bandes de Schumann-Runge augmentent alors d’intensité. 


de flint donnant une dispersion de 7,5 Â/mm vers 
7 000 À. L’analyse a été possible jusque vers 
& 900 À ; au delà, la dispersion devient trop petite 
en même temps que la superposition des bandes 
augmente, la parabole de Franck-Condon étant 


large. Nous avons utilisé alors un spectrographe 
à réseau échelle donnant une dispersion de 
2,5 À/mm. Dans ces conditions, l’analyse a pu 
être poussée jusqu’à 5 700 À. Nous donnons dans 
le tableau I la liste des bandes nouvelles analysées. 


TABLEAU I 


LONGUEURS D’ONDE DES ORIGINES DES BANDES DU SYSTÈME DE SCHUMANN-RUNGE QUE NOUS OBSERVONS. 
(Les valeurs entre parenthèses ont été calculées. Elles correspondent à des bandes faibles.) 


\v’ 0 1 2 3 A 5 6 
» 
20 4 422,9 
21 4 503,5 & 372,6 
22 4 586,6 4 455,6 
23 4 816,8 & 672,4 
24 5 064,8 4 905,3 (4 761,3) 
25 5 332,1 5 155,7 (4 996,8) 
26 (5 620,7) (5425,0) (5 249,3) (5 091,5) 
27 (5714,5) (5519,95) 5 345,61 
28 (5 671) (5 466,1) 
29 (5 719,8) 


Le tableau IT donne les longueurs d’onde et les 
fréquences des raies mesurées. Les valeurs des 
constantes moléculaires B’, B", D', D’ ont été déter- 
minées par la méthode graphique indiquée par 


où D, est calculé par 
les formules 


Dy = De + Bel + 1/2) 


&PBe 
Herzberg (1950) en portant (X + 1/2)? en abscisses De "et 
et K +1 ordonnées (colonnes 2 et 3 des ta- FER (ee _ nn) hd 
bleaux ITT et IV) et par la méthode de Birge (co- PU ee 
lonne 5, tableaux III, IV) qui donne Re 


Ge Oe < 
&By = E[A, F(K) + 8D,(K + 1/2)3] (K + 1/2) 
EX + 1/2)? 


Les valeurs obtenues pour D, figurent dans la 
colonne 4, tableaux III et IV. On a en outre déter- 
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miné AG (et x) par la méthode graphique décrite TABLEAU II (suite) 
par Herzberg en portant [R,,.(KX) — R,,.(K)] 


Où [Row (K — 1) — Ryy(K — 1)] en fonction de Æ VR vP A3 F” A, F° 
KÇK + 1) (tableau V). 57 898,30 728,70 173,60 242,55 
59 27,15 648,11 179,04 250,19 

TABLEAU II 61 759,92 

63 676,09 


LONGUEURS D’ONDE, 
NOMBRES D'ONDE ET DIFFÉRENCES DES TERMES 
POUR LES BANDES. (1,21), (2,21), (2,22), (2,23), (2,24), BANDE (2,22) 
(2,25), (3,23), (3,23), (3,24) 


1 21 799,64 21 795,70 3,94 
BanDe (1,21) 3 99,64 87,67 11,97 11,97 

ee : NE s des De 97,55 80,17 17,38 19,47 
Ke Si Le 2e 7 92,75 68,84 23,91 28,71 

1 9 85,39 85,66 29,73 37,09 

3 11 75,85 39,34 36,51 46,05 

5 13 63,62 21,35 42,27 54,50 

7 15 48,95 00,93 48,02 62,69 

9  22187,0 29 158,4 29,0 : 17 32,40 677,36 55,04 71,59 
11 76,6 181,7 34,9 15,7 100049 13,10 51,96 61,14 80,44 
13 65,8 121,6 &4 2 550 21 21 690,483 21 622,974 67,509 90,13 
45 51,5 102,1 49,4 63,7 23 665,710 592,898 78,812 97,585 
È 3n 8 078.6 56,2 729 25 639,027 560,346 78,681 105,364 
Fa ne is De D DEC 610,447 525,344 85103 113,683 
À TRE ee a 29 579,474 489,348 90,126 121,099 
æ 068.9 21 994 1 71.8 992 31 545,664 447,671 97,993 131,803 
Da 0611 960.3 80.8 1086 33 509,506 405,423 104,083 140,241 
2 0111 925.6 85,5 155 35 470,881 361,059 109,823 148,447 
Es 21 981 5 887.2 9% 3 1239 37 429,604 914,182 115,422 156,700 
9u6 3 847.0 993 1385 39 386,421 264,690 121,731 164,914 
x 908.3 8032 1051 2 ete 340,586 242,847 127,739 173,574 
e 871 4 755 8 115.6 1525 43 292,362 158,061 134,301 182,525 
à 830.3 708.5 121 8 160. DE 241,648 103,543 138,105 188,819 
D 786 1 652.6 1335 1777 47 189,422 044,682 144,740 196,966 
… 70.6 49 134,681 20 983,417 151,264 206,005 
# 692.0 51 076,582 920,450 156,132 214,231 
# 53 21 016,721 * 920 854,509 162,212 222,073 
47 
49 BANDE (2,23) 
51 1 

3 
BANDE (2,21) 5 20 756,27 

33 22 57,37 22 53,76 103,61 142,05. 7 20 752,86 20 729,20 23,66 27,07 
35 16,42 06,96 109,46 15041 9 20 745,54 20 716,27 29,27 33,59 
37 472,95 357,31 115,64 159,11 41 20 736,98 20 701,77 35,21 43,77 
39 26,82 04,93 121,89 168,02 13 20 726,02 20 684,04 41,98 52,94 
41 378,32 250,38 127,94 LOTO: 20 712,88 20 664,28 48,60 61,74 
43 26,90 193,57 133,33 184,75 17 20 697,18 20 642,27 54,91 70,61 
45 273,31 33,92 139,46 192,98 19 20 678,92 20 617,79 61,13 79,39 
47 17,15 071,81 145,34 201,57 21 20 658,82 20 591,20 67,62 87,72 
49 158,15 07,48 150,67 209,67 23 36,22 63,08 73,14 95,74 
51 096,84 21 940,49 156,35 217,66 925 11,07 31,96 79,11 104,26 
53 33,09 870,65 162,44 22619 27 20 584,25 20 499,12 85,13 111,95 


D9 21 967,25 98,59 168,66 234,50 29 595,26 63,75 91,51 120,50 


39 


43 


© I Où ww = 


VR 
23,57 
20 490,05 
53,83 
16,06 
20 375,95 
32,57 
20 287,55 
20,56 
20 190,77 
38,95 
20 085,22 
29,21 


19 688,72 
673,64 
654,81 
634,04 
611,09 

19 585,92 
558,71 
529,64 

19 498,67 
464,84 
429,02 

19 390,40 
350,70 
309,16 

19 265,41 
218,80 


18 454,8 
130,8 
393,8 
354,7 
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TABLEAU II (suite) 


VP 


26,45 
20 386,88 
44,85 
20 300,03 
20 253,85 
04,93 
20 153,85 
20 100,89 
20 045,75 
19 988,20 
28,79 


BANDE (2,24) 


19 690,07 
674,65 
656,27 
635,70 

612,39 

19 587,45 
560,55 
531,74 
500,72 
467,23 
431,36 
392,69 
353,68 
312,18 
268,67 
222,95 
175,28 
125,00 
073,69 


BANDE (2,25) 


18 349,6 
309,3 
264,8 
221,5 


BANDE (3,22) 


A, F' 


97,12 
103,17 
108,98 
116,03 
122,10 
127,64 
133,70 
139,67 
145,02 
150,75 
156,43 


53,02 
61,25 
67,36 
73,49 
79,35 
85,20 
91,48 
98,28 
105,98 
111,16 
116,84 
121,73 
127,75 
133,88 
140,41 
145,11 


115,2 
121,5 
129,0 
133,2 


76,33 

86,19 

94,26 
102,30 
110,37 
118,69 
127,35 
136,95 
144,99 
152,66 
160,35 
167,45 
175,42 
184,16 
191,72 


155,5 
166,0 
172,3 


K 
11 


13 
15 
17 
19 
21 
23 
25 
27 
29 
31 
33 
35 
37 
39 
41 
43 
45 
47 
49 


© 1 Ov © 


11 
13 
15 
47 
19 
21 
23 
25 
27 
29 
31 
33 
35 
37 
39 
k1 
43 
45 
47 
49 
51 
53 


VR 


22 414,8 
103,2 
389,1 


349,1 
327,6 
299,3 
273,9 


208,6 
172,7 
134,9 
093,3 
052,3 
006,2 
21 960,3 
905,0 
852,8 
798,5 
740,6 


21 376,6 
363,4 
350,6 
336,0 
316,0 
292,9 


249,6 
214,2 
181,4 
149,2 
113,9 
076,0 
033,9 
20 992,4 
945,3 
898,9 
848,4 
796,8 
741,1 
683,6 
623,0 


TABLEAU IT (suite) 


VP 


22 379,6 
360,1 
340,1 
316,7 
290,8 
259,0 


194,8 
158,4 
121,6 
078,6 
033,9 
21 986,4 
937,2 


A 
35,2 
43,1 
49,0 


58,3 
68,6 


BANDE (3,23) 


21 279,3 
253,9 
227,8 
196,3 
161,3 
130,4 
091,2 
055,2 
012,7 

20 968,2 
921,7 
872,7 
821,9 
767,0 
713,6 


56,7 
62,1 
65,1 


81,3 
83,8 
90,2 
94,0 

101,2 

107,8 

112,2 

119,7 

123,4 

131,9 

134,8 


N°1 : 


71,3 
82,1 
88,2 
96,6 


112,2 
123,0 
126,2 
136,5 
445,7 
154,3 
161,2 
171,0 
178,3 
185,3 
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TABLEAU II (suite) TABLEAU IV 


CONSTANTES DE ROTATION B% ; De 


v 9 1 
Baxne (3,24) POUR LE NIVEAU INFÉRIEUR À 92 


K INSNRE INQNRE » » » ” 
ES ve ns 2 2 BANDE Bo GRAPH. D, GRAPH. Des CALC. Bo BIRGE 
} Hs 1E ee Es LE 
3 (2,21) 1,105 5,3.10—6 4,95.10—% 1,104 
5 14 L4 LA 11 
7 Ba graph. D39 Da Dos 
9 (2,22) 1,091 h,9.10—6  4,95.10—8 . 1,085 
a1 (3,22) 1,085 
13 ” ” ” ” 
45 Pos Das Dos Bas 
L (2,238) 1,068 41.10—6 5,0 .10—% 1,069 
19 (3,23) 1,082 10—5 1,068 
21 By D Dos Ba 
M 0,120 266,0 (2,24) 1,055 6.10—5 5,0 .10—5 1,049 
25 243,0 20 166,3 76,7 99,7 (3,24) 1.055 
27 214,8 132,5 82,3 110,5 . : ; "1 
29 186,6 098,1 88,5 116,7 Ps Das D; Ps 
31 155% 061,0 94,7 125,6 (2,25) 1,033 
Ah) 122,4 020,9 401,5 134,8 
35 085,6 19 976,1 109,5 146,3 TABLEAU V 
37 048,2 929,5 118,7 156,1 & | 
39 009,7 887,8 121,9 160,4 ALEURS OBSERVEÉES ET CALCULÉES DE n + 1/2 
44 19 966,9 841,0 425,9 168,7 D'u— 2 ef CALCULÉ 
: 6,3 me = _ 
43 291 790,6 432,5 476, At 1 066,9 1 067,5 
NO 1 040,6 1 043,0 1 042,9 
A Go3,5 1 011,8 ANOMALIE 4 016,2 
Les bandes (1,22) et (2,25) n’ont pu être iden- 
TABLEAU III tifiées qu’en partie. La bande (1,22) a la plupart 
de ses raies « blend » avec celles d’autres bandes 
ConNSTANTES DE ROTATION Be ET Do et est assez faible. La bande (2,25) est « blend », 
POUR LE NIVEAU SUPÉRIEUR B SII en partie, avec la bande (5,27) 
Baxpe B,Grapn. Di Grapn. Ds cac. B: Bircr Prédissociation. — On n’a pas donné de bandes 
4 2. es ee >. RARE : RES me ; 
(2,21) 0,786 4,8.10—5 USE 4. Celles-ci sont en effet très difficiles à 


mesurer en raison de leur manque de définition. 


2,22 0,786 4,8.10—$ 0,784 8? : 

ivre re _ A EAN Ne Les bandes à #” — 3 sont beaucoup plus diffuses 

A et de as jm que les bandes à v’ — 2. Les bandes à #9’ — 4 
? ? à ? fr 7 


Fe sont si diffuses qu’on ne peut pas les mesurer. Un 

moy. 0,786 5,0,10—% moy. 0,785 résultat analogue a été trouvé également dans 
K x l’ultraviolet par P. K. Caroll (1959). D’après cet 
auteur, la prédissociation commencée à p' = 3, 

B3 cRapu. DS GRAPH. DS cac.  B;Birce est très forte à v’ — 4 et beaucoup moins marquée 
à pe" — 5. Nos résultats sont en parfait accord avec 


ER 0,769 ces remarques 

(3,23) 0,769 4,76.107° 0,778 De nombreux travaux ont été faits sur cette 

(3,24) 0,721 prédissociation. Flory (1936) avait signalé une 

(3,25) 0,780 prédissociation à p’ — 3. Feast (1949 ; 1950) avait 
no 07 Fan rejeté l’hypothèse de Flory sous le prétexte qu’il 


avait observé des bandes à p” — 3 en émission. 
Nous avons également identifié des bandes, en 
émission, à #” — 3 et p’ — 4 dans l’ultraviolet 
(Rakotoarijimy et autres, 1958). Ceci s’explique 
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aisément du fait que l’élargissement s’exprime en 
em, or, un élargissement de 10 cm1, par 
exemple, autour de 2 500 À correspond à 0,6 À 
alors qu’à 5 000 À, il correspond à 2,5 À. Dans nos 
conditions expérimentales, nous avions une dis- 
persion de l’ordre de 5 À par mm dans l’ultraviolet, 
ce qui donnait un élargissement de 1/10 mm. A 
5 000 À, nous avons une dispersion de 2,5 À 
par mm, ce qui donne un élargissement de 1 mm. 
Il est évident que, dans nos conditions expérimen- 
tales, les raies qui n’avaient aucune apparence par- 
ticulière dans l’ultraviolet sont, par contre, très 
élargies dans le visible. 


e Aydberg 


o Morse 
* Hulburt-Hirschfelder 


PA) 1e 207 15 ? c 


F1G. 2. — Courbes de potentiel 
pour les niveaux >, et %>, de O.. 


Cette prédissociation avait été prévue par 
P. G. Wilkinson et R. S. Mulliken (1957). Pour 
ces auteurs l'intersection de la courbe instable 
aurait lieu vers 3,5 alors que, pour Caroll, elle 


N°1 
aurait lieu sensiblement pour #’ — 4. Nos obser- 
vations semblent montrer que le point d’intersec- 
tion est au voisinage de p’ — 4, mais situé entre 
3 et 4, sensiblement vers 3,75. D'ailleurs, on 
devrait pouvoir préciser davantage cette position. 
D’après Caroll, le niveau provoquant la pertur- 
bation serait un niveau 3J],. La transition serait 
permise, ce qui expliquerait l’intensité de la per- 
turbation. 

Sur la figure 2, nous avons tracé les courbes de 
potentiel pour le niveau inférieur X 32, et le 
niveau supérieur B 327. Pour le niveau inférieur, 
nous avons utilisé la formule de Morse et la mé- 
thode de Hulburt et Hirschfelder. Pour le niveau 
supérieur, la formule de Morse et les résultats de 
Rydberg. On voit qu’il est vraisemblable que la 
courbe donnée par la formule de Hulburt — Hirs- 
chfelder soit plus approchée de la vraie courbe de 
potentiel que celle de Morse par comparaison avec 
les courbes relatives au niveau supérieur. On voit 
d’après cette figure qu’on doit s’attendre, par appli- 
cation du principe de Franck-Condon, à un 
maximum d'intensité pour les bandes (2,23) ; 
(3,24) ; (4,26) ; (5,27) ; (6,29). Nous avons égale- 
ment tracé, sur cette figure, l’état hypothétique 
répulsif 5IT, qui donnerait la prédissociation du 
niveau B 527 à v’ voisin de 4. Dans la région 
située autour de 5 400 À, d’après le tableau cal- 
culé à partir de la formule 


v = 49 357,55 + (698,86 v’ — 10,708 0? — 0,194 vi) 
— (1 568,33 2” — 11,983 v’? + 0,0517 v’® — 0,00143v’4) 


donnée par Curry et Herzberg (1934), les bandes 
possibles sont : (2,25) ; (3,26) ; (4,26) ; (5,27) ; 
(6,27) ; (6,28). Les seules raies distinctes s’inter- 
prètent comme (2,25) et (5,27). Les autres bandes 
possibles (3,26) et (4,26) se manifestent par des 
paquets flous répartis tout le long du spectre dans 
cette région. L’étude des régions de plus grande 
longueur d’onde se poursuit. 


Bandes atmosphériques. — En même temps que 
les bandes de Schumman-Runge, lorsqu'on aug- 
mente le débit d’oxygène, on obtient 3 séquences 
de bandes assez bien développées. On peut les 
identifier aux bandes atmosphériques 


DIS RCE 


de O, comme on peut le voir dans le tableau IV 
dans lequel on a porté les longueurs d’onde des 
têtes mesurées et les longueurs d’onde des têtes 
calculées. Les écarts entre les valeurs calculées et 
les valeurs mesurées sont de l’ordre des erreurs de 
mesure, la dispersion variant de 110 à 200 À par mm 
dans cette région. Les bandes en italique n’ont pas 
été observées précédemment. La limite du système 
à 8 698 À est purement instrumentale et due à la 
chute de sensibilité des plaques utilisées (ZN Kodak 
non hypersensibilisées) (tableau VI). 
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TABLEAU VI 


LONGUEURS D’ONDE MESURÉES ET CALCULÉES DES TÊTES DES BANDES ATMOSPHÉRIQUES b 12 — X 32 OBSERVÉES 
DANS UNE DÉCHARGE DANS L'OXYGÈNE SOUS PRESSION ATMOSPHÉRIQUE 


(HN ENCALE 2 OBS. (0 De A CALE MN oBS ol (0, D'IMAGE À OBS. 
(1,0) 6868 6866 | (0,0) 7589 7598 | (0,1) 8583 8598 
(2,1) 6955 6953 | (1,1) 7688 76890108 (1.9) 8685 8698 
(3,2) 7 046 7 043 (2,2) O1 
(4,3) 7 141 "141 | (3,9 7879 7876 
(5,4) 7 242 7 240 (4,4) 7 984 7 987 
(5,5) 8096 8103 
(1,0), (2,0), (3,0). Il est remarquable que dans le 
eo À] cas de X, on observe les raies interdites 5 577 
_—— et 2 972 de OZ. Dans l’oxygène pur, R. Herman 
——……— et S. Weniger (1950b) avaient trouvé un maximum 
se d'émission pour les bandes atmosphériques à des 
pressions de quelques dizaines de mm de mercure. 
Elles étaient accompagnées d’un spectre OZ et 
_- 6866 (1,0) d’un fond continu allant de 9000 à 2 500 À. 
6953 (2,1) Pour des pressions de l’ordre de 1 atmosphère, on 
MT eue ‘ n’avait plus que les raies de OZ et de OJZT et un 
2 | --7043 (8,2) fond continu d’origine moléculaire. On a pu mon- 
21-7140 (43) trer, dans ces expériences, que l’excitation était 
= 7240 (5,4) électronique et que les bandes atmosphériques 
| d et les bandes de Schumann-Runge étaient excitées 
par des chocs de faible énergie. Nous retrouvons 
également le fond continu signalé dans les expé- 
7598 (0,0) riences ci-dessus lorsqu’on augmente la densité de 
- 7689 (1,1) courant. 
i- 71774 Oy+(2,2) 
… 7876 (3,3) 
H,0 7937- PL- 7987 (4,4) 
£ - 8103 (5,5) 
H,0 8227 -4 
Or 8446- 
4 -- 8598 (0,1) 
_=- 8698 (1,2) 
Fic. 3. — Bandes atmosphériques observées dans une 


décharge dans l’oxygène sous la pression atmosphérique. 


On observe en outre, comme on peut le voir sur 


la figure 3, trois bandes de la vapeur d’eau 6 994 À, 
7 957 À, et 8 227 À. 

On sait que les bandes atmosphériques corres- 
pondent à une transition interdite et sont difli- 
ciles à observer en émission. Toutefois, Kaplan 
(1947) les avait obtenues dans un mélange 
N, +O,, R. C. Herman, G. R. Hornbeck et 
S. Silverman (1949) dans un mélange CO + O. 
R. et L. Herman (1949) ont observé, en présence 
de He et de X, que l’intensité relative des bandes 
atmosphériques augmente avec la pression. Les 
bandes observées sont (0,1), (1,2), (0,0), (1,1), (2,2), 


F1G. 4. — Schéma de l’arc haute fréquence montrant les 
différentes couches d’émission et d’absorption. 
Les raies OZ et OT sont émises dans la zone 41 de 
forte excitation ; les bandes de Schumann-Runge sont 
émises dans la zone 2 et réabsorbées dans la zone 3. 
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TABLEAU VII 


NOMBRES D’ONDE EN CM—1 DES NIVEAUX DE VIBRATION DES ÉTATS ÉLECTRONIQUES IZÿ ET °Zg 


pra Dr 
U Ty v F, AE 
0 13 908 0114819069 0 VE 
1 45.914301 10: 45909. PU0E 
2 16 690 | 414: 16,636. 54 
3. 418 039 | 42. 17 940 99 
& 19360 | 13 19221 139 


Diseussion. — Dans la décharge obtenue à l’aide 
du générateur R. F. on observe 3 zones représentées 
sur la figure 4 : la zone 1 parcourue par l’étincelle 
est celle où a lieu l’excitation et la décomposition 
des molécules (c’est cette région qui donne nais- 
sance aux raies OZ et OT) ; dans la zone 2 a lieu 
la recombinaison, on y observe surtout les bandes 
de Schumann-Runge ; dans la zone 3 a lieu la réab- 
sorption. 

Bates et Dalgarno (1953) ont montré que la 
réaction 


O(D)+O{X 527, v = 0)—O(8P,) + O,(b 1,0 — 2) 


peut avoir lieu avec une grande probabilité parce 
que le bilan d’énergie est très favorable. Or, nous 
avons ici une grande concentration d’atomes 11) 
mais pas de concentration importante de molé- 


DE O» 

IT Sa | 

a) 15; U Fe, AE | 
5 20653 | 14 20480 + 139 
6210018 4500210716 1738 
7 0198 156 | 16:22 928 202 
822565 MEME M11 248 
9 195 547 | 18 25 283 264 
10026 700710192026 49% 276 


cules O, (b 12+ v — 2). Il est donc peu probable 
que cette réaction ait lieu. Par contre, lorsqu'on 
compare les énergies des niveaux de vibration des 
états b 12# et X 827 on constate que, sur un 
grand intervalle, les niveaux de vibration de ces 
deux états sont assez voisins, la distance étant 
inférieure à 0,03 eV, comme on peut le voir sur 
le tableau VII. On constate que la meilleure con- 
cordance (inférieure à 0,02 eV) a lieu entre 9 et 14 
pour les quanta de vibration de l’état fondamental. 
Ces niveaux sont en effet très peuplés dans nos 
conditions expérimentales, il est donc possible que, 
par les chocs nombreux qui peuvent se produire, 
il yait passage de l’état 3277 à l’état b 12f, ce qui 
expliquerait l’extension inhabituelle des séquences 
dans les bandes atmosphériques. 
Manuscrit reçu le 2 août 1960. 
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SUR LE CALCUL VARIATIONNEL DES PROPRIÉTÉS MOLÉCULAIRES 


Par Françoise DUPONT-BOURDELET, Jacques TILLIEU et Jean GUY, 
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et Groupe de Chimie Théorique de l’École Polytechnique. 


Résumé. — Un formalisme général est exposé pour le calcul variationnel des propriétés molé- 
culaires (électriques et magnétiques principalement), de manière à pouvoir tenir compte de plu- 
sieurs perturbations simultanées. Les énergies sont calculées explicitement jusqu’au 4e ordre et 
l’on montre qu’elles ne font intervenir que les fonctions perturbées jusqu’au 2e ordre. Moyennant 
certaines hypothèses, les équations et expressions obtenues sont décomposables (souvent sous 
forme simplement additive) soit par rapport aux divers opérateurs perturbateurs, soit par rapport 
à des groupes indépendants de particules (électrons et noyaux, par exemple). 


Abstract. — À general formalism is set up for the variational calculation of molecular pro- 
perties (especially, electrie and magnetic) in order to account for several simultaneous per- 
turbations. Energies are explicitly calculated up to the fourth order and it appears that the 
formulae only contain the perturbed functions of the first and second orders. By means of cer- 
tain hypotheses, the equations and matrix terms can be separated (often in a simply additive 
form) either to different perturbating operators or to different independent groups of particles 


TOME 22, JANVIER 1961, PAGE 9. 


(e.g. electrons and nuclei). 


I. Introduction. — Le nombre des applications 
de la méthode variationnelle au calcul des diffé- 
rentes propriétés des atomes et des molécules aug- 
mente continuellement. Outre un article général 
déjà ancien de Gans et Mrowka [1] qui semble être 
passé inaperçu, citons, sans prétendre faire une 


bibliographie complète, quelques articles sur le 


calcul des polarisabilités électriques [2], des suscep- 
tibilités magnétiques [3] et des différents effets 
intervenant en résonance paramagnétique nu- 
cléaire [4]. Dans ces conditions, il nous a semblé 


intéressant de présenter cette méthode sous une 


forme assez générale pour être utilisable rapidement 
dans de nombreux cas, même plus compliqués que 
ceux envisagés jusqu'à maintenant et de discuter 
diverses approximations que l’on est amené à uti- 
liser au cours des applications. Partant d’un sys- 
tème moléculaire dans un état stationnaire non 
perturbé décrit par la fonction d’onde (f. 0.) 4 
ou de manière équivalente, par le ket [0 >, nous 
le soumettrons à diverses actions perturbatrices, 
par exemple un champ électrique et un champ 
magnétique extérieurs, ou encore des termes négli- 
gés en première approximation comme le champ 
dû à un moment magnétique nucléaire. Dans ces 
conditions, l’hamiltonien de la molécule perturbée 
peut être écrit sous la forme : 


H = H(N + eH(9 + e? H(?) 


où « est un paramètre introduit pour la commodité 
des calculs et servant à repérer « l’ordre » de la 
perturbation. 

Dans cet article, nous n’aurons pas besoin de 
préciser plus la forme des hamiltoniens pertur- 
bateurs H® et H®, Remarquons simplement que 


nous supposons que JT est au plus du 2° degré en &. 
Conformément à la théorie des perturbations, 
nous écrivons la f. 0. perturbée sous la forme : 


ge ER ep) et rene) 

ou, en utilisant les vecteurs kets de Dirac : 
Né NO ue > het Re Der Je 
(27) 
Si nous voulons que la fonction 4 soit normalisée 


à tous les ordres en €, nous devons satisfaire aux 
conditions successives 


290 RS AE (3) 
202 AMP AL 
Il est commode de prendre la première condition 
sous la forme plus restrictive : 
OR CA 02 "0 (3°) 


c’est-à-dire de choisir ® orthogonale à 4). De 
même il peut être plus facile de choisir 4® ortho- 
gonale à 49 en imposant : 


PR DR Te UN (4) 
On peut passer d’ailleurs de 4?’ satisfaisant 
à (4’) à une fonction 4‘? satisfaisant à (4) par : 


PR ES ga —< A A7 h00 


et on verra que les relations où figure 4% ne dé- 
pendent pas de la condition utilisée. 
De même l’energie perturbée est écrite 


E = E0 + EU) + et EU + ... + em Et + ... (5) 
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Il est d’ailleurs assez rare que l’on désire aller 
au delà de la perturbation d’ordre 2 pour l’énergie. 
Toutefois, des travaux récents, notamment sur 
l’effet Kerr et les « hyper-polarisabilités » [5], font 
intervenir des énergies d’ordres supérieurs à 2 et, 
dans ce qui suit, nous expliciterons les expressions 
jusqu’à l’ordre 4. 


II. Caleul de l'énergie perturbée. — Nous vou- 


lons déterminer les kets [1 >, [2 >, ... [n > 
de manière que l’énergie : 
__ < yAlg 7 : 
Fe Ne : 


soit minimum. Il revient au même d’annuler l’ex 
pression : ; 


F=<YHY>—E< 9927 (7) 
qui serait effectivement nulle pour la fonction 
exacte, et la vraie valeur de Æ. Nous pouvons 
utiliser ce fait, dans le cadre de la méthode des per- 
turbations, en annulant F aux différents ordres et 
en minimisant ensuite les termes correspondants 
de l’énergie. En développant 4, et E par rapport 
aux diverses puissances de €, en tenant compte 
de l’équation de Schrôdinger pour le système non 
perturbé et de l’hermiticité de A9, soit 


H(®0 > — E%|0 > 
et 
< OHO]j > = < jjHO(0 >* 
= E0 < jl0 >* = EO < 0|j > 


on obtient finalement l’expression : 


F = e[< 0[H@(0 > — Et} + e?[< 1]AH(0 — EN > 


1 41H — EI) > Lol EN 

+ < 0JA® — E%]0 >] + [< 1]H0 — E09 > 

+ < 2/A(0 — Et0f > + < O[H( — EN > 

+ 2H = EU) > + AH EUR 

+ < 0[A® — EN > + <'1[H® — E]0 > -— EG] 
+ #[<A[H E0/3.> + < 3|H(9 — Eh > 

+ << 2]H(0 — Et92 > + < 0JHO — E0B > 

+ < 3/40) — E]0 > + < 1/40 — E02 > 

+ < 2{H0 — EE > + < 0JH® — E2 > 

+ < 2[H — E%10 > + < 14% — EN > | 
— EB(< 0H > + < 110 >) — Et]. (8) 


Si nous voulons que } soit nulle quel que soit e, 
nous devons avoir : 


. EU = < 0|H00 > (9) 


N°1 


Et = < 1[H0 — EN > + < 1H ERI0 > 

+ < 0JHO — EN > + < 04810 > (10) 
EG A HN ENS RS MINS EURE 

+ < 0[HO — E02 > + < 2[H00 — EI0 > 

+ <1[H0 — EU > + < 0[H(® — EN > 

+ <1| H® — E®]0 > (11) 
EU = <14|H(0 — Et0]3 > + < 3[H(9 — Eh > 

+ << 2/H(0 — E02 > + < 0[H® — E03 > 

+ < 8JH0 = EU0 > +'< 1[H0 — EM > 

+ < 2/H0) — EU > + < 0[H® — F2 > 

+ < 2]H% — E%10 > + < 1[H® — Eli >. (12) 


A propos des formules (9) à (12) (que l’on géné- 
ralise facilement pour un ordre n quelconque), nous 
pouvons faire les remarques suivantes : nous re- 
trouvons le résultat bien connu que l’énergie per- 
turbée d'ordre nr ne dépend que des fonctions per- 
turbées jusqu’à l’ordre n — 1 ; (9) est la formule 
habituelle pour l’énergie du premier ordre, Le 
4e terme de (10) est la perturbation du 1% ordre 
correspondant au terme du second degré H‘? de 
l’hamiltonien perturbé. De même, les deux derniers 
termes de (11) correspondent à la perturbation du 
2e ordre pour A, etc... 


III. Calcul de l’énergie du premier ordre. — Le 
calcul de E®% par la formule (9) ne présente pas 
de difficulté particulière puisqu'il ne nécessite que 
la connaissance de 4. Il convient d'introduire 
maintenant certaines notions et approximations 
qui présentent un grand intérêt physique et qui 
serviront également pour les ordres supérieurs. 


A) ADDITIVITÉ PAR RAPPORT AUX OPÉRATEURS 
DE PERTURBATION. — Explicitons la forme des 
hamiltoniens perturbateurs A® et H® afin de 
pouvoir envisager plusieurs perturbations simul- 
tanées. Nous écrirons, de manière convention- 
nelle : 

? 
HO = AAPEN HE (131) 
u=1l 
— 


Le ù 
ù Ày DH ae (A Q à).H 0. (13-2) 


U,vy=1 


Ho) = 


Le dernier membre de (13-1) représente le pro- 
duit scalaire de deux vecteurs appartenant à un 
espace à p dimensions, somme directe des p espaces 
unidimensionnels relatifs aux divers opéra- 
teurs HP ; le 2e membre de (13-1) représente le 

— 


: ER 
produit doublement contracté du tenseur H® par 
le produit tensoriel À @ A. 


Nat 
— 


— 

Puisque H® n'intervient dans l'expression de 
l’énergie que par son produit doublement contracté 
avec le tenseur symétrique À @ À, on peut se 
borner à le considérer comme un tenseur symé- 
trique. En outre, on ne diminue pratiquement pas 
la généralité en prenant des ® réels. 

Les définitions générales (13) sont destinées à 
tenir compte de la présence éventuelle de deux ou 
plusieurs opérateurs de perturbation scalaires, vec- 
toriels ou tensoriels, réels ou complexes. Dans ces 
conditions, les à, doivent, au même titre qu’e, 
être considérés comme des paramètres de pertur- 
bation indépendants les uns des autres. Il convient 
de remarquer que certains /®) peuvent être nuls 
sans que les À correspondants le soient. Par 
exemple, pour un champ électrique nous pren- 
drons : 


re = Éy 
a 


À er ; 


k = 6 À—6 
HO) — — er 


er a = — ey AD —= — 3 


quels que soient u et v. 
Nous pouvons alors écrire l’énergie perturbée 
du 1*% ordre sous la forme 


DTA EN ABC (14) 
He 
(Il convient de ne pas confondre des termes 
comme £® et E), En tenant compte de (13-1), 
l’expression (9) devient 


ED = XX < 0JAD 0 > (15-1) 
LH 


d’où, en comparant à (14) 


ED = < 0|AŸ10 > 


Il en résulte que l’énergie perturbée du 1° ordre 
est la somme des énergies dues aux différentes 
perturbations. 

Très souvent, les £Ÿ sont nulles : a) soit pour 
des raisons de parité : lorsque A anticommute 
avec l’opérateur parité x et que 4% a une parité 
bien définie ; b) soit pour des raisons de réalité : 
lorsque H{ est un opérateur imaginaire pur et 4 
une fonction réelle pure. Cette dernière condition 
est fréquemment réalisée, car la quasi-totalité des 
molécules possèdent un état fondamental non dégé- 
généré qui peut être décrit par une fonction réelle. 


B) ADDITIVITÉ PAR RAPPORT AUX PARTICULES 
COMPOSANTES. — Nous allons particulariser la fonc- 
tion d'onde 4° en l’écrivant sous la forme d’un 
produit de fonctions que nous appellerons des 
«fonctions partielles », soit : 


ñn 
go = D po. (16-1) 


j=1 
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Les différents 4{” sont supposées normalisées ; 
elles peuvent dépendre des variables dynamiques 
de plusieurs particules, noyaux ou électrons, et 
leur choix (ainsi que leur nombre n) dépend du 
mode de description choisi pour la molécule. Par 
exemple, une fonction partielle décrira le mou- 
vement des noyaux et une autre celui des électrons, 
si l’on a séparé ces deux mouvements à l’aide du 
théorème de Born et Oppenheimer. La fonction 
partielle électronique pourra elle-même être écrite 
de bien des façons : si l’on utilise la méthode des 
orbitales moléculaires, on prendra un produit (que 
nous supposons d’abord non antisymétrisé) de 
fonctions partielles monoélectriques ; si l’on utilise 
la méthode des paires («valence bond method »), 
on prendra un produit de fonctions partielles bi- 
électroniques 4% décrivant une liaison entre deux 
atomes R ets. 

Nous supposerons également que H® et H(? 
peuvent être mis sous la forme d’une somme 
«d’opérateurs partiels », soit 


n 

(1) 

HO = Y Hf 
j=1 


Ht2= ù 4®, 
fl 


(16-2) 
(16-3) 


Les opérateurs A et H® agissent uniquement 
sur les fonctions partielles de même indice 7. Cette 
décomposition est au moins valable pour les per- 
turbations extérieures au système de particules 
considéré. 

A l’aide des définitions précédentes, on obtient 

n 
EDS IIA MINES ER (17-1) 
1 


ei 
avec 


EP = < 0ÿ|A® 0; >. (17-2) 

Le ket [0j > est évidemment le « ket partiel » 
associé à la fonction partielle 4°. Il ne faut pas 
confondre les sommations sur 7 qui portent sur 
différents groupes de particules, et celles sur y qui 
portent sur différents opérateurs. 

(17-1) indique que l’énergie perturbée du premier 
ordre est la somme des énergies perturbées relatives 
aux divers groupes de particules. Si la fonction 
partielle 4° décrit l’ensemble de toutes les parti- 
cules identiques du système, par exemple tous les 
électrons ou tous les noyaux semblables, on trouve 
facilement que le terme Æf° correspondant est 
additif par rapport à tous les corpuscules de l’en- 
semble 7, et conserve la même valeur que la fonc- 
tion soit symétrisée ou non dans le cas de bosons, 
antisymétrisée ou non dans le cas de fermions ; 
en particulier, pour des électrons décrits par la 
méthode des orbitales moléculaires, on obtient le 
même résultat avec un produit antisymétrisé (AP) 
qu'avec un produit simple (SP). 
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IV. Caleul variationnel de l’énergie perturbée du 
second ordre. — À) EXPRESSIONS GÉNÉRALES. — 
L’expression (10) de E® est une fonctionnelle dé- 
pendant de la fonction inconnue 4% (ou du ket 
inconnu |{ >) et nous devons chercher le minimum 
de E®. (Dans ce processus, le dernier terme 
< 0[H®10 > n’interviendra pas.) Utilisons la 
représentation g de la mécanique ondulatoire. En 
explicitant les intégrales de (10), nous obtenons 


EW = fun (HO) — EU) go) 
L pO* (HE — ED) YO E pO* (HD = EU] Ya] dr 


En imposant des variations arbitraires 341 et 
SY®* aux fonctions D et 4%, nous obtenons 


SE = fév (HU — Et9) qth 
+ dE (AC — Et) 54 + S pH (HO — ElH) (0) 
+ QOX (A0) — Et) S 4] dr 


ou en tenant compte de l’hermiticité de H® et 
H® 


SE?) = fByr* (HO — ES) 


+ 8 QD (H(OX — E(0) QN% + S LN* (HN) — EU) g(0) 
+ à pt (HUN*X —— Et1) W(0)*] dr. (18) 


Nous voulons avoir SE® — (0 quelles que soient 
les variations 94 et 54%, ce qui impose à la 
fonction 4% les deux conditions suivantes, équa- 
tions d’Euler du problème variationnel considéré, 
soit 

(HO) — Et) g9 + [HO — EU) Yo = 0 

(HO% — EU) Q* Æ (HOD* — EU)) YO* = 0. 


Ces deux équations sont complexes conjuguées 
l’une de l’autre. Dans ce qui suit nous supposerons 
que l’on a toujours H\%* — HO c’est-à-dire que 
l’hamiltonien non perturbé est un opérateur réel 
pur. Des équations (19) nous pouvons déduire 
immédiatement une relation que nous utiliserons 
plus loin ; en multipliant (19-1) par 4%, (19-2) 
par d® et en intégrant, on obtient 


< 4[H(9 — Et) > + < 1]H0) — E00 > 

< 4/H(0 — Et] > + < 0[H0) — EU > = 0 
d’où l’on déduit 

< 1140 — EN > 


I 
(=) 


= < O[HH— EM > (20) 


ce qui prouve que ces éléments matriciels sont 
réels purs. 

Dans les problèmes moléculaires, on pose fré- 
quemment : 


LOUE) = f(æ) WE) (21) 


NO 


ou | 

M > = f(æ)0 > (21”) 
en désignant par æ l’ensemble des coordonnées 
(d'espace et de spin) de tous les corpuscules de la 
molécule. Si l’hamiltonien non perturbé H® ne 
fait pas intervenir les spins, nous aurons simple- 
ment 


HO Qu = ES A — MU) + VU 
h? 
nr 
+ 2 gradz ie g(0)] + Pf go 
h? 
= fHC® 0) —È HAE [4( Az f 


+ 2 gradz f.gradz (0) 


et les équations (19) deviennent 


2 
(HG) — EU) go) — 2 2e se g(® Az f 


<È es Ge grade VO.gradx f = 0 (2-1) 
(HO De EU) cd 0)* Éd AV f* 
h? 
—5 3} grade UN grade j* — 0. (2222) 


Si Y® est solution rigoureuse des équations (19), 
nous pouvons mettre l’énergie perturbée du second 
ordre sous les différentes formes en utilisant 
soit (18) et (20), soit (19-1), soit (19-1) et (21). 


EU = < 1[H0 — EON > +2 < 1|H00 > 


+ < 0[H(3]0 > (23-1) 
= — < 4[H0 — EU]1 > + < 0[H0 >  (23-2) 
= < 14H00 > + < 0[HEI0 > (23-3) 

fYt0]H]0 > + < O]H0 >. (23-4) 


Si l’on utilise seulement une fonction appro- 
chée, comportant un certain nombre de paramètres 
ajustables, il faut utiliser l’expression (23-1) et la 
minimiser par rapport à ces paramètres. 


IV-B) CARACTÈRE RÉEL OU IMAGINAIRE DE Q, 
CAS PARTICULIERS FRÉQUENTS. — Nous allons uti- 
liser le fait obtenu-en (20) que le terme matriciel 


< 4[HG — E00 > [= < 1]H00 > en vertu de (3’)] 
est réel pur, et étudier particulièrement le cas très 
fréquent où ®® est réelle pure. Mettons tout 


d’abord en évidence les parties réelles et imaginai- 
res par les décompositions 


C0) = ton + à Quoi) 
QD = qÜr + à pub), 


(241) 
(24-2) 


No1 


L'élément matriciel en question peut alors être 
écrit 
< 11H60 > = < 147] HU [Or > + < 1i[HO0i > 


+ i< 4rlHO0i > — < 4 HO [Or >]. (25) 

Les termes du 2° membre sont des intégrales 
calculées sur l’espace de configuration de la molé- 
cule et portant sur les parties réelles ou imaginaires 
des fonctions envisagées. Distinguons maintenant 
deux cas suivant que A est réel pur ou imagi- 
naire pur. On couvre ainsi toutes les possibilités 
car un opérateur complexe peut toujours être con- 
sidéré comme la somme d’un opérateur réel pur 
et d’un opérateur imaginaire pur. 


PH RÉEL PUR. — H = HO% 

En utilisant les décompositions (24), nous trans- 
formons les équations (19) par addition et sous- 
traction en 


(HG) — Et) (On E (H(0) — E(0)) dur) — 
(H0) = Et) (0)  (H(0) — E(0) LU) = 0, 


St YO est réelle pure, le premier terme de(26-2) 
est nul et l’on prendra 44 — 0, 44 devenant 
réelle pure. 

On peut vérifier que ce choix est nécessaire 
lorsque l’état non perturbé n’est pas dégénéré ; 
d’autre part, il est le plus commode lorsque l’état 
est dégénéré et que néanmoins les diverses fonc- 
tions 4° adaptées à la perturbation sont réelles, 
car la valeur de << 1[H%10 = est indépendante 
de la forme de d‘® compte tenu que celle-ci ne 
peut être qu’une combinaison linéaire des 4°. 

D’autre part on a souvent E® — 0, éventuel- 
lement par un choix convenable de l’origine des 
coordonnées, de sorte qu’il reste uniquement l’équa- 
tion portant sur des quantités réelles 


HU g(0) + (H(0 — EC) 4 = 0 (27) 
ou, avec la fonction j introduite en (21) 
h? 


87? mx 


— 2 


k 


HO) go Y g(0) Az f 
k 


Fe gradz ÿ(0).gradg f — 0. (27°) 
mx 


4T? 
. On peut également envisager directement l’ex- 
pression (25) qui se réduit à 


< 41400 > = < 4r|H]0r > — à < 11H07 > 


et l’on remarquera que cet élément matriciel doit 
être réel, ce que l’on peut effectivement obtenir en 
prenant 49 réelle pure. 

Si Y0) est complexe, ce qui est lié à une dégéné- 
rescence, on retrouve, par suite des équations (26), 
le fait que l'expression (25) est réelle pure ; 4 est 
généralement complexe et ne peut se réduire à une 


CALCUL VARIATIONNEL DES PROPRIÉTÉS MOLÉCULAIRES 13 


quantité réelle que si ®% est fonction propre 
de 4% pour la valeur propre EM. 


2. H® IMAGINAIRE PUR. — HO% = — H®, 
En procédant comme ci-dessus, nous pouvons 
mettre les équations (19) sous la forme 


LH) (0) — Et {Or + (H(0) —— EC) din = O0 (28-1) 
HU) on — EU) qi + (HO —— E(0)) QU) = 0. (28-2) 


La première équation est réelle, la seconde imagi- 
naire. 

Si YO est réelle, on sait déjà que E"® = O0 et 
on voit facilement à partir de (27-1) qu’il convient 
ici encore de prendre ÿ47 — O0, c’est-à-dire de 
choisir une fonction de variation imaginaire pure 
satisfaisant à l’équation unique (28-2) qui devient 


HU) (9 + (H(0 — El) gi = 0. (29) 


Si YO est complexe, on voit que la fonction Y® 
l’est également. En vertu des équations (28), l’élé- 
ment matriciel 


< 4/40) — E00 > 


est toujours réel comme il se doit, et l’on peut faire 
des observations analogues à celles indiquées pour 
H® réel. 


IV-C) DÉVELOPPEMENT PAR RAPPORT AUX OPÉ- 
RATEURS DE PERTURBATION. — Nous allons em- 
ployer les expressions .(13-1) et (13-2) des hamil- 
toniens perturbateurs et nous écrirons la fonction 
perturbée du premier ordre sous la forme 


“4 
go = À An WP = Ad (30) 
u=1l 
ou, pour les kets correspondants 

h> = > Alt, ul > = AY (30’) 

et l'énergie perturbée du 2€ ordre sous la forme 

4e 
EG = Yu ED = (À @ À) É®. (31) 
LL, Y 


Les équations d’Euler (19) deviennent 
(H(9) — Et) Ad + AH — EU) 40 — 0 
(H(9) — Et) AQO* + ACHO — E0)) 4% = 0, 


(Les À, peuvent toujours être considérés comme 
réels.) Ces équations doivent être vérifiées quels 
que soient les paramètres indépendants À, et cha- 
cune se sépare en un système de p équations indé- 
pendantes déterminant les p fonctions indépen- 


dantes 4, soit 
(HO FAO) EEE CH EE) LE (32-1) 


Le Fo) LOS (HP EE) NES 066270) 
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La formule (10) devient 
E® = (AQ À).< 11@ (HO — E1 > 
+ (1Q À) [< 1[Q H00 > 


— 

+ < 0[HO QI > + < 0j 0 >] (33) 

où nous avons introduit, par exemple, le tenseur 
< 1] Q (HO — EO)|1 > 
de composantes 
< Au[H(0 — EN > 
et tenu . de la relation générale 
(Au) (Av = (2 je uu) (È à v) 
== 2 Âu À Uyu UV, —= (À ee À) (u ?). 


En comparant (31) et (33) on voit immédia- 


tement que 
_ 
E® = < 11Q (HO — E%)|1 > + < 1|@ HO > 


> 
L < oJHO @IL > + < off >. (34) 


=> 


> 
En remarquant que le tenseur Æ® est symé- 
trique (E%, — E?) on peut écrire ses composantes 


Es = <'IHO — EU v > + <Au/AV10 > 


AH IOE EE CNAB] 

Si les fonctions 4® et 4 sont solutions des 
équations d’Euler « partielles » (32), nous aurons 
plus simplement d’après (23-2) 


ED = — <'Au|AHO — EM v> + < 0[H2]0 >. (34°) 


Certains des £® peuvent être nuls, pour des 
raisons de symétrie géométrique ou de réalité, en 
particulier si 4° est réelle, HŸ? réel tandis que H? 
est imaginaire pur, on sait que D est réelle, 
49 imaginaire, de sorte que le 1% terme de E®) 
est nul. 


IV-D) ADDITIVITÉ PAR RAPPORT AUX PARTICULES 
COMPOSANTES. — Nous considérons maintenant le 
cas où la molécule est décomposée en plusieurs sys- 
tèmes de particules, c’est-à-dire décrite par la 
f. 0. (16). Nous cherchons une fonction perturbée 
d’une forme analogue à (16), soit 


(35-1) 


en prenant pour chaque fonction perturbée par- 
tielle, l'expression 


(0) , ( 2 ,1,(2) Qc 
di = dj” + cd ) + €2 dj" + ... (35-2) 


En portant (35-2) dans (35-1), nous obtenons 
ÿ=IT +eZ IL WW 


2 +e[21 II 2 go +E',I 


2 IL JL ps gi we |: (36) 


En comparant (36) et (2), nous pouvons donc 
poser 


yo (D) 
u-n#(5S)-vw(s#) 


î i 
et il est commode d’écrire, par analogie avec (21) 
gi = = f; go) (38) 
d’où immédiatement 


go = pin) ( DE (37°) 
î , 


Conformément aux hypothèse utilisées dans dif- 
férentes méthodes de la chimie théorique et à la 
forme (16-1) considérée comme valable pour 4, 
nous supposons que l’hamiltonien H® peut être 
pratiquement décomposé en hamiltoniens par- 
tiels AO tels que l’on ait 


HO = Y H (39-1) 


HN po) SE ES go (39-2) 
Il faut bien remarquer que les égalités précé- 
dentes ne peuvent être qu'approximatives. Dans 


ces conditions, l’équation d’Euler (19-1) devient 


EG? — 50) quo f, + (HP 


+2 (mp 


j,k 


SF) gen] 
= EN) YO = 0 
et les termes de la seconde somme sont nuls par 
suite de (39-2). Compte tenu de (16), il reste donc 
D (IL 4) [OR — 2859) 49 + (HP — EP) 91 = 0 
j \kÆi 

Pour que cette équation soit vérifiée, il faut que 
le soient également les équations partielles 

(4 — Ef) ÿ} = 0, . Hi] 


comme on le voit en appliquant la méthode de 
séparation des variables. L'expression (23-4) de- 
vient alors 


EH) y +(H— 


ED = XI f; WOAM lo; > + < 0j|HP 0; >] 
HET< 5 WI0 7 .< HP I0R >. 
2,4 


Les termes de la deuxième somme sont nuls en 
vertu d’une relation analogue à (3”) valable pour 


- les fonctions partielles US et 4 de sorte que 


nous obtenons, de manière générale, une somme 


No1 


d'intégrales portant sur les ensembles partiels de 
particules, soit 


HU=DES (41) 
|] 
avec 
EP = < 1j|Hÿ) [0j > + < 0jÿ}HP 0j >. (1°) 


Supposons maintenant que nous ayons affaire à 
un système de » électrons. À partir d’un produit 
simple (SP) de la forme (35-1), nous obtenons 
l’additivité suivant (41), avec des énergies indivi- 
duelles £® (SP). On sait qu’une meilleure fonction 
est obtenue en tenant compte du principe de Pauli 
et en prenant un produit antisymétrique (AP). On 
verra facilement que la relation d’additivité (41) 
est encore valable (de manière beaucoup plus for- 
melle et moins pratique) avec des E@ (AP) ; mais 
en général on a EF? (AP)  E@ (SP), car l’équa- 
tion d’'Euler totale n’est plus séparable en équa- 
tions d’Euler individuelles, portant chacune sur une 
fonction 4. Ce résultat a déjà été constaté dans 
des calculs sur la polarisabilité électrique [2c] et 
est corroboré par un récent travail d’Allen appli- 
quant la méthode de Hartree-Fock aux systèmes 
perturbés [6]. 


V. Caleul de l’énergie du 3° ordre. — V-A) Ex- 
PRESSIONS GÉNÉRALES. -— Reprenons l’expres- 
sion (11), soit 
EG) = < 1|H00 — Et0f2 > + < 2]H(0 — EN > 

+ < 0[H0 — E02 > + < 2/H0 — EU > 

+ <1[H0 — EON > + <0[H%NZ + <1/49102 


la fonction 49 a été déterminée par le calcul rela- 
- tif au 2€ ordre. Quelle que soit Ÿ?, les termes de 
la première ligne sont nuls en vertu des équa- 
tions (19) et il reste seulement 


EG) = < 1|H0) — Eu > 
+ < 041 > + < 11H10 >. 


(42) 


L'énergie perturbée du 3€ ordre ne dépend donc 
que des fonctions ÿ® et 4, 


V-B) DÉVELOPPEMENT PAR RAPPORT AUX OPÉRA- 
TEURS DE PERTURBATION. — En tenant compte 
des formes (13-1), (13-2) et (30), nous pouvons 
écrire 
EG) = (1Q À@Q À) [< 11Q (HO — EN) Q|1 > 

> > 
> > : 
+ < 0[AD QI1> + < 1,Q H(A|0 >] 


= (1@ A@ ÀA)EM (43) 


en convenant que le symbole E représente un 
tenseur de rang n en même temps qu'il indique 
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l’ordre de la perturbation. Ici E désigne le ten- 
seur de rang 3 : 


EG — < 1,Q (H — El) Q/1 > 
> > dé 
— = 
+ < 0H QI > + <11Q HO > (44) 
de composantes 
Eve = tu|HW + ED} _ 
+ < 0[AP 6 > + <'iulAP)]0 >: (449 


Il convient d’ailleurs de remarquer que ces com- 
posantes n’interviennent dans (43) que par des 
combinaisons symétriques en pu, v et p. 

Certains de ces termes pourront être nuls par 
suite des symétries géométriques ou du caractère 
réel imaginaire des opérateurs et des fonctions ; en 
particulier, si 4‘ est réelle et si 1 ou 3 des in- 
dices u, v, b correspondent à des opérateurs AH, 
imaginaires (x — um, v ou p), donc à des fonc- 
tions 40 imaginaires, la somme des 6 termes Ef?,., 
sera nulle (u', v’, p’ permutation de u, v, p). 


V-C) ADDITIVITÉ PAR RAPPORT AUX PARTICULES 
COMPOSANTES. — En tenant compte de la forme 
(16-1) de Y® et des définitions (16-2) et (16-3), 
(42) devient 


go = D < f, YOlA — HE go > 
j,k,l 


2 


es D L< o[4f, go > à <f, vo] HE 


0 >] 
DINAN EP ENT 0 A 


1 


+ < 4/40; >] 


+ >, [<i> < 0k[Aa = Ex > 
. 


+ < 15105 > < 2 h1r > 
E < oi 2 <AkAPIOK > 


+ < 0jftj > < |A — EP|0r > 


+ '< 1510; > < 0422 — E|0r > 
+ 2 <Ajloj > < OA 2 EDor > < oi >. 
j,k, 
Les 2e et 32 sommes s’annulent soit par suite 
de (17-2), soit par suite des conditions 
TONER 0: 


Donc E% est additive par rapport aux ensembles 
partiels de particules et nous pouvons écrire 


go — ÿ K® (45) 
avec = 
PNR EM 
RAP > + <Aj]4Ÿ 0; >. (459 
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Si nous tenons compte de l’additivité par rapport 
aux Opérateurs, nous pourrons écrire des termes 
de la forme 
Elps = < 1, j|HV3 — ESJI1, 6, 1 > 


LE <oHD 1, p,1> + <iu, HP 40 >. (45’) 


VI. Calcul variationnel de l’énergie du 4€ ordre. 
— VI-A) EXPRESSIONS GÉNÉRALES. —— Partons de 
l’expression (12) 

EG = < 1]H(0 — EU0j3 > + < OH — EU > 
+ < 8[H(00 — EN > + < 8[H0 — EU0 > 
+ < 2]H(0 — H0|2 > + < A[H) — EU2 > 
+ < 0[A) — Et3}2 > + < 2H) — EG > 
+ < 2H) — E|0 > + <1[H — EN > 
__ E®(< 0h > + <1[0 >). 

La première ligne est nulle si ÿ"® vérifie les 
équations d’Euler (19); il en est de même du 
dernier terme par suite de la normalisation. La 


seconde ligne dépend de 4® et les conditions de 
minimum donnent les deux équations conjuguées : 


(HG — EXO) ga* + (HO — EU) go 
+ (HO% — E®) OX — O (46-1) 
(HU) 2 Et) 40) + (HO) — EU) gU) 
+ (H® — EU) go = 0. (46-2) 


Si Y® vérifie ces équations, il reste simplement 
EU) = — <.9|H(0) — E0f2 > + << 1|[H— EN >. 
(47) 


L'énergie perturbée du 4€ ordre ne dépendra que 
des fonctions Y® et pt? 


VI-B) DÉVELOPPEMENT PAR RAPPORT AUX OPÉ- 
RATEURS DE PERTURBATION. — Par une décompo- 
sition analogue à celle utilisée pour 4, nous 
aurons 


ae 
Yo = Ah dés = (À @ À) (48) 
L,Y=1 
ou encore 
— 
> 


R> = SA 2, uv > = (AQA)LZ. (48) 
L,V 


En tenant compte de (13-1), (13-2) et (30), nous 
déduirons de (46) une série d'équations qui déter- 
minent, indépendamment les unes des autres, les 
composantes 4 


EY) 4 pe 


+ (HS — 


LH EU) JP LE (pe 


ER) g9 — 0 (49) 


ou encore puisqu'on peut prendre une fonction 
symétrique en u et v 


Cao po) g@ + : (HD — EN) 


1 
+3 (HP — EN) qu + (His — Eu) YO = 0. (49°) 


Quant à l’énergie £, elle peut, avec cette nota- 
tion être mise sous la forme 


Lä Ê 
EU = (1@ A@ AQA)l—<2|Q(H0— EN) > 
sms 
+ <11@ (40 — EDIQIL à 


= (1Q AQ À1@ À) Et" 


E® est un tenseur du 4€ rang ; 
sont 


(54) 


ses composantes 


Es = — < 2, uv|H(0) — Et02, po > 


+ <iujAŸ — Eflio >. (50) 


Lorsque 4® est réelle, on voit facilement, à 
partir de l’équation d’Euler (49) que 4% doit 
être prise réelle pure si A et HŸ ont le même 
caractère de réalité, et qu’elle doit être imaginaire 
pure dans le cas contraire ; en outre, si un nombre 
impair des indices u, v, p, « correspond à un opéra- 
teur H® imaginaire pur, la somme des 24 termes 
Pos est nulle (u’, v’, p’, o’ permutation de y, 

V, 0, & 

On peut également faire, au sujet des symétries 
géométriques, le même genre de remarque pour E® 
que pour E. 


VI-C) ADDITIVITÉ PAR RAPPORT AUX PARTICULES 
COMPOSANTES. — La fonction d’onde non perturbée 
étant sous la forme (16-1), les opérateurs sous les 
formes (16-2) et (16-3), nous avons défini f par la 
relation (21). De même, en égalant à 4? le terme 
en €? de (36), nous définirons la fonction g telle 
que 


gen = go = ZX IT Yo dé +Z' IT WW 
î ki î, 


k ik 
(51) 
(0 Ut (0 ! dan a ? 
poz yo poz 0 EU (51°) 


Si nous posions, de manière analogue à (38), 
pour chaque groupe 7 de particules 


Ph gr Uf (51”) 
il viendrait 
gta = vo lz g +Z'h he | (52) 
#2 j,k 
ou, en comparant (51”) et (52) 
= 2 +X' ff. (53) 
j 3, 


No1 | 


ee #| 


NET 


La fonction g ne pourrait être considérée comme 
additive. En réalité, le développement qui suit 
montre qu'il est équivalent de prendre la fonction 
additive 


s— De 
(2) = È au 


(53 
(53") 


puisque chaque g; obéit alors à la même équation 
d’'Euler et que l’on aboutit aux mêmes formules 
qu'avec le choix (53) que nous utilisons ci-dessous. 

Remarquons en passant que les conditions de 
normalisation (3), (4) ou (3’), (4’), écrites pour les 
fonctions globales V® et 4? peuvent être égale- 
ment appliquées respectivement aux fonctions QU, 
V2 où YD, E. 

“L équation d’Euler (46-2) peut s’écrire, à l’aide 
de (53) 


Z | — Ef9) qio (ss + Eh) 


+ (AP — EF) qi f, + CA — EP) YO — 0 


ou encore, en séparant les termes pour lesquels : 
et 7 sont égaux 


>> [U4” — Ef9) go ( gi + fi . fe) 


+ (49 — EN) Uo f; + (HP — EF?) je 


+2 (A 20) yo [e +h2 fe) 


+ (AP — EÙ) go hi] = 0. 

Les termes de la seconde ligne sont nuls en vertu 
de (39-2) et (40) supposées vérifiées. En multi- 
pliant scalairement les termes de la première ligne 
par He, y et en tenant compte des conditions de 


an, on obtient l’équation déterminant 
la fonction partielle g, SOIT 


(AY — Ef) 9 gi + (HE — ES) WE fi 

+ (HP — ES?) 9 = 0 (54) 
ou 
(H$9 — Ef9) 9 + (HP — EN) W° 

+ (HP = - Ef?) y 0 0. (54°) 


de forme analogue à (45). 
En tenant compte de (38) et (52), nous pouvons 
mettre £% sous la forme 
n 


EU) — 2 
j,k,l1=1 


[< vo, |A — EP)| go) f, Fe 
—< (8 +f 2 im) lHE— Eÿ] 4 {a+ h Eh)? |. 


En utilisant les conditions de normalisation, 
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l'équation de Schrodinger partielle (39-2) et en opé- 
rant comme pour E®%, on obtient la nouvelle 
expression 


EU) = D [= EEE he <o AMEN 


_. 2 <1jl1j > [< 044 — EP ]0r > 


— < 144 — ED l18 >]. 
En multipliant scalairement l’équation d’Eu- 
ler (40) par Y, l’équation d’Euler (54) par 49, 
et en retranchant les deux expressions obtenues, 
on trouve que : 
< 0j|4P — Elo; > 


— AA = Eh; > = 0 


c’est-à-dire que l’on peut écrire pour E® 


n 
EU) = Y Ef (55) 
j=1 
avec 
Eÿ9 = <1j|4P — Eh; > — <a2j[a9 — Eh; > 
(55') 


L'énergie du 4 ordre, comme celles des ordres 
inférieurs, est additive par rapport aux différents 
groupes 7 de particules. 


VII. Remarques finales. -— Résumons rapi- 
dement les principaux résultats du formalisme gé- 
néral exposé ci-dessus. Développant l’énergie et la 
fonction d’onde du système étudié en série des pa- 
ramètres de perturbation, nous avons obtenu les 
équations d’Euler (19) et (46) qui déterminent les 
fonctions perturbées du 1° et du 2€ ordres. Il 
convient d’ailleurs de remarquer que ces équations 
sont identiques à celles de la théorie des pertur- 
bations [7] ; notre méthode est donc soumise aux 
mêmes réserves quant à la convergence des expres- 
sions obtenues. 

En outre, nous n’avons pas, la plupart du temps, 
mentionné explicitement le fait que 4® peut cor- 
respondre à un état dégénéré ; le problème est 
alors posé de savoir de quelle manière il convient 
de lever la dégénérescence, en utilisant de préfé- 
rence la perturbation totale, ou éventuellement la 
perturbation la plus intense. 

Nous avons également obtenu les expressions 
donnant les énergies perturbées jusqu’au 4€ ordre. 
Les équations et les expressions dont il est question 
peuvent d’ailleurs se rapporter soit à la pertur- 
bation totale, soit aux perturbations partielles su1- 
vant lesquelles on peut décomposer la pertur- 
bation totale ; dans ce dernier cas, les expressions 
des énergies font apparaître, à partir de l’ordre 2, 
des termes « croisés », c’est-à-dire correspondant à 


2 
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des interactions possibles entre les différents opéra- 
teurs perturbateurs. Ceux-ci d’ailleurs peuvent cor- 
respondre à des ordres de grandeur expérimentaux 
différents, que notre méthode abstraite ne met pas 
en évidence ; c’est seulement dans des applications 
précises que l’on devra s’en préoccuper. Il paraît 
particulièrement intéressant de noter que l’énergie 
perturbée du 32 ordre EŸ est calculable à partir 
de ®® et 4, et l’énergie perturbée du 4€ ordre 
E® à partir de 4 et 4? seulement [10]. 

D'autre part, moyennant certaines hypothèses 
simplificatrices sur la f. o. totale, on peut séparer 
les équations d’Euler en équations portant sur des 
groupes de particules (éventuz2llement, sur une 
seule particule) et les énergies obtenues sont addi- 
tives par rapport aux différents groupes de parti- 
cules ainsi distingués. Nous l’avons établi assez lon- 
guement, mais nous pourrions nous contenter de 
remarquer qu'il s’agit là d’une conséquence bien 
connue de la forme sous laquelle, grâce à diverses 
hypothèses, nous avons écrit l’hamiltonien per- 
turbé total, soit 


n 
H = Ÿ H;. 


jæ1 


NS 


Dans ces conditions, une f. o. perturbée peut 
être écrite sous la forme 


ÿ = IT 4, 
3 
avec 
Hj ÿ5 = E; Ÿ 
et 
E =YE,. 


1 

Les divers développements, caractéristiques de la 
méthode des perturbations, peuvent alors être 
introduits à propos de chaque indice j et conduisent 
aux résultats obtenus plus haut. 

Remarquons, pour finir, que les grandeurs cal- 
culées dans cet article, sont des grandeurs micro- 
scopiques, se rapportant à une molécule. Pour obte- 
nir des grandeurs macroscopiques, il est nécessaire 
d'effectuer une moyenne statistique sur tous les 
états de la molécule et sur toutes ses orientations. 
On peut trouver des indications sur cette question 
dans l’ouvrage de Van Vleck [8] ou dans [9]. 


Manuscrit reçu le 21 septembre 1960. 


BIBLIOGRAPHIE 


[1] Gans (R.) et MrowxA (B.), Kônigsberger Gelehrte 
Gesellschaft Nachr., 1935, 12, 1. 
[2] a) Szater (J. C.) et Kirkwoop (J. G.), Phys. Reo., 
1931, 37, 682. 
b) Guy (J.) et HarranD (M.), C. R. Acad. Se., 1952, 
234, 616 et 716. 
c) ABBOTT (J. A.) et BorTon (H. C.), Proc. Roy. Soc., 
1953, À 216, 477. 
[(8]-a) Tizzieu (J.) et Guy (J.), C. R. Acad. Sc., 1954, 289, 
- 1203 et 1283 ; J. Chem. Phys., 1956, 24, 1117. 
Tizzieu (J.), Ann. Physique, 1957, 2, 471 et 631. 
Das (T. P.) et Bersoun (R.), Phys. Rev., 1956, 104, 
849 ; Phys. Rev., 1959, 115, 897. 
b\.STEPHEN (M. J.), Proc. Roy. Soc., A, 1957, 242, 
264 et 274. 
[5] a) Courson (G. A.), Maccozz (A.) et SuTron (L. E.), 
Trans. Farad. Soc., 1952, 48, 106. 
b) BucxixcnaAm (A. D.) et Popze (J. AÀ.), Proa. Phys. 
Soc., 1955, À 68, 905. 


) 
[4} a) 


c) BuckiNGHAM (A. D.) et Porre (J. A.), Proc. Phys. 
Soc., 1956, B 69, 1133. 
d) BucxinGxaAm (A. D.), Proc. Phys. Soc., 1957, B TO, 
758. 
e) BucxincHAm (A. D.), Quart. Reo., 1959, 18, 183. 
[6] ALLEN (L. C.), Phys Rev., 1960, 118, 167. 
[7] Cf. par exemple KRAMERS (H. A.), Quantum Mechanics 
(North Holland Publ. Co), 1956, 188 sq. ou Messtax 
(A.), Mécanique quantique (Dunod), 1956, tome II, 
p. 586 sq. 
[8] Van VL£Ecx (J. H.), The theory of electric and magnetic 
susceptibilities (Oxford, Univ. Press), 1932. 
[9] Axprews (A. L.) et BuckiNGxaAm (A. D.), Mol. Phys., 
1960,3, 183. 
[10] Pour une généralisation de ce résultat, cf. Duponr- 
BourpeLer (K.), Tizzieu (J.) et Guy (S.), J. Phy- 
sique Rad, 1960, 21, 776. 


4 


É 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


TOME 22, JANVIER 1961, PAGE 19. 


RÉSONANCE FERROMAGNÉTIQUE DANS UN ALLIAGE Fe-Ni 
SOUS FORME DE COUCHES MINCES 


Par A. VAN ITTERBEEK, J. SMITS, G. FORREZ et J. WITTERS, 


Institut voor Lage Temperaturen en Technische Fysica, Leuven (Belgie). 


Résumé. — Utilisant un spectromètre à ondes centimétriques nous avons étudié les raies de 
résonance de plusieurs films obtenus par évaporation du permalloy C, sous vide, sur un support 
en verre. Le facteur g et l'induction sont calculés à partir des grandeurs mesurées. On a observé 
l'influence d’un traitement à haute température sur les facteurs suivants : la forme de la raie, la 
valeur du champ magnétique de résonance et la valeur de l’induction. 


Abstract. — Using a microwave spectrometer we have studied the ferromagnetic resonance of 


films, obtained byvacuum evaporation Permalloy C, on a glass layer. 


tion are calculated from the measurements. 


The g-factor and the induc- 


The influence of annealing is studied on the follo- 


wing factors : the lineshape, the value of the resonance field and the value of the induction. 


{. Introduction. — Poursuivant nos recherches 
sur la résonance ferromagnétique d: couches minces 
dont les premiers résultats furent publiés récem- 
ment dans ce même journal [1] nous avons perfec- 
tionné notre appareillage de sorte que la déter- 
mination de la valeur du champ de résonance est 
plus précise. Ceci nous a permis de calculer les 
valeurs du facteur g et de l’induction pour les diffé- 
rentes couches. Dans nos mesures antérieures nous 
avions observé l’apparition de raies doubles de réso- 
nance. Dans le présent travail nous en avons étudié 
les causes. 

Les films sont obtenus par évaporation du per- 
malloy GC, mis à notre disposition par le « Vacuum 
Schmelze Aktiengeselschaft » de Hanau (Alle- 
magne). Après une première mesure, les films sont 
soumis à un échauffement de 300 °C à 500 0C 
pendant des temps différents. Ceci se fait également 
sous vide. Les mesures taites après ce traitement 
nous donnent une indication importante concer- 
nant l’influence de cet échauffement sur les pro- 
priétés magnétiques de l’alliage. 


2. Dispositif expérimental. — La figure 1 repré- 
sente le schéma du dispositif expérimental. Le 
klystron fournit une fréquence de 9 230 MHz, sta- 
bilisée par un stabilisateur IF de Pound dont l’oscil- 
lateur à 9,5 MHz est piloté par un crystal ther- 
mostatisé. Après détection le signal de basse fré- 
quence (64 Hz) est appliqué à l’entrée d’un ampli- 
ficateur sélectif de Bruel et Kjaer, qui précède un 
détecteur de phase. 

Le champ magnétique constant (41) est livré 
par un électro-aimant Oerikon à basse résistance, 
alimenté par une dynamo, dont l’excitation pro- 
vient d’une deuxième dynamo. Celle-c1 est placée 
sur le même axe et nous en réglons le courant 
d’excitation au moyen d’un hélipot de 200 ohms 
et A0 tours. 


La disposition est telle que Æ, est perpendi- 
culaire au champ magnétique HF. Le champ alter- 
natif (k cos Qt) de modulation est livré par un 
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alternateur de 64 Hz. Les films sont toujours 
placés parallèlement au champ magnétique HF, 
tandis que l’orientation par rapport au champ ma- 
enétique constant est variable. L'erreur relative 
des mesures est de 3 %. Elle provient essentiel- 
lement de l’incertitude de la relation entre le cou- 
rant et le champ magnétique produit par l’électro- 
aimant. 


3. Les échantillons. — Les couches métalliques 
sont déposées dans un appareil construit au labo- 
ratoire et qui fut décrit antérieurement [2]. La 
pompe à diffusion fournit un vide de 1076 mm Hg ; 
pendant l’évaporation, le vide est supérieur à 
1075 mm Hg. Une lamelle en verre sert de support. 
Elle a les dimensions suivantes : 15 X 3 X 0,6 mm 
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L’épaisseur des films est déterminée par la méthode 
de Tolansky. Elle varie entre 1 000 et 2 160 À, c2 
qui correspond donc à une quantité de métal de 40 
à 85 microgrammes. 


Les résultats. — 1. VALEUR DU CHAMP DE RÉSO- 
NANCE. — Le tableau I donne les résultats expéri- 
mentaux quant au champ magnétique de réso- 
nance dans deux orientations du film : 77, est le 
champ de résonance pour la position parallèle au 
champ magnétique constant, 4, pour la position 
perpendiculaire. 


TABLEAU I 

Fizm N° H; (Œ) H, (Œ) us 
il 4 042 10 923 2 160 
2 TAROT 11 108 2 160 
3 1 018 11 065 2 160 
& 1 068 10 915 . 2 160 
5 1 054 10 865 2 160 
6 998 10 580 — (*) 
7 1 004 10 686 — (*) 
8 999 10 789 — (*) 
9 1 076 nes) 1 030 
10 1 014 _ 1 000 
11 1 006 = 1 000 
12 1 023 = 1 000 
13 1 073 => 1 030 
14 1 109 ee 1 030 


(*) Ces films furent oxydés pendant le traitement à 
haute température et il était donc devenu impossible d’en 
déterminer l’épaisseur. Une estimation basée sur la quantité 
de métal évaporé a donné : 1 000 À. 

(**) La signal est trop faible pour déterminer la valeur 
du champ de résonance. 


Les valeurs du champ magnétique de résonance 
sont données par les deux conditions de réso- 
nance [3] correspondant aux deux orientations 
(4 est suivant l’axe z, qui correspond à l’orien- 
tation du champ magnétique constant). 


o = YVHi(Hi + 47 Mi) (1) 
© = y(H;—4rMi) (2) 


wf27r — v: fréquence correspondante à une 
absorption maximale ; 

MW, : magnétisation maximum dans la direction 
du champ magnétique constant Æ, (resp. H;). 

La dispersion des résultats dépasse la précision 
des mesures et doit donc être attribuée à des diffé- 
rences dans la structure des films. 

Dans le cas où il n’y a pas de raison spéciale qui 
donne lieu à une déformation de la raie, sa largeur 
mesurée entre les points d’inflection est de 50 Œ. 


Pour un film de fer pur, que nous avons mesuré 
pour comparaison cette largeur est de 350 Œ. 


2. CALCUL DE g ET DE 4m M,. — A partir des 
deux valeurs du champ magnétique A, et ,, on 
peut calculer g et 4x M, au moyen des formules 
suivantes : 

H; +2H; —V5H? + 4H, H;, 


4TMs = ; 3) 


= of(Hy— 4x Ms) — gef2me (a) 


g : représente comme d’habitude le facteur de 
Landé ; 

y : rapport gyromagnétique. 

On obtient ces formules en éliminant « entre 
les formules (1) et (2), sous condition que la magné- 
tisation M, atteint, dans les deux cas, sa valeur de 
saturation A, ce qui est certainement vrai. Pour 
le calcul nous avons pris les valeurs des constantes 
physiques suivantes [4] : 

Le magnéton de Bohr : 


ñ 
B — de — 0,9271 x 10—% erg/Œrsted 


ñ — 1,0542 X 10—?7 erg.seconde 
de sorte que 
eme — 0,879 A0. 
La fréquence est : 
9 227 Mhz 


donc : 
o — 57 995 X 106 s—1, 


Le résultat de ces calculs est donné dans le 
tableau II. 


TABLEAU II 
Firm n° 4m Ms (Œ) £ 

1 7 875 2,16 
2 0 077 2,18 
3 8 030 2,17 
4 7-832 2,14 
5 7 805 2,15 
6 7 644 2,24 
7 7 725 2,28 

7 820 2,22 


La valeur de g correspond aux valeurs données 
dans la littérature (en moyenne g varie entre 2,07 
et 2,20) [4], ce qui confirme la fidélité de nos 
mesures en tenant compte des limites des erreurs 
d’observation. Si l’on compare les valeurs de 47 M, 
ou B,; aux données disponibles pour les alliages 
Fe-Ni (+ 8 700 — -- 10 800) [5], on remarque des 


hi te deg = sl Le - 


No 


différences. Ces différences doivent être attribuées 
à des défauts de structure et à des tensions dans 
les films. Ce sont ces facteurs qui, entre autres, 
ont une influence sur B,. 


3. INFLUENCE DU TRAITEMENT A HAUTE TEMPÉ- 
RATURE. — Le traitement à haute température que 
l’on opère pour obtenir les alliages dits « per- 
malloy » a nettement pour but d’annihiler les 
défauts et d’augmenter ainsi la valeur de B,. Ceci 
est confirmé expérimentalement, comme montre le 
tableau ITT. 


TABLEAU III 


Fizm N° TRAITEMENT AH, (Œ) H,(ŒE) 4rM (Œ) g 
1 7 h à 500 °C 798 13 882.10 835 116 
2 7 h à 370 °C 868 12 763 927182 2,17 
3 7 h à 500 °C 7994 43 666 10-642 2,18 
4 7 h à 500 °C 817 ASS NT 076 "17 
5 14 h à 500 °C 805 AGO OL GPS 2 TS 
6 2 h à 350 °C 966 11 286 So 20 2 20 


On remarque que les valeurs ainsi obtenues 
pour 4x4, se rapprochent d’autant plus de celles 
du permalloy que la température est plus haute 
et que la durée du traitement est plus longue. 


arbitraires 


1000 1500 


2000 Œ 

Fi1G. 2. — Spectre d’absorption d'ondes centimétriques 
dans un film de permalloy sur un support en verre ; 
épaisseur du film : 2 160 À). 
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4. ANISOTROPIE. — Quand l’état de tension n’est 
pas le même à travers tout le film, ceci peut 
donner lieu à des conditions de résonance diffé- 
rentes pour la surface libre du film et pour celle 
qui est en contact direct avec le support. Une 
tension peut apparaître par suite de la différence 
entre les coefficients de dilatation du métal évaporé 
et du support en verre. D’autres chercheurs ont 
également trouvé une différence entre la condition 
de résonance pour la surface libre du film et pour 
celle qui est attachée au support. Ainsi pour du 
nickel, Griffiths [6] trouve une différence du champ 
magnétique de résonance de 400 Œ. Dans nos 
expériences les deux régions se trouvent en même 
temps dans le champ magnétique HF et nous 
observons simultanément deux raies de résonance. 
Dans la figure 2 ce phénomène est illustré pour le 
film n° 2. Le maximum le plus élevé correspond, 
comme montre la valeur du champ magnétique 
(1 011 Œ), à la résonance principale, tandis que 
l’autre raie de résonance est causée dans la région 
en contact direct avec le support. La distance 
entre les deux maxima est de 324 Œ. 

Pour contrôler cette hypothèse nous avons fait 


= 
1509 


500 1009 Œ 


F1G. 3. — Raie de résonance après traitement de 7 h 
à une température de 370 °C, 
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-_ quelques expériences sur le traitement à haute tem- 
pérature. Après un traitement de 7 h à 370 ©C, 
DRE suivi d’un refroidissement très lent (fournace 
arbitraires 5 ; : : 
cooling) la raie de résonance devient celle de la 
figure 3. On remarque la diminution de la raie 
supplémentaire et le déplacement vers les champs 
plus petits. Après un deuxième traitement de 7 h 
à 500 °C avec le même procédé de refroidissement, 
la courbe a perdu son deuxième maximum (fig. 4). 
À la fin de cet article nous tenons à remercier 
PUnion Minière du Haut-Katanga, ainsi que le 
Centre Interunivaersitaire de Recherches Nucléair2s 
pour leur aide financière au cours de ces recherches. 
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D ce Mn ance Manuscrit reçu le 14 octobre 1960. 
après un deuxième traitement de 7 h à 500 °C. 
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SUR L'ÉNERGIE D'EXCITATION DU PREMIER NIVEAU 2+ DES NOYAUX PAIR-PAIR 


Par C. YTHIER et R. VAN LIESHOUT, 


Instituut voor Kernphysisch Onderzoek, Amsterdam. 


Résumé, — La présence d’une structure fine dans la variation de l’énergie d’excitation du pre- 
mier niveau 2+- des noyaux pair-pair en fonction du nombre de neutrons dans la région 
90 < N < 116, suggérée par Chupp e. a. [1], semble confirmée. Une structure analogue semble pré- 


sente dans la région 30 < N < 46. 


Abstract. — The variation of the energy of the first excited state in even-even nuclei, as a 


function of the neutron number, shows a fine structure. 


The existence of this effect has already 


been reported by Chupp et al. for the region 96 < N < 126 and it is confirmed by more recent 
data. An analogous effect seems to be present in the region 30 < AN < 6. 


L'énergie d’excitation du premier niveau 2 + des 
noyaux pair-pair est connue dans un très grand 
nombre de cas avec une précision de quelques 
pour-cents. Dans la région où le nombre de neutrons 
présents dans les noyaux varie entre 90 et 116, 
c’est-à-dire dans la région de déformation perma- 
nente, cette précision est même beaucoup plus 
grande dans de nombreux cas, là où le rayon- 
nement de désexcitation a été mesuré à l’aide d’un 
spectrographe à cristal courbé. Chupp e. a. [1] ont 
remarqué que la courbe représentant pour chaque 
élément la variation de Æ(2 +), en fonction de NW, 
bien qu'ayant l’allure générale prédite par la théorie 
de Bohr et Mottelson [2], présente une structure 
fine : par exemple IS0Hf a presque la même éner- 
gie E(2+), que !8Hf, comme si le moment d’iner- 
tie du noyau ne diminuait pas lorsque NV passe de 
106 à 108 ; il semble d’autre part exister un mini- 
mum dans la variation de £ aussi bien à N = 98 
(164Dy) qu’à N = 104 (174Yb). 

Afin de mettre mieux en relief les particularités 
que l’on pourrait qualifier de structure fine, 1l nous 
a paru intéressant d’examiner à nouveau, en tenant 
compte des données expérimentales obtenues 
depuis l’observation de Chupp e. a. [1], la varia- 
tion de E(2+), dans la région 90 < N < 116. La 
figure 1 montre qu’en dépit de l’erreur probable 
assez importante de ces déterminations spectros- 
copiques récentes [3] il semble exister, à côté des 
discontinuités signalées plus haut, d’autres dis- 
continuités encore : à N — 108 pour le tungstène, 
comme le montre la position inattendue de I$#0W ; 
à N — 98 pour l’erbium, comme il résulte de la 
valeur trouvée récemment pour 154Er et si l’identi- 
fication de 17Er est correcte [1]. L’anomalie du 
hafnium à V — 106 est soulignée par la position 
de 17%4Hf le minimum dans les courbes du dys- 
prosium et de l’ytterbium (si l’identification de 
176Yb est correcte [1]) est soulignée par la valeur 
trouvée pour 160Dy et pour 168Yb et 170Yh. 

Chupp e. à. [4] ont d’autre part appelé l’atten- 


tion sur la linéarité remarquable de la variation 
de E(2+), pour le tungstène entre N — 108 et 
N — 112. Comme le montre la figure 1, cette varia- 
tion semble être également linéaire, dans un certain 
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FiG. 1. — Variation de l’énergie d’excitation. du premier 
niveau 2+ des noyaux pair-pair en fonction du nombre 
de neutrons entre N — 90 et N —116. Les valeurs de 
l’énergie portées sur la figure sont prises dans les « Nu- 
clear Data Sheets » (3). — Erratum : au centre et à 
droite de la figure, il faut un trait plein entre :8Hf 
et 1H. 
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intervalle, pour le dysprosium (N — 94 à 98), pour 
lPerbium (N = 98 à 102) si l'identification de 17Er 
est correcte [1], pour l’ytterbium (W — 100 à 104) 
si l'identification de 176Yb est correcte [1] (l’incer- 
titude sur la valeur de l’énergie dans le cas de 
168Yh ne permet pas de dire si cette linéarité 
s'étend à N = 98),et probablement pour le hafnium 
(N = 102 à 106). Pour les éléments W, Hf, Yb 
et Er, cette linéarité est brusquement interrompue 
par les discontinuités relevées ci-dessus. Pour 
mieux mettre en évidence cette situation, il est 
commode de comparer, pour chaque élément, la 
moyenne des énergies trouvées pour ces valeurs 
de N, par exemple 108 et 112 (W), à l’énergie 
trouvée pour la valeur intermédiaire de A, ici 
N = 110. Le tableau I montre le résultat d’une 
telle comparaison. 


e 
——— 


+ E (MeV) 


Nombre de neutrons N ——+ 


FiG. 2. — Variation de E(2+), pour les noyaux pair-pair 
de la région 28 < N < 50 (cf. légende fig. 1). 


En dépit du caractère incomplet et, dans quel- 
ques cas, encore imprécis des données sur £(2+), 
dans cette région, son examen semble indiquer que 
la variation de Æ n’a pas l’allure continue qu’on 
attendrait d’après la théorie de Bohr et Mottelson. 

Afin de savoir si la présence d’une telle structure 
fine est caractéristique des noyaux à déformation 
permanente, l’analyse a été étendue aux autres 
domaines de MW, à l'exception des noyaux légers, 
car on en connaît trop peu de même Z et de N 
différents, et des noyaux pour lesquels N corres- 
pond à la fermeture d’une couche de neutrons [5] 
ou aux limites de la région de déformation perma- 
nente [6]. (La variation en fonction de Z offre peu 


d'intérêt car on connaît trop peu de noyaux de 
même N et de Z différents. Une comparaison de la 
position des premiers et seconds niveaux 2 + est 
difficile en raison de la dégénérescence des seconds 
niveaux excités [6, 31.) 

Dans la région 30 < N < 46, la variation 
de Æ(2 +), pour un élément donné semble bien être 
linéaire entre certaines valeurs de NV, comme le 
montre la figure 2. Le tableau IT permet de vérifier 
cette linéarité : l’écart ne diffère de zéro de façon 
significative que dans deux des cas étudiés et, 
même dans ces deux cas, sa valeur relative est 
inférieure à 2 %. Pour tous les éléments examinés 
ici, la variation de Æ a même allure : l’intervalle où 
elle est linéaire est brusquement interrompu par 
une discontinuité très marquée lorsque N atteint 
la valeur 34 (Ni), 38 (Zn), 42 (Ge et Se) et 46 (Kr) ; 
entre deux valeurs consécutives de cette suite, la 
variation de Æ a le même sens pour tous les élé- 
ments ; ce sens diffère d’un intervalle au suivant. 

La linéarité n’est pas rompue pour NV — 40, 
nombre qui doit correspondre à la clôture d’une 
sous-couche de neutrons ; on observe seulement un 
relèvement d’ensemble de la courbe qui culmine à 
N — 38. Ce fait mérite d’être comparé aux obser- 
vations de Way, Kundu, MecGinnis et Van 
Lieshout [7], qui ont signalé l’importance du 
nombre N — 38 dans la classification des spins des 
états fondamentaux des noyaux. La position très 
élevée des noyaux de nickel dans la figure 2 est en 
relation avec la valeur de Z (Z = 28). On voit que 
l’augmentation de Æ par rapport aux éléments 
voisins, par exemple le chrome (Z — 26), est du 
même ordre que celui observé pour N = 28. 

Dans la région 50 < N < 82, l'influence de la 
fermeture d’une couche de neutrons à N — 82 est 
déjà ressentie pour V — 72 ; en outre, non seule- 
ment l’étain, mais aussi le cadmium et le tellure 
manifestent l'influence de la fermeture d’une couche 
de protons à Z — 50 [5] ; dans ces domaines l’éner- 
gie d’excitation E(2 +), est beaucoup plus élevée : 
aussi serait-il utile, pour y relever une éventuelle 
structure fine, de connaître la valeur de Æ avec 
une très haute précision. 

La figure 3 montre que les seuls détails remar- 
quables dans cette région sont d’une part une dis- 
continuité marquée (à M = 56 pour Z = 42, à 
N:=: 60 pour. Zi 40, AN" 62spour 7 =; 
à N — 64 pour Z — 50), d’autre part une varia- 
tion remarquablement linéaire dans certains inter- 
valles (NW = 60 à 64 pour Z = 46, N — 62 à 66 
pour Z =48;/N= 68 à 72 pour Z —=.50). 

On ne relève aucune ressemblance entre la varia- 
tion de Æ dans cette région et celle observée dans 
la région 30 < N < 46. Il n’en est pas de même 
si l’on examine la variation de B, au lieu de celle 
de Æ (B, est le paramètre de masse du mouvement 
de vibration). Cette variation a été examinée dans 
le cas du cadmium, de Pétain et du tellure par 
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Alkhazov e. a. [8] ; la figure 6 de cette référence 
suggère l’existence d’une structure fine anslogue 
à celle signalée dans la région 30 < N < 46. 


+ 


Nombre de neutrons N ——+ 


F1G. 3. — Variation de Æ(2+), pour les noyaux pair-pair 
de la région 50 < N < 82 (cf. légende fig. 1). 


* 
* * 


En conclusion, les observations faites par 
Chupp e. à. [1, 4] peuvent être étendues à quelques 
autres cas dans la région 90 < N < 116. Il existe 


Nos 


d’autre part une structure remarquable dans la 
région 30 < N < 46. Enfin quelqu:s particularités 
peuvent être signalées dans la région 50 < N < 82. 

S'il n’est pas possible à partir de ces observations 
d'élaborer une règle générale, le rôle joué par le 
nombre de neutrons A dans la variation de E(2 +), 
ressort néanmoins, par exemple dans la tendance 
manifestée par cette variation à être linéaire dans 
des intervalles dont la largeur est le plus souvent 
ANNEE 

Il apparaît également que la structure fine revêt 
des aspects différents selon le domaine de N consi- 
déré. 

Aucun de ces aspects ne paraît avoir été prévu 
par la théorie. Non seulement aucune théorie n’a 
conféré une importance particulière à des paires de 
paires de neutrons, mais encore on ne voit pas, 
par exemple, comment expliquer dans le modèle 
en couches l’apparition de discontinuités aussi bien 
à N — 34, 42 et 46 qu'à N — 38 pour les élé- 
ments Ni, Zn, Ge, Se et Kr. Si un maximum à 
N — 64 dans la variation de E(2+), pour l’étain 
n’est pas surprenant, puisque l’importance de ce 
nombre a été signalée à propos d’une classification 
des spins des états fondamentaux des noyaux [7], 
on est étonné d’observer un maximum analogue 
pour le cadmium à N = 62. Enfin les particularités 
observées dans la région des noyaux à déformation 
permanente contraste avec la variation régulière 
que laisse prévoir le modèle de Bohr et Mottelson. 

Le présent travail, exécuté dans le cadre du pro- 
gramme de la Fondation pour les recherches fonda- 
mentales sur la matière (F. O. M.), a bénéficié de 
l’aide de l’Organisation néerlandaise pour la re- 
cherche pure (Z. W. O.). 


Manuscrit reçu le 28 septembre 1960. 
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LA FUSION DES CORPUSCULES EN THÉORIE FONCTIONNELLE 


Par Jean-Louis DESTOUCHES, 


Institut Henri-Poincaré, Paris. 


Résumé. — Étude en théorie fonctionnelle relativiste de l’onde barycentrique d’une partie d’un 
système physique et des ondes relatives. Définition d’une partie « fondue » (une partie est fondue 
si les ondes relatives de ses corpuscules sont des constantes), Expression de l’équation d’évolution 
et des équations de condition pour l’onde up d’une partie fondue. Étude de cas particuliers : fusion 
de 2 et de 3 corpuscules, fusion de la totalité du système. Étude de deux fusions successives, équi- 
valence à une seule fusion. Cas particulier de corpuscules de spin 3/2 ou 2 obtenus par fusion 
directe de 3 ou de 4 corpuscules de spin 1/2 ou par fusions successives. Étude de plusieurs fusions 
successives : plusieurs fusions équivalent à une seule fusion. Les phénomènes nucléaires avec 
Nr nv et disparition ou transformation de particules peuvent être décrits par des processus 

e fusion. 


Abstract. — Barycentric wave for a part of a physical system and relative waves in relativistic 
functional theory. ‘‘ Melted ” part of a system : a part is meltedif the relative waves are constant, 
Evolution-equation and condition-equations for the up wave of a melted part P. Particular 
cases : fusion of 2 and 3 particles, fusion of the whole system. Successive fusions, proof of equi- 
valence with only one fusion ; particular cases : particle with spin 3/2 and particle with spin 2 
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obtained by direct fusion of 3 or 4 particles of spin 1/2 or by successive fusions. Several fusions 
are equivalent to only one fusion. Nuclear phenomena with créations and annihilations or 
change of particles can be described by fusion processes. 


1. Introduetion. — M. Louis de Broglie [1] 
en 1934 à introduit sous le nom de fusion un pro- 
cédé formel permettant de former des corpuscules 
de spin arbitraire en mécanique ondulatoire usuelle 
en partant de corpuscules de spin 1/2 ; en parti- 
culier, il obtenait le photon par la fusion de deux 
corpuscules de spin 1/2. Mme F, Aeschlimann [2] 
dans un mémoire récent a montré qu’en théorie 
fonctionnelle, la fusion de deux corpuscules appa- 
raissait comme un processus physique (et non plus 
seulement formel) correspondant à un état de mou- 
vement du système dans lequel il y a évanescence 


. du mouvement autour du barycentre. Nous nous 


proposons ici d'étendre ses résultats au cas d’un 
système de n corpuscules. 


2. Équations fondamentales. — Nous admet- 
trons que les N corpuscules d’un système S sont 
numérotables au point de vue de leur description 
théorique (voir [3]) et que chaque corpuscule C; est 
représenté par une fonction w; ayant des compo- 
santes z;, en nombre fini. La fonction uw, obéit à 
une équation de la forme 


où £; est un opérateur linéaire pour uw; dépendant 
des autres fonctions u, du système et Q; une expres- 
sion analytique non linéaire en w; dépendant aussi 
des autres fonctions w.. L'opérateur £; est de la 
forme 


0 
L; —a 5 + ts; 


& , ue ù 
où $; est un opérateur linéaire ne contenant pas 5 


En outre l’équation (1) est supposée satisfaire au 
principe de relativité restreinte (invariance de 
forme par rapport aux transformations de Lorentz); 
cette condition est en particulier vérifiée si l’opéra- 
teur 6; est l’opérateur hamiltonien d’une équation 
de Dirac : u; et Q; sont alors deux spineurs. 

La fonction uw; ayant des composantes w;, nous 
poserons 

Nate 0e 


et les équations (1) si on explicite leurs composantes 
prennent alors la forme 
« d le À 
ME = (Nja — 65) Ua (2) 
3. Onde barycentrique. — Soit P; une partie 
de "7 corpuscules d’un système S. Nous appel- 
lerons onde barycentrique de la partie P; la fonc- 
tion up, définie par le produit des ondes de chaque 


corpuscule de la partie P%, soit pour une compo- 
sante X : 


UPuose «UP» t) =a ul Uj,aj(P, t). (3) 
Rx 


où P,t désigne un point variable dansl’espace-temps. 

Considérons p parties de S deux à deux dis- 
jointes, P,, ..., P, et soit P la partie réunion 
Hé RER PES 


D 
PEAU) 
k=1 


28 « 


nous appellerons onde barycentrique de la réunion 
des parties P, la fonction uw} définie par 


? 
UPou,. . LiB1- . Br. . An 4 un UPyf:. Br" (4) 


En comparant les formules (3) et (4), on voit que 


D. 
up = II /II Uj\ = IT Uj 
k=1 \56P,, jeP 


et que uw? est ainsi l’onde barycentrique de la 
partie P, donc l’onde barycentrique d’une réunion 
de parties disjointes est identique à l’onde bary- 
centrique de tous les corpuscules appartenant à la 
réunion. Ceci correspond à la propriété classique : 
le barycentre de barycentres partiels est identique 
au barycentre de tous les éléments. 
De (3), en dérivant, nous avons 


dUPpsan.. .@q ÿ ui dUi,a; 
————————— — j,.&. 
dt Py,c,...%q 4 x; dt 


A 


d’où à partir de (2) l’équation de l’onde barycen- 
trique de la partie P; : 


dup 
» Her» « +. g At 
ns RP. . (Nix; — Uj,a; Dj Uj,æ;). 


(5) 


4. Ondes relatives. — Nous appellerons onde rela- 
tive d’un corpuscule C; par rapport au barycentre 
d’une partie P contenant C;, la fonction w,,; définie 
par 

es 
M 


LAS A PORT nt MER ere re (6) 


où mn; est un coefficient positif constant laissé arbi- 
traire pour le moment (sans liaison avec la masse 
au repos de C;) et M la somme des coefficients m; 
des corpuscules de la partie P : 


M =Y My. 
jeP 
De là, on tire 
ut 
M 
Ua = ri. gps. y" (7) 


L’équation (5) peut alors être transformée en 
utilisant la définition (6), c’est-à-dire en remplaçant 
les u;,, par leur expression (7). Posons pour sim- 
plifier l’écriture : 


mMj de) 


Lo R 0j A Up. CY P;.u Pi». 207 (8) 
et (5) prend la forme 
7 Pros» a ne ? 
ME ToT FSr: Ter O3». . sg ieP 
(Na; Er à DR fa erGg Te 84). (9) 
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No1 
En décrivant (6), on a 
7: 
1 : 
dur ja. 0 Lis a: dUj,a; 
dt AN DE 
M; 
Lt Mj re Net DUP a... 
Me Pa. . je dt 


En tenant compte de (2) et (5), puis de (7), on 
obtient l’équation des ondes relatives pour la par- 
tie P; : 


… DUr,j,o.. -&q 


(P) 
LR — © = \N;,4; — 0 12 PNEU PR 
dt ( 1%; p° 7,0 q 9,0% q 
M; 
Tr Ur. fon. 0 > 
M FE CSP 


—1 (P) 
ENS TT Uri. .,o op, y Uri. su (10) 


(L'indice op placé à gauche de X signifie que 
l’opérateur $ contenu dans J opère sur tout ce 
qui est écrit à sa droite, donc ici opère sur 
Uri,oy...œg" ) 

En effectuant le produit des ondes ;,, d’une 
partie P, en tenant compte de (7) et de la défini- 
tion (3), on obtient la relation 


= 


I Creer 


(11) 


5. Partie fondue. — Nous appellerons mouvement 
d’un système S une solution w,, 4», ..., un des 
équations (1) du système S. 

Nous dirons qu’une partie P d’un système S est 
semi-fondue pendant un intervalle de temps (4,4) 
d’un observateur, si le mouvement du système est 
tel qu'introduisant l’onde barycentrique de là par- 
tie P et les ondes relatives pour les corpuscules C; 
de cette partie P, ces ondes relatives u,, se rédui- 
sent à des constantes dans cet intervalle de temps 
(c’est-à-dire que les w,; sont des fonctions cons- 


- tantes du point variable P,t de l’espace-temps dans 


l'intervalle de temps (t, t,) d’un observateur). 


(7,6 PSS EN 


DR Cia. .%n 


La partie P sera dite fondue si les constantes 
Cia. Sont indépendantes des indices ou bien 
nulles, soit 

Cio..m _ Cika.on 


OÙ Xa…2, St un nombre égal à 0 ou à À selon les 
valeurs des indices &,,..æ4,, et il existe au moins 
certaines valeurs x, égales à 1 (composantes 
Ciao, non toutes nulles). 


Selon les interactions, les perturbations subies, 


les conditions initiales fixées pour les équations (1) 
du système, la partie P peut être fondue pendant 
intervalle (4,, {,) et ne plus l’être à des époques 
antérieures ou postérieures à cet intervalle, ou 
bien continuer à rester fondue. Les circonstances 
du mouvement peuvent être extrêmement variées. 


Ale No! 


of DAMES 


Lorsque la partie P est fondue, les équations (9) 
et (10) se simplifient car les wa, étant dans 
ce cas des constantes, elles commutent avec les 
opérateurs. Posons pour simplifier l’écriture : 


NEA ro à (ENS po RÉ. ay) (13) 
où l’indice F placé à gauche de N° signifie que 
dans W on a utilisé la formule (7) et remplacé 
les Ua. Par leurs valeurs constantes. Avec 
cette notation, on obtient à la place de (9) et (10) 
dans le cas de fusion de la partie P : 


dup œ1 (2 
. As. .sXg (P. 
LR TRE = Ua. 0g: Do. eg (14) 
F CP) Mj (P) 
Na; A CRU te 7 Sera = 0 (15) 


Péquation (14) provient de (9) compte tenu de la 
fusion de la partie P, et constitue l’équation d’évo- 
lution de la partie fondue P. Chaque équation (10) 
. pour la partie P se réduit à une équation (15) ; on 
obtient ainsi g équations de condition qui doivent 
être remplies à tout instant de l'intervalle (44, t) 
pendant lequel il y a fusion de la partie P constituée 
par g corpuscules. On remarquera que ces g équa- 
tions sont liées par l'identité : 

ZE Nas — mare) — Sang = 0 (16) 
jeP e 

obtenue en ajoutant membre à membre les équa- 
tions (15) et en tenant compte de la définition (13). 
_ Iln’y en a donc que g — 1 qui sont indépendantes. 
Pour les corpuscules du reste du système, c’est-à- 
dire ceux de la partie S — P, on conserve les équa- 
tions (2) pour chaque corpuscule de cette partie, 
en tenant compte de (7) et (12) pour les fonctions w 
des corpuscules de la partie P qui y figurent. On 
voit que dans une partie fondue on ne peut pas 
distinguer des corpuscules la constituant. La partie 
se comporte comme un tout, en somme comme un 
corpuscule complexe. Pour qu’un système S 
admette une partie P fondue, il est nécessaire que 
l’équation (14) et les équations (15) soient compa- 
tibles (s’il y a g corpuscules dans la partie P). Si 
ces équations ne sont pas compatibles, c’est-à-dire 
si l’on ne peut pas avoir une solution non cons- 
tante Up des équations (14) et (15), alors il 


est impossible que la partie P soit fondue. 


6. Systèmes contenant des parties fondues. — 
Au lieu qu’il n’y ait qu’une partie P fondue, il 
peut y avoir dans le système S des parties dis- 
jointes P;, ..., P, qui soient fondues pendant un 
intervalle (5, t) ; pour chaque partie fondue ?%, 
on écrit une équation (14) et autant d’équa- 
tions (15) que la partie P;4 contient de corpuscules ; 
pour les corpuscules n’appartenant pas à une partie 
fondue, c’est-à-dire les corpuscules de la partie 
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S— UP;, on conserve les équations (2). Si le 
k 


système S contient N corpuscules, au lieu de N 

équations de la forme (2), lorsqu'il y a des parties 

disjointes fondues P,, ..., P; comprenant respec- 
k 


tivement n3, ..., rù corpuscules, les Y n; équa- 
i=1 


tions (2) des corpuscules appartenant à des parties 
fondues sont remplacées par les k équations (14) 
& 


k 
et les 2 ñr équations (15) dont ZX 7"; — k sont 
1 i=1 


indépendantes. Le système S se trouve alors régi 
par N équations indépendantes ; # équations (14), 


2 n; — k équations (15) indépendantes, et 
j=1 


k 
N — Y n; équations (2). 
i=1 


Lorsqu'on arrive aux extrémités de l’intervalle 
(to; t1), ou bien les fusions demeurent, ou bien cer- 
taines fusions cessent d’avoir lieu. On remarquera 
qu’une fusion est un événement objectif (c’est-à- 
dire une propriété intrinséque de la partie P) et 
invariant (c’est-à-dire indépendant du système de 
référence considéré). 


7. Cas particuliers de fusion. — On remarquera 
aussi que les formules (14) et (15) sont valables 
quel que soit le nombre q de corpuscules constituant 
la partie fondue P. En particulier, pour q = 1, 
l’équation (14) se trouve identique à (2) tandis 
que (15) est identique à l'identité (16), compte tenu 
de la définition (13) : on a bien réduction à l’équa- 
tion d’un seul corpuscule. 

Pour g = 2, les deux expressions (15) donnent 
deux expressions opposées en accord avec (16) et 
(13). Les équations (14) et (15) s’explicitent alors en 


nou F F 2 ; 
FREE = U 2)a@( Not Nes — Le En Ke Sue) (142) 


0 — MEN x — x) — mm Nas — am) 


(152) 

Ces équations, à la notation près, sont identiques 
à celles données par M°° F. Aeschlimann [2], elles 
sont identiques aux équations d’un corpuscule de 
spin maximum 1 (spin O0 ou 1). Pour q = 3, les 
trois équations (15) sont liées par l’identité (16), 
deux sont indépendantes et on a comme équations 
de la partie fondue : 


OUT 
in Qt = Wsañr Saby (143) 
Ne. 3e) 7h SU) es 
0 = Na; — RGBY 74, aëv» U=1,2,3). (153) 


Un autre cas particulier est celui où tous les NW 
copuscules du système S sont fondus, dans ce cas, 
P —S, on a une seule équation d'évolution (14) 
pour le système S, et N équations de condition 
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liées par l’identité (16). Si ces équations sont com- 
patibles, on a S 


dUN) C2 œ 

: 41». € NT (NW) 

HSE = UN ),an,. s0n Dose ON (14) 
(N) Mj Q{N) 

0 — N'a; — jo, ON En My Soure.En (15) 


8. Fusions successives. — Considérons un sys- 
tème S admettant des parties fondues deux à deux 
disjointes P,, P,,...,Pz,..., P, composées respec- 
tivement de n,, No, ..., lp, -.., y COrpuscules et 
supposons que ces parties viennent à se fondre en 
une partie P. Chaque partie P, obéit à des équa- 
tions (14) et (15). Si certaines de ces parties P; 
étaient formées d’un seul corpuscule, alors l’équa- 
tion (14) se réduirait à (2), et (15) serait une iden- 
tité (car elle s’identifie alors à la définition (13)) ; 
les équations demeurent valables dans ce cas par- 
ticulier. L’onde barycentrique upo..eg 4€ là par- 
tie P formée par la réunion de parties P, est donnée 
par la formule (4). D’après la propriété de l’onde 
barycertrique énoncée au $ 4, la fonction wa... 
est le produit de toutes les fonctions u,.,, pour tous 


les corpuscules C; appartenant à la réunion des 
parties P,. En dérivant (4), et en tenant compte 
de (14), on obtient pour chaque partie P; 
sut 17) 
Chaque partie P, étant fondue, mais les diverses 
parties P, n’étant pas en général fondues entre 
elles, on doit introduire des ondes relatives pour 
chaque partie fondue P,, soit : 


—My/M 
Url, md Poires D)Oise te sg 

En dérivant par rapport au temps et en omettant 
les indices &,, ..., & qui figurent dans la com- 


posante considérée, on a 


dUr,Py _My,1M dUPy 


E, M; _(MyM)-1 Ddup 
ot L A $ 


QU a rl = 


En tenant compte de (14) pour la partie P, et 
de (17), on obtient : 


dUr,Py ü 
DA 


Mr ar) 
—7 È s( 2}; (18) 


— My,1M ,(P 
up k/ un (sr 


Si les parties P, sont fondues entre -lles pour 
former une partie fondue P, on doit, par définition 
de la fusion, poser que les ondes relatives w;,p, Sont 
des constantes dont les composantes sont indépen- 
dantes des valeurs des indices, mais certaines peu- 
vent être nulles, soit : 


Ur Pps01. 04 Cr ts... age 


Noï 


Dans ces conditions, les équations (17) et (18) 
deviennent : 


rs F Pa) 
ih eo = UP. à 4q > Sa 4 (19) 
F (Pr) Mr À ra{Pr) 
0 — Se han 2 Bo - sg" (20) 


À ces équations, il faut encore ajouter les équa- 
tions de condition provenant de la première fusion : 


GAL KES (Pz) Mi e(Py) 
ES Ce 


(21) 
et y effectuer la seconde fusion. 


9. Comparaison des équations d’évolution. — 
Comparons l’équation d’évolution (19) avec celle 
que donnerait une fusion directe de tous les cor- 
puscules de la partie P : 


dUP,a,..…. N7 


.(P) a 
t Dt UP 109 Sas. C7 (22) 
dt ; (P) LG) Ces. D 
0 = :PNc, —K;.. y M Sa %q (23) 


Selon la définition (8), on a après la première 


fusion et en introduisant l’onde barycentrique up : 


milM} 


—milM 
pi ÿ; UPy,01,. 07 


Pise 2 LE) 


3elPr) —=q 4 


ME MmilM}.My1M & mj/M} 


do. mj/My. My,1M 
D Ur Pioee, UP, ee, 23 Ur, Pros. s%g UP 


d’où 


P — mil M 
Je) + i1My 


(Pr Mr 
JC; 4, +. LPO : 


Ur:P},,œ.. 


a 
done »Xg° 
Si maintenant nous effectuons la seconde fusion, 
nous aurons des fonctions 4, p, qui seront devenues 


des constantes, d’où commutation et, 


13 
ef. Fr p — (24) 


"kg 
D’après la définition (13) pour la partie P, et 
après la seconde fusion, nous aurons 


F e(P#) 
Sr. %g 


SRE PP (A) DE 
EN es) (5 


en désignant par 7% e nombre de corpuscules fon- 
dus en Be Soit n — Z; nr le nombre de corpuscules 
de la partie P, nous avons, d’après (13), en cas de 
fusion directe : 


ÿ 59 F 
Re — > ( PNa; 


P) 
DE 
q jeP Es 2) 


(26) 
nor 


arguments w, pour C4; € P par wy.up" xIM puis, 1 
de la fusion directe, en rendant constant le fac- 
teur wy3. Au contraire FN; S’obtient en rem- 


mais s'obtient en remplaçant dans W,, ; les 


plaçant les arguments uw, pour C4 eP; par u,,z. Hero 


à 5 Tah NS eg 06 


FT ECS 


ee ane de 


Non 


puis lors de la première fusion, en rendant constant 
le facteur &,,4 et °N,,, s'obtient en remplaçant 
dans NV, les arguments wy par up. Ur be puis, 
lors de la seconde fusion, en rendant constant le 
facteur &r,P. Dans ces conditions, P.N;,, s'obtient 


en remplaçant dans NV, les arguments w, par 

M . 
Pme et EN; Sobtient en remplaçant 
mLIM, 


d’abord dans N’;., les arguments u, par Cry UP 
EE, 


tout revient à remplacer les u, par C;,x.up ‘puis, 
vu 

comme up, est remplacé par C;p,.ur! tout 
M M 

revient à remplacer les wy par CCE Epee. 


Mais on a 
mel M MylM\mxk/M 
1e Ur RU. CE Urik-M Pi Up À 
or cette relation subsiste même si les w,4 devien- 
nent des constantes, par suite 


Cr = Cri: CPS 


par suite les arguments substitués aux w, dans N°, 

; F x 
pour obtenir FN, et FN; sont les mêmes, 
d’où la relation 


Ne + NE RS (27) 
En comparant (25) et (26) compte tenu de cette 
relation (27), on voit que l’on a 


P Le P ; 
so ..4q >> PgLEx Xg° (2 &) 


k=1 
Cette relation nous montre que les équations (19) 
et (22) sont identiques, d’où : 


THÉORÈME. — l'équation d'évolution pour une 
partie fondue P — U P, est la même que l’on fasse 
directement la fusion de tous les corpuscules de la 
partie P ou que l’on effectue d’abord la fusion des 
corpuscules des parties P,, P,, ..., P, puis la 
fnstwmdes parties Pi Pa. P,. 


10. Comparaison des équations de condition. — 
Les équations de condition (20) compte tenu de (28) 
s’expriment par 


CE 547 2 XX" (29) 


Les équations de condition (21) pour la pre- 
mière fusion se transforment après la seconde fu- 
sion, d’après (24), en 
xl?) 


(30) 


FF M; F 
OT Ne NT S(Px). 


k 

Après les deux fusions successives, les équations 
de condition sont (29) et (30). Au contraire, dans 
la fusion directe, les équations de condition ont la 
forme (23). 
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Remplaçons dans (30) la quantité FSP» par son 
expression tirée de (29), nous obtenons 
nie LE 


FF 
0 = V';,0; 


M elP) 
TE Sa. (31) 


Œs Xy 
ces équations sont identiques aux équations (23) 
par suite de la relation (27). Donc les équations de 
condition (20) et (30) provenant des fusions succes- 
sives entraînent les équations de condition de la fusion 
directe (23). 

Ajoutons membre à membre les équations (23) 
ou (31) pour 7 e P,, on obtient 


0 — > (Fr, ,@j 
À 


Re (32) 


CENT ET 
or la sommation, d’après (13), n’est autre 
que S®», d’où les équations (20). Remplaçons 
alors dans (31) la quantité S® par son expression 
tirée de (20)/(32), soit : 


= M Fs(Px) 
MX 


s(P) 


nous obtenons les équations (30). Donc les équa- 
tions de condition (23) provenant de la fusion directe 
entraînent les équations (20) et (30) provenant de 
fusions successives. 

Comme ceci a lieu aussi pour l’équation d’évo- 
lution, on arrive à ce théorème : 


THÉORÈME. — Que l’on effectue la fusion des cor- 
puscules de parties disjointes P,, P,, ..., P; d'un 
système S puis la fusion de ces parties entre elles 


? 
(P — U P,) ou que l’on effectue directement la 
k=1 


fusion de tous les corpuscules de la partie P, on est 
conduit dans les deux cas aux mêmes équations 
d'évolution et de condition. 


11. Cas de corpuseules de spins 3/2 et 2. — 
Comme exemple d’application du théorème, pré- 
cédent, considérons une partie P (d’un système S) 
constituée par la fusion de trois corpuscules de 
spin 1/2. Si nous effectuons la fusion directe, nous 
obtenons comme équation de la partie fondue P 
une équation d’évolution (22,), et trois équations 
de condition (23,) dont deux sont indépendantes, 
soit en tout trois équations indépendantes (22,) et 
deux équations (233). 


AVE 6 
iñ Et — — Ug)aby Saôx (223) 
8 mi QG 
0 — Ni — gr — 57 Sabr (233) 
où 
Xy — ;, ta = P, gs — Y, ]=4,2,38 


Si les opérateurs 5; sont ceux de corpuscules de 
spin 1/2 (opérateurs $; issus d'équations de Dirac), 
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la partie P notée (3) se comporte alors comme un 
corpuscule de spin maximum 3/2 (états de spin 3/2 
et 1/2). 

Mais au lieu d’effectuer la fusion directe, on 
pourrait d’abord effectuer la fusion d’une partie P, 
constituée de deux corpuscules de spin 1/2 qui se 
comporte comme un corpuscule de spin maximum 1 
(états de spin 0 et 1), d’où comme équations après 
la première fusion 


ROLE DS ti eue SU 01) 
dt ù dt 
— FN: HSE —— (2 
0 CEE 1 ram JC1,a8 my + M SC (232) 


puis on effectue ensuite la fusion de la partie P, 
notée (2) avec une partie P, notée (1) constituée 
d’un seul corpuscule de spin 1/2 ; d’après le théo- 
rème précédent, on obtient ainsi le même résultat, 
partie P fondue, que dans le cas de fusion directe. 
Dans ce cas, comme équation de la partie P, on a 
l’équation d’évolution (19,) et deux équations de 
condition (20,) et (29;,) identiques l’une à l’autre 
et qui prennent la forme (29,) et deux équations 
de condition (30) identiques l’une à l’autre résul- 
tant de (23,) de la première fusion : 


iñ du 3 :(2 F 
= = Us).agrl" S6Ëy + Sr) (193) 
0 = HART — (my +_M3) Sax (292) 
_ m ? 
0 Pre NT More Tes A es A Poe SR (302) 


My + Mo 


On a donc trois équations indépendantes, une 
d’évolution (19,) et deux de condition (29,) et (30) 
qui sont équivalentes à celles de la fusion directe 
d’après le théorème précédent, ou encore d’après 
les formules (24) et (28) qui se réduisent ici à : 


Fa l2) me Fr sa 


CO) (8) __ Fiç(2) 
HaBr — j,apy ; Sr == "Sc 


Dans (30,) remplaçons S® par 


M + Mi cu) | M3 eu 
DEA CA 
et remplaçons le second terme par son expression 
en 8% tirée de (29,), introduisons S® par (28;), 
d’où (23;) pour 7j = 1. Pour j — 2, on aurait une 
équation: non indépendante de la précédente etissue 
de (30,). Dans (29,) ajoutons et soustrayons 


3 FS remplaçons S2 par sa définition et intro- 
nr p 


duisons S% par (28.), on obtient (23:) en) Por 
D'où les équations (23) à partir de (29,) et (30). 
Inversement ajoutons (23,) membre à membre pour 

— À et 7] = 2, on obtient (29,) compte tenu de la 
définition des 8. Remplaçons S% par son expression 
tirée de cette équation et portons dans (233) pour 
j = 1, on en tire (30). Les deux systèmes d’équa- 
tions de conditions sont bien équivalents. 


N°9 


Donc il n'y a qu'une seule façon d'obtenir un cor- 
puscule fondu de spin maximum 3/2. Si Ci, Ca, C3 
désignent les trois corpuscules constituant la par- 
tie P à fusionner, et si l’indice F marque la fusion, 
d’après le théorème précédent il y a quatre manières 
équivalentes d'obtenir la partie P fondue : 


F(C;, Co; C3) EE F(F(C;, Ca), C3) 


= F(F(C2, Ca), C5) = F(F(C3, Ci), Co) 

De même une partie fondue formée de 4 corpus- 
cules et constituant un corpuscule de spin maxi- 
mum 2 (en particulier un graviton) peut être for- 
mée par fusions successives ou par fusion directe, 
on obtient toujours le même résultat, soit : 


F(C;, Co, Ca, Ci) = F(F(C;, C2), F(C3, Ci)) 
+ E(F(CY, C2, C3), C5) + F(F(F(C:, Co), Ca), Ci) 


ainsi que les cas obtenus par permutation circu- 
laire des indices. Le dernier cas ne résulte pas direc- 
tement du théorème précédent puisqu'il y a trois 
fusions successives, mais d’après ce théorème on a 


F(C1 Cas C3) = F(F(C3, Co); Ca) 


donc on peut effectuer en deux étapes la fusion des 
trois corpuscules, puis encore une fusion 
F(F(C1, Ca, Ca), Ca), donc on obtient encore le même 
résultat après trois fusions. Les quatre cas précé- 
dents correspondent aux quatre décompositions 
de 4 par l’addition : 


=A1+1+1+1=2+92—=38+1—=2+1 +1. 

12. Cas de plusieurs fusions successives. — Sup- 
posons que l’on effectue plusieurs fusions sueces- 
sives de corpuscules et de parties d’un système S. 


Soit F5, une première fusion de corpuscules en des 


parties P,, ..., P,, ..., P, ; soit F, une seconde 
fusion sur les parties Pis. PEL ETES 
parties P°, ..., P;, ..., P, ; d’après le théorème 


précédent appliqué à chaque partie P;, on a, en 
désignant par & la conjonction de fusions succes- 
sives dans l’ordre de gauche à droite, et l'indice é 
marquant que l’on fusionne des corpuscules non 
fondus (ou corpuscules élémentaires) : 


Fay El, pas. (33) 
On a également d’après le même théorème : 
FAN CRRE=EMEES (33;) 
il en résulte que l’on a 
(Fe,1 & Fo) & F3 = Fée (333) 


c’est-à-dire qu'une première fusion F4, suivie d’une 
fusion F, puis d’une fusion F; des parties P1, ..., 
P; PnscPi enr Pr DER identique 
à la fusion directe des corpuscules en les parties 
P:, -.., Pi, ..,, P,. Plus généralément, en omet. 
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tant d'écrire les parenthèses indiquant les fusions 
déjà effectuées on a : 


Fe & Fa & ... & Fr = Fe (33%) 


en effet, supposons que la formule soit valable 
pour # — 1, soit (334_;), d’après le théorème pré- 
cédent, on aura 


Fer & Fr = Fer (335) 


en remplaçant F4x_, par son expression tirée 
de (334), on obtient (33;) ; on a donc 


H (33,) et + (332) & (332) — (332). 


Comme on a établi (33,;) et par suite (332), la 
formule (334) est valable pour #4 — 2 et k — 3. 
En vertu du principe d’induction complète, la for- 
mule est établie pour toute valeur de Z, soit : 


H Vk.(33z). 
d’où ce théorème : 


THÉORÈME. — Que l’on effectue successivement les 
fusions F3 de se en des parties P:, ..., P,, 
puis Fa des parties P,,... P, en parties Pen 
puis F, des parties P;, .. + CRE LEP ON LE etc. 
enfin une fusion F, des parties pers pH 
en des parties P®, , PS le résultat est Le Dee 
que si l’on avait effectué directement la fusion Fix 
des corpuscules en les parties PP, ..., P® (chaque 
partie étant contenue dans une partie PŸT9 et cela 
pour tout j entre O et k — 1 et tout m), les équations 
d'évolution des parties P#, ..., P® et les équations 
de condition sont les mêmes dans les deux cas. 

En résumé, plusieurs fusions successives équi- 
valent à une seule fusion. 

Revenons aux relations (33,) et (33,) ; suppo- 
sons que l’on effectue d’abord la fusion F4, des 


corpuscules en les parties P,, ..., P;, ..., P, puis 
nétusonnt.. desiparties P,,::.,.:P;,....,.2, 
directement en les parties P:, ..., Pi, ..., Py. 
On aura d’après (33,) 
Fé,1 & Fos = Fé,s (3377) 
mais d’après (332) et (33;), on aura 
(Fe, & F2) & F3 = Fes & Fos — Fes (34) 


Fs. est la première fusion à effectuer ; ensuite on 
peut donc d’une manière équivalente, effectuer 
successivement F, et F, ou bien une seule fusion 
F;,,3. On peut alors poser par définition 


F; & F3 —a F2,s 
car F, n’a pas Ge cens si elle n’est pas précédée 
de Fs1et F,, de même F,,, n’a pas de sens si elle 
n’est pas précédée de F1, mais on peut donner 


un sens à F, & F, par la définition précèdente ; avec 
cette notation (34) devient 


(Fé,1 & F2) & F; = Fé,1 & (Fo & F3) (35) 


ce qui établit une relation d’associativité entre plu- 
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sieurs fusions successives. On peut vérifier direc- 
tement cette propriété en écrivant les équations des 
parties fondues. 

Plus généralement, soit une suite Fe, F2, ..., 
Fi, ..., Fy, ..., F, de fusions successives et dési- 
gnons par F;4 la fusion faisant passer directement 
de F;_, à F,. Montrons que l’on a : 


Hé Poe Cible CAP Er OR Tine 
ce Fi& Fat... & Fin & Fix, & Fix & ... & Fine (36) 
D'abord 
LPC CM RCA A 


et les dernières fusions étant les mêmes, il faut 
établir que 


Fe & Fi &... & Fe Pays & Fire (87) 
Tout d’abord, d’après (35), on a 
(Féi—s & Fi) & Fi4s = Fées: & (Fi & Fi+s) (38) 
et on peut définir F;;,, par la parenthèse du 
second membre, soit 
| Fit = Fi & Fit (39) 
Supposons que l’on ait établi 
Fois & Fi & ... & Fins = Fi & Fi 
Fi; = Fiji & F; 
alors on aura d’après (35) 
(Feis & Fi) & F5 = Fei—s & (Fi, & F5) 
et on peut poser + 
Fi, = Fiji & F; 


donc si la formule est valable pour 7 — 1, elle est 
encore valable pour 7, mais la formule est valable 
pour j = & + 1 d’après (38) et (39), donc par récur- 
rence elle sera valable pour J=1 +2, puis pour 
j—=i+3, etc. jusqu'à 7 — k, donc valable 
me k arbitraire (k > i), ce qui établit donc la 
formule (37). : 

On peut remplacer dans (37) Fée; 1 par la 
suite Fes, & Fo .. & F;_, de fusions qui lui 
est équivalente, ce qui établit 


PNG & Eee Ne Gil 


De ES al 


de plus on peut faire suivre la fusion F, ou la 
fusion F;4 (qui a abouti au même résultat) par la 
même suite de fusions Fy41, Frs0, ..., F,, Ce qui à 
partir de l’égalité (40) donne la formule géné- 
rale (36), d’où ce théorème : 


THÉORÈME. — Si on effectue plusieurs fusions 
successives F,, ..., F,,...., F,, on peut remplacer 
des fusions consécutives Fi, Fi,s, ..., FES 
par une seule fusion F;3. 
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13. Conclusion. — Nous arrivons donc à ces 
résultats : 

19 Un système physique S peut avoir des parties 
fondues pendant un intervalle de temps (46, t:). 
Une partie P peut ne pas être fondue avant t, ou 
après Éy. 

20 La fusion est un processus physique réel (et 
non un procédé formel purement mathématique). 

30 Les équations d’une partie fondue sont obte- 
nues déductivement d’une façon parfaitement cohé- 
rente en respectant les exigences physiques. 

40 Si on effectue plusieurs fusions successives, 
on peut remplacer des fusions consécutives par une 
seule fusion qui est équivalente à ces fusions consé- 
cutives. 


No1 


50 Si P est une partie fondue, elle est objecti- 
vement fixée par les corpuscules qui la composent, 
quelle que soit la façon dont on l’a obtenue par des 
fusions successives. 

Les phénomènes nucléaires dans lesquels des par- 
ticules de divers types apparaissent ou dispa- 
raissent, ou sont transformées en d’autres parti- 
cules, peuvent être décrites par de tels processus de 
fusion. Ainsi, en théorie fonctionnelle, ces phéno- 
mènes sont décrits d’une manière différente de celle 
de la théorie quantique des champs usuelle et cette 
nouvelle manière est probablement beaucoup plus 
directement physique. 


Manuscrit reçu le 30 mai 1960. 
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Résumé. — En généralisant les idées développées dans les publications [1], [2], on est arrivé 
ici au formalisme mathématique de la mécanique quantique non relativiste. On a montré que les 
relations d'incertitude, ainsi que les expériences mentales qui sont liées à leur interprétation, 
ont toujours un sens statistique et que l’on ne peut pas en tirer de conséquences sur la conduite 


d’un corpuscuie unique. 


Abstract. — As a summary of the ideas developed in the previous works [1], [2], the author 
has received in this work the mathematical formalism of the prerelativistic quantum mechanics. 
It is shown that the relations of uncertainty as well as the attempts for their interpretation hav 
always a statistical sense and from them one cannot draw a conclusion for the behaviour of a 


separate particle. 


Dans le travail [1] nous avons montré qu’en 
faisant une certaine hypothèse physique sur la 
structure du champ auquel nous avons attribué 
un support matériel (subvac), on est arrivé à une 
image qualitative du mouvement d’un micro- 
corpuscule x (électron ou autre). Dans le travail [2] 
nous avons montré que la description probabiliste 
du mouvement de uw à l’aide d’une fonction 
w(x, y, z,t) > 0 (densité de probabilité) se ramène 
à l’étude d’une fonction (x, y, z, t), pour laquelle 
on à |f[? — w. On a trouvé pour f une équation de 
probabilité [2, (60)], identique à l’équation de 
Schrôdinger. Dans le présent travail nous mon- 
trerons qu’en élargissant dans le même espri‘ 
l'hypothèse physique émise, on peut arriver au 
formalisme mathématique général de la mécanique 
quantique non relativiste. Nous donnerons aussi 
ici un aperçu général des résuitats trouvés en les 
comparant à ceux connus de l'interprétation 
orthodoxe de la mécanique quantique et en étu- 
diant aussi la question de la relation entre l’obser- 
vation et la réalité dans le microcosmos. 


A. Formalisme mathématique général. — 
D'après l'hypothèse émise [1], le champ de poten- 
tiel U(x, y, z), se composant des corpuscules AS, 
peut se décomposer en formations ®,(k4 — 1,2, ....). 
A cause des fluctuations des corpuscules AS et des 
formations ®,, un microcorpuscule p aura un mou- 
vement compliqué non classique. En cas de mou- 
vement lié, le corpuscule x} passe par des états 
stationnaires A3, à chacun desquels correspond 
une densité de probabilité w,(x, y, z) de fonction fx 
correspondante et une valeur moyenne déter- 
minée E4 de l’énergie. Ces valeurs quantifiées de 
l’énergie et les fonctions f, peuvent être considérées 
comme valeurs propres et fonctions propres de 


l’équation de probabilité [2 (57)] (identique à 
l’équation d'amplitude de Schrôdinger). D’après 
la théorie générale des équations différentielles 
elle peut être écrite sous la forme opératoire sui- 
vante : 


Hy — Ey 


où H est l’opérateur linéaire suivant : 


(1) 


ñ? À 
7, + U(x, y; z) 


H— (2) 
et où E est un paramètre. 

Mais le mouvement classique d’un corpuscule x 
possède une suite de caractéristiques mécaniques : 
position r, vitesse v, impulsion p, énergie ciné- 
tique 7, moment cinétique K etc. (En tenant 
compte de Fimportance connue en mécanique 
classique des variables canoniques menant aux 
équations de mouvement de forme symétrique, 
nous nous en servirons ici). En considérant le 
mouvement quantique du corpuseule y lié à l’image 
physique donnée, on doit admettre que les carac- 
téristiques mécaniques de x possèdent aussi des 
variations chaotiques dans le temps, quoiqu’elles 
aient à chaque instant des valeurs bien déter- 
minées. (Ainsi varient par exemple la position r et 
l'énergie Æ de x, où la valeur probable de la posi- 
tion de x et les valeurs fixes moyennes de Æ sonv 
liées à la solution de l’équation (1). Indiquons par B 
une quelconque de ces caractéristiques mécaniques. 
Nous admettrons que la distribution des valeurs 
de B suit aussi des lois statistiques et que dans 
chaque état stationnaire A4 de x elles ont des 
valeurs fixes moyennes, qui possèdent des proba- 
bilités déterminées. Afin d’étudier ces propriétés 
statistiques nous nous occuperons d’abord de la 
question des valeurs moyennes. 
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Comme nous avons vu 
AW = wir, 2, dde dydz =|\Fl" dr, five dei 16) 


représente la probabilité pour que le corpuscule w 
se trouve dans le volume dr, l’intégrale dans la 
condition de normalisation étant prise dans tout 
le domaine de variation des coordonnées x, y, z 
Il s’ensuit, en accord avec la théorie classique des 
probabilités, que la valeur moyenne d’une fonc- 
tion F(x, y, z) des coordonnées de y sera donnée 
par l’expression 


F= frwP ar = [Ft y, 3 


Mais comment sera donnée la valeur moyenne 
d’une fonction G(p) de l’impulsion P(P+, Py, P2), 
ou des coordonnées et de l'impulsion ? Pour 
répondre à cette question, considérons d’abord le 
problème à une dimension suivant OX. 

Dans le cas de mouvement mécanique de y sur 
OX, si Test le temps de parcours d’un intervalle 
Z1, 2 dans un sens déterminé, nous avons la 
«probabilité classique » de présence, donnée par 
[2, (11)], que nous écrirons ici à l’aide de la variable 


Ps; 
CANIN) FRET RD SO CE SEC EA Le 


*{æ, y, 3,1) F(x, y,z,t)dr. 
(4) 


(5) 


SR  Trole 
m 


1 = [7 w(a) de = © — — 
T1 %1 Pax 


On a donc pour les valeurs moyennes classiques 
de deux fonctions données F{(x), G(x, P+) dans 
l'intervalle x,, x, (les intégrales qui suivent sont 
toutes prises dans l'intervalle x,, x) 


mu 


%3 F (x) dx = M ff G{x, Px) At 

G=— : 6 

1P Li Paz Ë 1 L1 Pa (6) 

On aura plus spécialement pour l’impulsion p, et 
pour sa nième puissance 


m 

Pz Dr dæ = 7 (te — %1) (7) 
m _ 

pà = T pi de (7°) 


et de même pour l’énergie cinétique Æ; — p£/2m et 
l’énergie potentielle Æ, — U(x 

+ 1 = U{x) da 
Ee = 3 f pa de, En = % D. (73) 


D’un autre côté nous avons trouvé pour la 
solution f de léquation de la probabilité en cas 


d’approximation classique [2, (37)] 
C' exp (2 : ji TVomtE 0) da) 


W2m(E — U) 
> C’ exp CGJ Pa da) à 
VPe 


(C' est une constante), où, d’après [2, (59)] en 
prenant en (8) le in + pour fixer les idées, 


HER 7e P (2 if Pa Ax — 21 “). (9) 


Îx = 


La condition . normalisation (5) donne 


1 = fret de = cf 


d’où, en comparant avec (5’), on trouve pour la 


constante C” 


CZ CE MIT. (10) 


Cherchons maintenant dans le cas de mouve- 
ment quantique la valeur moyenne d2 l’impulsion, 
en nous servant du principe de correspondance. 
La formule (4) donnant la valeur moyenne de la 
fonction F(x) suggère la voie à suivre pour trouver 
la valeur moyenne de p,. En effet la formule (4) 
appliqué: à l’approximation classique #(9) montre 
immédiatement que si l’on y remplace le facteur F 
par la dérivée par rapport à æ appliquée à 
l’un des deux facteurs W* ou Ÿ, la grandeur p: 
disparaît sous l’intégrale, et la formule trouvée 
devient identique (à une constante multiplicative 
près) à la formule classique des valeurs moyennes 
(7). Alors, en utilisant (9) avec |C’[? tirée de (10), 
calculons l'intégrale suivante 


fr Que 

dx 

1 fa à fe 
exp (-;/ Px de) à exp Gf De da) 

_ . = 

VPx d VPz 
mn dx | 1. dpy ARE Pa + 
VPz 2 Vpà d& hi V px 


(11) 


On voit sur le dernier terme de l’égalité (11) 
que si À — 0, on aura le cas classique et le 
premier membre sous l'intégrale peut être 
négligé devant le second. On voit de plus que la 
formule limite déduite de (11) serait tout à fait 
identique à la formule classique (7) si l’on 
avait remplacé sous l’intégrale en (11) F par 
(A/1) (d/0x) au lieu de 2/04. Par conséquent on 
aura en approximation classique la formule connue 
(7) représentant la valeur moyenne classique pe, si 
l’on prend pour valeur moyenne quantique de 


Non 


Pa la formule suivante, pareille à (4) (les intégrales 
qui suivent sont prises entre les limites — co et 


+ co) 
Pa = [#* P, Ÿ dx (12) 


où, P, est l'opérateur suivant, correspondant à la 
grandeur mécanique p;, 


P, — (ñJi) (D[dx) = — iA(d dx). (13) 


D'un autre côté on vérifie immédiatement que 
si l’on avait interchangé en (12) les places de W* et 
de Ÿ et calculé l’approximation classique de la 
valeur moyenne, on aurait trouvé la même p, (7) à 
condition de remplacer en même temps P, par PX, 
c’est-à-dire on aurait 

a = fr PXP* dx. (12°) 

Les formules (12) et (12’) suggèrent la pensée 
que P, est un opérateur hermitique. 

Il est facile de voir qu’une formule de la forme de 
(12) est valable aussi pour la valeur moyenne de 
p”", c’est-à-dire. 


pa = Jw+paw de, = cn dE. (14) 

En effet en effectuant en approximation clas- 
sique un calcul semblable à (11) on voit qu’en ne 
retenant dans le développement sous l’intégrale 
que les membres d’ordre le plus bas par rapport à % 
(les membres proportionnels à #", puisque À — 0), 
on trouve la valeur classique p°(7;). 

De même on trouve à l’aide du principe de cor- 
respondance qu’une formule pareille à (12) ou à 
(14) est valable pour les valeurs moyennes £,, E, 
(en utilisant (7,), et aussi pour l’énergie totale 


E(x, Pz) = Ee + E,. On aura ainsi pour E 
É = fr E Ÿ dr, 


22 d? 


RSR U(x). (15) 


1 L 
E= Pa + U(c) = 


On vérifie immédiatement la même formule aussi 


pour la valeur G de chaque fonction rationnelle 
entière G(p;) de p, (et de même pour G(x, Pr) en 
remplaçant dans (12) p, par l’opérateur G(P;). 

Les dernières formules peuvent être généralisées 
facilement pour le cas de deux ou de trois dimen- 
sions (ainsi que pour plusieurs corpuscules) en se 
servant de l'indication naturelle que si le nombre 
des dimensions se réduit à un, les formules géné- 
ralisées doivent passer en les formules données 
plus haut. Mais on peut aussi y arriver directement 
en partant du cas à trois dimensions à l’aide du 
principe de correspondance. Considérons par 
exemple le mouvement tridimensionnel du corpus- 
cule x sur une trajectoire 1, décrivant une fois 
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pendant le temps 7 la distance sur 1 entre les 
points Moto, Yw 20), M{(x, y, 2). Prenons par 
exemple la formule de la moyenne classique 
[2, (11°)] pour l’énergie cinétique Æ, 


= a eue M TA 


M j 
E, P di. (16) 


Ji omTU ET Ji, 

D'un autre côté on trouve sans peine dans le 

cas de trois dimensions pour l’approximation de la 

fonction [2 (59)] comme pour [2 (37)] l’expres- 
sion 

C rt [M an 

Pig eyi ts, 1) ARE (2 à [, pa). (16°) 

Alors en utilisant pour Æ une formule sem- 

blable à (15), le calcul à l’aide de (16”) pour # — 0, 

comme dans le cas à une dimension, donnera la 


moyenne classique E, (16), si maintenant on utilise 
les trois opérateurs suivants, généralisations de (13), 


P, — 200 loal, CP 200 0 Pl. 
(17) 


On aura ainsi la formule générale donnant Æ, 


E, he E,Ÿ dr,  E, — P?/2m. (18) 
Dans chaque cas concret semblable au précédent 
on vérifie une formule semblable à (18). 
Admettons maintenant qu’une grandeur méca- 
nique déterminée B est donnée par une fonction 
entière rationnelle B(x, y, Z, Pa, Puy, P:) des &,... 
P», -... D’après l’hypothèse faite au début, la 
grandeur B en chaque état de mouvement prend 
au cours du temps plusieurs valeurs chaotiques 
soumises à une distribution statistique. -Dans 
chaque état stationnaire elles se groupent autour 
des valeurs moyennes déterminées b, qui auront 
des probabilités déterminées. Nous admettrons 
que la valeur moyenne de B dans le cas général 
peut être trouvée de la façon suivante, suggérée 
des exemples concrets considérés et qui, dans le 
cas classique amènera, comme pour la fonction 
G(x, px) en (6), à la valeur moyenne classique LB. 
En remplaçant P>, Py, pe dans B par les opéra- 
teurs (17), on trouve un opérateur linéaire 
B(z, ..., P,, ...). La valeur moyenne B sera 
donnée par l’expression suivante, semblable à (18), 


B = fy* Bo PEUT (9) 


Nous supposons, comme cela a été vérifié dans 
des cas particuliers, que l’opérateur B est aussi 
hermitique. Donc la valeur moyenne B définie 
par (19) sera réelle, comme il est naturel de sy 
attendre. 

Mais comment pourrait-on calculer les valeurs 
moyennes Ÿ, dans un état stationnaire donné 
À;, c’est-à-dire ces valeurs de la grandeur B aux- 
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quelles nous avons admis qu’elles se réalisent en 
effet et qu’elles seront par conséquent trouvées 
dans l’expérience correspondante ? Le cas de 
l’énergie (B — ÆE) qui dans un problème donné a 
des valeurs déterminées représentant des vaieurs 
propres de l’équation différentielle (1) suggère 
l'affirmation suivante : soit B l’opérateur corres- 
pondant à B. Formons l’équation suivante du 


type (1) 
Bq = bq 


où g(x, y, z) sont des fonctions déterminées dans 
un domaine correspondant et b — un paramètre. 
Alors les valeurs propres à, de l’équation (20) sont 
les valeurs que prend la grandeur B dans les états 
stationnaires correspondants (en réalité #4, repré- 
sente la valeur moyenne de B dans l’état station- 
naire A4, les différentes valeurs possibles de B 
dans cet état étant groupées autour de #;). On peut 
appuyer l’affirmation précédente par les raisons 
supplémentaires suivantes : | 

4. L'hypothèse fondamentale dans le présent 
travail d’après laquelle un problème classique se 
transforme en problème quantique possédant des 
états stables et des énergies quantifiées, consistait 
en la décomposition du champ classique continu U 
en formations ®, composées par des corpuscules AS. 
Mais sur la dynamique intime de cette décomposi- 
tion, la dynamique interne des ®;, mêmes, ainsi que 
sur l'interaction entre les corpuscules AS et les ®, 
nous ne savons rien, et à cause de cela l’idée 
émise était formulée de manière très générale. J1 
est donc naturel d’élargir la supposition initiale, à 
savoir, que grâce aux facteurs encore inconnus, il 
existe un type plus compliqué de « quantification 
interne », propres aux formations ®, mêmes, qui 
impose des valeurs moyennes déterminées, c’est-à- 
dire des valeurs quantiques, aussi aux autres gran- 
deurs mécaniques liées au corpuscule. 

2. La solution de l’équation (1) en coordonnées 
convenables (par exemple en coordonnées polaires) 
se décompose dans certains cas en deux ou trois 
équations de une ou de deux variables (séparation 
des variables). Ces équations sont de même type 
que (20) et le paramètre à quantifier, qui y déter- 
mine les valeurs propres, correspond à une cer- 
taine grandeur mécanique appartenant au corpus- 
cule u. Il est naturel de penser qu’une telle équa- 
tion pour la grandeur correspondante peut être 
directement écrite sous la forme (20). 

Les probabilités mêmes avec lesquelles se réa- 
lisent les valeurs b, de B, peuvent ètre trouvées de 
la manière suivante. Admettons pour simplifier 
que l’équation (20) a un spectre discontinu sans 
dégénérescence, b4 et gx étant ses valeurs propres 
et fonctions propres. Les 9, forment un système 
complet, comme cela est connu par la théorie des 
opérateurs linéaires. Alors la fonction Ÿ, solution 
de J’équation de probabilité [2, (60)], peut être 


(20) 
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développée en série suivant les à, sous la forme sui- 
vante 


Firiy ai 2 mt) mr Y; 2). (21) 


En utilisant la condition de normalisation (3) et 


la formule des valeurs moyennes (19) on peut 
écrire 


1= fiyiar= E dite du(e) J qqn d7 = E lan(o)l? (22) 


Passat) dut) [am B qn dr = D bmldm(0)l2. (23) 


mn 


La comparaison de ces deux formules aux 
formules correspondantes de la théorie des proba- 
bilités montre que la constante |[d,(1)|? = |d,(0)/? 
peut être considérée comme la probabilité de 
trouver la grandeur B en un état correspondant à 
l'indice m, de valeur b,. (C’est naturellement, un 
résultat connu de la mécanique quantique, mais 
qui dérive logiquement de nos considérations pré- 
cédentes.) 

Appliquons encore la formule (19) à l’énergie Æ 
correspondant à un état stable fixe 4, du corpus- 
cule u. On aura, en utilisant (2) (B= E,E = H) 


E= [ve vrdr — ff E far = 25 fut de = Fr 


(24, 


donc l’énergie Æ, dans l’état 4, est la valeur 
moyenne de l’énergie totale de x, ce qui est d’ac- 
cord avec notre hypothèse de départ sur le carac- 
tère physique de la notion d’état stationnaire. 

Les résultats exposés jusqu’à maintenant sur le 
mouvement quantique d’un microcorpuscule w 
peuvent être résumés dans les affirmations géné- 
rales suivantes, qui peuvent être acceptées comme 
des positions de départ et sur lesquelles peut être 
fondée la caractéristique quantique totale de la 
conduite du corpuscule p. 

À chaque grandeur mécanique B(x, y, z, 
Pz, Pr, P-), fonction de la position et de l'impulsion 
d’un microcorpuscule 4 dans un champ U(x, y, z) 
on fait correspondre un opérateur B(x, ...,P, ... 
formé par les opérateurs x, ..., P,, ... (17), en le 
remplaçant formellement dans la fonction B. 
B est un opérateur linéaire hermétique. A l’aide 
de B on forme l’équation (20), dont les valeurs 
propres représentent les valeurs quantifiées de B. 
Ces dernières se réalisent avec certaines proba- 
bilités données par les nombres |d,(t)|? dans le 
développement (21). La probabilité de présence 
de y dans le volume dr est donné par |Ÿ{?4r, où la 
fonction Ÿ est une solution de l’équation de pro- 
babilité [2, (60)]. La valeur moyenne de B est 
donnée par l'expression (19). 

Nous avons ainsi montré qu’en partant des 
hypothèses physiques exposées en (1), (2) s’ap- 


? 


pe 


ane, Zap 
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puyant au fond sur l'hypothèse d’un support 
matériel du champ (subvac) possédant une certaine 
structure « atomique », on arrive, pour la descrip- 
tion du comportement d’un microcorpuseule à pos- 
sédantà chaque instant une position et une vitesse 
déterminées, au formalisme mathématique de la 
mécanique quantique (donné ici en représentation 
«æ»). On en déduit, par des moyens mathéma- 
tiques connus, la théorie générale de la repré- 
sentation, mais puisque son développement n’est 
pas lié à des moments physiques nouveaux et 
puisqu'elle est en réalité équivalente à la repré- 
sentation « æ », nous ne nous en occuperons pas ici. 


B. La solution Ÿ de l’équation de probabilité 
et l’onde Ÿ attachée au corpuseule. — Dans les 
considérations précédentes, en cherchant la pro- 
babilité de distribution w et la fonction corres- 
pondante f ou Ÿ caractérisant le mouvement d’un 
(ou de plusieurs) microcorpuscules , nous avons 
abouti à l’équation de probabilité pour f [2, (57)] 


_ ou pour Ÿ [(2, (60)], qui sont identiques à l’équa- 


tion de Schrüdinger correspondante (équation 
d'amplitude ou équation temporelle). Pourtant les 
raisons physiques données ici pour construir: les 
équations de probabilité [2, (57)], [2, (60)] sont 
tout à fait différentes de celles qui sont à la base 
de la mécanique ondulatoire de de Broglie-Schrô- 
dinger. Voilà pourquoi il ne serait pas sans intérêt 
de faire ici une comparaison de ces deux points de 
vue dans quelques cas typiques. Ils auront 
évidemment un traitement mathématique identi- 
que puisqu'ils conduisent aux solutions d’une 
même équation, mais d’un autre côté ils seront 
liés à des interprétations différentes. 

Dans le cas à trois dimensions avec un champ 
U(x, y, z) ne dépendant pas du temps, la solution 
de l’équation de probabilité [2, (60)] a la forme 


L 


La fonction (25) a exactement la forme d’une 
onde qui se propage, d'amplitude g variable avec 
Pendroit et de phase 


(x, y, 2,1) = = Lvx, y, 3) — Et]. (25) 
e Ainsi la fonction Ÿ introduite dans (2) comme un 
moyen mathématique approprié pour la descrip- 
tion statistique de la conduite quantique du cor- 
puseule en mouvement y, est identique à la fonc- 
tion d’onde Ÿ de l’onde de Schrôdinger attachée 
au corpuscule. 

Admettons que le champ est constant et posons 
pour simplifier U — 0. Alors la solution de l’équa- 
tion de probabilité prend la forme (C, constante, 
r=r(x, y; 2)) 

Sy, = CiÉéxp rt |(KT) vi] (26) 
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Donc, en cas de mouvement uniforme de pu la 
fonction Ÿ donnant la description statistique du 
mouvement quantique de x est identique à la 
fonction d’onde de l’onde associée de Broglie de 


longueur À — h/mv. 


Considérons # comme un petit paramètre. Alors, 
comme nous l’avons montré, on trouve l’appro- 
ximation classique de la probabilité w. A l’aide de 
la substitution Ÿ — exp(iS/#) et en développant 
S'en série, S — Y(1%)" S,, on trouve comme pour 


(16) (n = 0,1) 


Pl, y, 2,1) = VE exp! (So —"Et), 
Sol, y, 2) = fPViRE— v) di (27) 


où la fonction S, vérifie l'équation de Jacobi 
1 %, 
(gradS)j+U=E. Me 


L'apparition ici de la fonction S, de Jacobi satis- 
faisant à (27) est une conséquence simple du 
principe de correspondance appliqué à l’équa- 
tion de probabilité [2, (57)]. La même fonction S, 
apparaît comme une phase dans l’onde attachée à u. 
On sait le rôle fondamental qu’à joué la mise en 
correspondance de la phase y de l’onde en cas de 
validité de l’optique géométrique et la fonction S, 
de Jacobi pendant la création de la mécanique 
ondulatoire. 

De même tous les problèmes stationnaires qui 
ramènent, comme nous l’avons montré, la déter- 
mination de la probabilité w aux fonctions propres 
de l’équation de probabilité [2, (57)], peuvent être 
interprétés dans le language de la théorie ondula- 
toire comme un résultat de phénomènes ondula- 
toires, où les états stationnaires 44 d’énergie Æ, ap- 
paraissent comme conséquence mathématique d’in- 
terférence des ondes correspondantes, ne possédant 
pas un sens physique concret (phénomènes de 
résonance). Les cas de passage des corpuscules 
par des orifices (par exemple, expérience de Young, 
réseaux de diffraction, etc.), ainsi que les différents 
cas de propagation avec des obstacles présents 
(diffraction sur un cristal, diffusion sur la surface 
d’un milieu semi-transparent etc.), qui peuvent 
être traités en cherchant la probabilité de distri- 
bution correspondante w — |W{? décrite par l’équa- 
tion de probabilité [2, (60)], trouvent l’interpré- 
tation simple et naturelle d’interférence des ondes, 
comme cela se fait en général en mécanique ondu- 
latoire. 

Considérons le cas d’un mur de potentiel traité 
dans [1] d’une manière qualitative (pour simplifier 
nous considérons le mouvement à une dimension 
suivant OX, le mur étant au point æ — 0). En 
remplaçant la variation brusque de U(x) au point 
x = 0 par une variation rapide mais continue entre 
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les points x, < 0 et x, > 0 et puis en faisant 
Îx — x,] — 0 on trouve (comme il est connu) que 
f(x) et f(x) sont continues pour x — 0. Cela permet 
de trouver la solution de l’équation de probabi- 
Hité [2, (60)] pour chaque x et d’étudier complète- 
ment la distribution de la probabilité de présence 
w(x) de uw (d’une façon analogue au cas d’une 
barrière de potentiel). Le même problème, lié au 
même traitement mathématique de l’équation de 
Schrôdinger reçoit dans le langage de la théorie des 
ondes la représentation élégante bien connue de 
la théorie des barrières de potentiel. 

En donnant en [1] une image qualitative du 
phénomène physique dans le subvac pendant le 
mouvement du corpuscule u, nous avons introduit 
un domaine Ÿ entourant x, où il existe un phéno- 
mène représenté par deux ondes (p,,_), ou bien 
une onde générale +. C’est elle qui permet au cor- 
puscule de connaître quel est le caractère du 
champ aux environs de sa trajectoire et d’influencer 
d’une certaine façon le mouvement de u. De prime 
abord, cette onde a des propriétés et un rôle 
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pareils à ceux de l’onde-pilote de de Broglie. Mais 
ces deux ondes ne sont pas identiques. En effet, on . 
a essayé d'identifier l’onde-pilote avec l'onde 
attachée au corpuscule, ce quiarencontré des diffi- 
cultés insurmontables, tandis que l’onde o(o}, p_) 
correspond à un phénomène physique réel qur est 
un des facteurs responsables de la création de 
l’image statistique du phénomène décrit par la 
même fonction d’onde. 

Ainsi la description mathématique de la -distri- 
bution des probabilités des différentes grandeurs 
liées au mouvement quantique d’un microcor- 
puscule est-elle donnée par des fonctions de même 
nature que celles décrivant en physiqu? classique 
la propagation d’ondes dans un milieu réel. Cette 
analogie merveilleuse évoque la grande unité ir- 
trinsèque des phénomènes appartenant à des do- 
maines différents de la natur2:. Pourtant, comme 
chaque analogie, elle possède aussi un domaine 
de validité dont on devrait indiquer les limites. - 
Nous reviendrons sur cette question à la fin de ce 
travail. 
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Résumé. — Une théorie générale des excitons est formulée du point de vue de la théorie de 
groupes. On considère en particulier le cas de formation optique des excitons dans le CdS. Dans le 
cas le plus simple chacune des bandes de valence donne lieu à une série d’excitons. Si l’on considère 
la structure complexe des bandes au voisinage du centre de la zone de Brillouin, on peut montrer 
que chaque état fondamental est séparé en un groupe de niveaux parmi lesquels le niveau l', et 
les 2 niveaux l'; dans chaque groupe peuvent être créés optiquement en lumière polarisée ayant 
un vecteur électrique parallèle ou perpendiculaire à l’axe cristallin respectivement. 


Abstract. — À general theory of excitons is formulated from point of view of the group theory. 
The particular case of optical exciton formation in CdS is considered. In the simplest case, each 
of the three valence bands gives one exciton series. If the complex band structure near the center 
of the Brillouin zone is considered, one can show that each exciton groundstatesplits in a groupof 
levels, among which the l, level and the two l'; levels can be created opticallyinlight polarized 
parallel or perpendicular to the crystal axis respectively. 


I. Introduction. — Dans les semiconducteurs et 
les cristaux ioniques, on observe souvent à basse 
température des raies ou bandes étroites au voisi- 
nage de la limite d’absorption. Ce phénomène 
général a été rencontré dans un grand nombre de 
cristaux [1]. En particulier, les spectres d’absorp- 
tion pour le CdS, publiés par Gross et Jakobson [2] 
indiquent un nombre considérable de raies ou 


bandes (en tout 17) ; pour celles-ci, il y a dichroïsme! 


comme pour la limite d'absorption 

Ces raies sont souvent attribuées à la formation 
d’excitons : il en est de même parfois des raies de 
fluorescence que l’on observe aussi au voisinage de 
la limite d'absorption. Or, il semble qu’en fait, dans 
bien des cas ces raies doivent être attribuées à des 
impuretés du cristal ainsi que l’ont déjà remarqué 
Hopfield et Thomas [3]. 

En effet, il est très difficile de préparer des cris- 
taux exempts de tout défaut ou impureté. Par 
ailleurs, il existe, quelques critères qui permettent 
de distinguer les raies d’impuretés des raies d’exci- 
tons. 

Dans cet article, nous nous proposons d’étudier 
les caractères généraux des niveaux d’excitons 
dans un cristal et de montrer comment la théorie 
des groupes permet de classer ces niveaux et de 
prévoir la forme des fonctions d’onde associées aux 
excitons. Le cas du sulfure de cadmium sera traité 
à titre d'exemple. Nous en tirerons les règles de 
sélection relatives à la création d’excitons par exci- 


(1) Actuellement School of Science and Engineering, La 
Jolla, Californie. 


tation optique. De plus, nous examinerons com- 
ment la théorie doit être modifiée si l’on veut 
tenir compte du fait que les extrema des bandes de 
valence et de conduction ne se trouvent pas au 
centre de la zone de Brillouin [4]. 


II. Caractères généraux des excitons. — Avant 
d’aborder l’étude des excitons, il est bon de préciser 


.les critères qui peuvent permettre de distinguer les 


raies d’impuretés des raies d’excitons. 

1) Les niveaux d’impuretés sont plus profonds 
que les niveaux d’excitons, les raies d’excitons 
doivent donc apparaître au voisinage immédiat des 
limites d'absorption. En effet plaçons-nous par 
exemple dans l’approximation hydrogénoïde et 
soient m. et m, les masses effectives relatives aux 
bandes de valence et de conduction. La masse 
réduite x associée aux niveaux d’excitons est 
donnée par la formule 1 /u = 1/me + 1/m, et par 
conséquent u < meet u < m,. Les énergies étant 
proportionnelles aux masses, les excitons seront 
moins liés que ne le sont des électrons autour des 
impuretés. 

2) Lorsqu'il y a photoproduction d’excitons les 
raies d'absorption correspondantes doivent être 
très intenses, si la transition est permise ; dans tous 
les cas, elles doivent toujours être plus intenses que 
les raies dues aux impuretés. Le rapport des inten- 
sités varie en gros comme le rapport du nombre 
d’impuretés au nombre total des atomes du cristal. 

3) Les raies d’excitons pourraient se trouver 
élargies lorsqu’elles correspondent à des transitions 
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à partir d’une bande de valence profonde ou jusqu’à 
une bande de conduction élevée. En effet, dans ce 
cas, les niveaux sont métastables. 

4) Enfin parfois on obtient des séries de raies 
d'absorption de caractère plus ou moins hydro- 
génoïde. Si ces raies sont dues à des excitons, les 
intensités varient rapidement avec le nombre quan- 
tique n ; en gros, comme l’a montré Elliott [5], elles 
diminuent comme 1/n°. Au contraire, si les raies 
sont dues à des impuretés, les intensités doivent 
varier lentement, simplement il peut y avoir élar- 
gissement des raies lorsque le nombre quantique n 

croît. 
= Par conséquent, si l’on observe en absorption des 
séries quasi-hydrogénoïdes comprenant un bon 
nombre de termes, on doit les attribuer à la pré- 
sence d’impuretés dans le cristal. 

Ces remarques permettent dans bien des cas de 
séparer les raies d’excitons des raies d’impuretés 
bien qu’il n’existe pas de critère absolu pour les 
distinguer. Les résultats observés dans le Cds, 
s’interprètent en supposant que les raies d’excitons 
sont dues à des transitions permises correspondant 
à des vecteurs d’onde de petit module (centre de la 
zone de Brillouin). 


III. États des exeitons. — Les photons inci- 
dents qui excitent les électrons hors de la bande 
de valence ne communiquent à ceux-ci qu’une 
quantité de mouvement très faible. On peut donc 
supposer que dans la transition les vecteurs d’onde 
sont conservés, le vecteur d’onde global de l’exciton 
formé sera donc nul. 

Naturellement, il n’en est pas de même si dans 
la réaction, il y a émission ou absorption d’un 
phonon. Toutefois, les probabilités de transition 
sont alors faibles et ne peuvent donner lieu à des 
phénomènes observables que lorsque les transitions 
directes sont impossibles. Il en est ainsi si les 
extrema des bandes de valence et de conduction ne 
correspondent pas au même vecteur d’onde de la 
zone de Brillouin. Dans ce cas, il peut y avoir 
encore formation d’excitons, mais ceux-ci ne se 
traduisent pas par des raies mais par des bandes ; 
à chaque type d’exciton correspond alors une limite 
d'absorption bien définie (c’est le cas du germa- 
nium et du silicium). Bornons-nous à l’étude des 
transitions directes sans émission ni absorption de 
phonons. 

Les excitations élémentaires du cristal, produites 
par action de la lumière, peuvent être classées en 
deux groupes : celles qui forment un spectre con- 
tinu et qui correspondent à la création de porteurs 
libres et celles qui forment un spectre discret et 
qui correspondent aux excitons. Les fonctions 
d’onde associées aux niveaux d’excitons de moment 
total nul forment donc des représentations du 
groupe des symétries du cristal. 

Un exciton étant associé à une transition entre 


No 


un maximum d’une bande de valence, et un mini- 
mum d’une bande de conduction, il est très utile 
afin de classer les niveaux d’excitons, d'étudier les 
relations entre les propriétés de symétries atta- 
chées aux fonctions d’onde des excitons et les pro- 
priétés de symétries relatives aux transitions. Exa- 
minons cette question en détails et précisons 
d’abord la forme d’un état d’exciton. 


a) REPRÉSENTATIONS GÉNÉRALES DES ÉTATS. — 
Soit D, l’état de base du cristal. Soit Ck et Chr les 
opérateurs de création et d’annihilation relatifs à la 
bande de valence de type n et soit de même C#z 
et Cr les opérateurs associés à la bande de conduc- 
tion de type m. 

Un état d’exciton pourra être représenté par la 
fonction 


à DE [Zmne Fmnx Che Cnk] ®,;. (1) 


La fonction F,w est une fonction d’onde dans 
l’espace des k. La normalisation de l’état ® 
implique alors la condition 


Ennt| Fmnr|? = 1. (2) 


Remarquons que lorsque nous écrivons la fonc- 
tion ® sous la forme (1) nous ne nous plaçons pas 
dans: l’approximation de Hartree-Fock, valable 
pour un système de particules indépendantes. De 
toutes façons, les excitons étant des excitations 
collectives comme les plasmons, ne peuvent être 
décrits dans cette approximation. Au contraire, 
nous supposons que les états de la bande de valence 
et de conduction sont relatifs à des électrons et des 
trous habillés. à 

Cet habillage des électrons et des trous peut 
théoriquement être réalisé en effectuant une trans- 
formation canonique sur les fonctions d’onde rela- 
tives à l’approximation des particules indépen- 
dant?s. Cette transformation canonique qui tient 
compte des phénomènes de corrélation a pour but 
d'éliminer de l’hamiltonien d'interaction entre élec- 
trons et trous, les éléments de matrice correspon- 
dant aux transitions virtuelles suivantes : création 
ou absorption simultanée de deux paires électron- 
trou, création ou absorption d’une paire électron- 
trou et diffusion simultanée d’un électron ou d’un 
trou. 

L’hamiltonien résultant conserve alors séparé- 
ment le nombre des électrons et le nombre des 
trous (sauf ce qui concerne les transitions réelles) 
c’est ce qui explique que l’on puisse écrire les fonc- 
tions ® sous la forme (1). 

L’interaction entre vraies particules habillées se 
traduit alors par l’existence d’une constante diélec- 
trique. Ce processus est équivalent à la renorma- 
hsation des masses et des charges en électro- 
dynamique quantique. Cet habillage des particules 
ne modifie d’ailleurs pas les propriétés de symétrie 
des fonctions d’onde. Les excitons sont alors des 


à rs pins 
sa a de 
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excitations collectives représentées par les fonc- 
tions ® qui appartiennent au spectre discret. 
Remarquons que les plasmons sont définis de façon 
analogue et peuvent aussi être représentés par des 
fonctions de la forme ®. Excitons et plasmons ne 
se distinguent que par les propriétés suivantes : les 
niveaux d’exeitons sont situés au-dessous de l’en- 
semble continu des excitations individuelles, alors 
que les plasmons sont situés au-dessus [6] (fig. 1). 


Excitons Excitations individuelles 


Plasmon 


Bande interdite 


Examinons les propriétés de symétrie des exci- 
tons, cette étude permettra de préciser la forme 
des fonctions Fr. 


b) TRANSFORMATION DES FONCTIONS Émnk. — 
Les fonctions ® forment des représentations irré- 
ductibles du groupe des symétries du cristal. La 
théorie des groupes va donc nous permettre de 
classer les différents niveaux d’excitons et de 
prévoir la forme des états correspondants. Pour 
plus de clarté, nous résumons dans l’appendice les 
notions et les principaux résultats de théorie des 
groupes, qui sont utilisés dans cet article. 

Soit T l’une de ces opérations de symétrie du 
cristal ; elle transforme les bandes n en n' ; men m' 
et le vecteur k en k’, ®, reste naturellement inva- 
riant puisque l’état de base est supposé non dégé- 
néré. 

TD = Ennk Fnrk Car’ Cr D, (3) 
mais nous pourrons aussi bien écrire : 


Le on * 
TD re ne CmiiCnr De (1) 


Posons alors 
Drux = Cmk Cuk D 0+ (5) 


Nous voyons qu’au lieu de transformer les 
états Dune On peut tout aussi bien transformer les 
fonctions Fm. Autrement dit, il est équivalent 
d'effectuer une opération de symétrie sur le eristal, 
ou d'effectuer la transformation inverse sur les axes 
de coordonnées. Faisons alors tendre Æ vers zéro, 
en ayant soin de garder au vecteur kÆ une direction 
iixe: 

Cette direction sera définie par le vecteur uni- 
taire u (k — uk). 

Les états Dune et Dmnp lorsque k — 0 vont 
tendre vers deux états Dmn(u) et Prrn(u) dégé- 
nérés, dans l’approximation des bandes, puisqu'ils 
se déduisent l’un de l’autre par symétrie. 

De même les fonctions à un électron dmx(r) et 


Unx(r) vont tendre respectivement vers des fonc- 
tions Ym(w)(r) et Yn(w(r) lorsque k — 0. 

Les états mn(u) se transforment comme le pro- 
duit des fonctions à un électron 


Ymteo(r) Ynaw(r). 


En effet les opérations Cmk et Cnx se trans- 


forment comme Ymk(r) et Une(r) respectivement [7]. 
Les représentations associées aux états Dmnn(u) 
sont donc obtenues comme le produit des repré- 


sentations Î', et [, associées respectivement à la 
bande de valence et à la bande de conduction 
pour k = 0. 

Remarquons toutefois que les fonctions Dmn(u) 
dépendent encore de la direction du vecteur k. 
Choisissons alors dans les représentations Te et l, 
(k — 0) une base unique de fonctions ® : soient 
respectivement +, et +, les fonctions de base. Ce 
choix fixe le sens de l’expression ®,,. 

Nous aurons : 


Dhv) = Ce Cv,o Ds. (6) 


Nous définirons alors de nouvelles fonctions ®, x 
en fonction des anciennes Dmnx par les relations : 


Dur = mn umlU) Eva(u) Dnnk- (7) 


Les coefficients a et b sont choisis de telle ma- 
nière que D, tende vers D, lorsque k& — 0. 

Remarquons que les bandes m et n étant définies 
par continuité dans l’espace des k, l’état Dnnx peut 
s’interpréter comme un état obtenu en faisant 
passer un électron de vecteur d’onde k de la 
bande de valence », dans la bande de conduction m. 
L'interprétation de ®,,4 est moins directe, D, 
étant essentiellement fixé par le choix d’ailleurs 


- arbitraire de la base +,, w,. Les indices w et v ne 


correspondent plus alors exactement aux bandes 
d’énergies, mais ceci importe peu. Par ailleurs ®,, 
doit rester normée. 

Ce calcul montre clairement que lorsque plusieurs 
bandes de valence ou de conduction sont dégé- 
nérées au centre de la zone de Brillouin, l’état 
excité doit comporter des contributions de diffé- 
rentes bandes. Dans ce cas, on ne peut plus faire 
correspondre à chaque état d’exciton une transi- 
tion simple entre une seule bande de valence et 
une seule bande de conduction ; l’exciton mélange 
les fonctions associées aux porteurs de même type 
et de masses effectives différentes. 


c) CARACTÈRES DE SYMÉTRIE DES FONCTIONS 
D'ONDE DES EXCITONS. — Nous allons maintenant 
nous attacher à exprimer la fonction d’onde de 
l’exciton formé par le mélange de différentes bandes 
de façon analytique en nous basant sur la théorie 
des groupes. Cette étude nous permettra en parti- 
culier de montrer que les représentations l, des 
différents niveaux d’excitons sont obtenus comme 
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le produit des représentations I4 et l', associées 
aux extrema de la bande de conduction et de la 
bande de valence pour k = 0 et des représen- 
tations l, que l’on peut associer aux « fonctions 
d'onde spatiales ». Physiquement, l’existence de 
ces fonctions d’ondes spatiales est due au fait que 
dans l’exciton, l’électron et le trou forment un 
ensemble lié représenté par une fonction d’onde 
localisée, alors que la polarisation due à cet ensem- 
ble est une onde plane qui se propage librement 
dans le cristal. Dans l’espace réciproque, ces fonc- 
tions d’onde spatiales d}, sont de simples fonctions 


du vecteur d’onde k et ne dépendent pas des 
symétries des bandes. Leurs propriétés de symétrie 
intrinsèques découlent uniquement de leur forme 
analytique en k. 

En résumé, nous allons montrer l’équivalence : 


SEEN se PC STE: (8) 


Utilisons alors la nouvelle représentation dans 
laquelle les états propres des excitons sont expri- 
més comme une combinaison linéaire des états ®,,%. 
Un exciton associé à la représentation [°, est défini 
par la fonction F#%, qui se transforme comme le 
vecteur de base €, appartenant à la représen- 
tation [. Nous allons chercher à l’exprimer sous 
la forme d’une série, chaque terme de la série 
étant le produit d’une fonction des coefficients x 
et v et d’une fonction de vecteur d’onde k. Tout 
d’abord, nous pourrons construire à partir des 
états D, des combinaisons linéaires (a, «) qui 
se transforment comme les variables x, attachées à 
la représentation T avec , € Fe X I. 

Si la représentation l, est contenue n fois dans 


le produit I. X T,, on peut obtenir » groupes de 


combinaisons linéaires indépendantes (a, «),. 


Dz(a,, 4). Dans la suite, pour simplifier nous écri- 
rons simplement ®,(a, «) sans oublier que ce sym- 
bole peut englober plusieurs groupes d’états asso- 
ciés à des représentations isomorphes le. 
DE” = D fév Dyvr. (9) 
Les coefficients /,% sont tout simplement égaux 
aux coefficients de Clebsh-Gordon généralisés. 
fus = Cla,a; eu; 7%, v). (10) 


Les fonctions €#*% sont évidemment normées et 
par conséquent : 


D frs ?= 1. 
LY 


(11) 


Nous introduirons par ailleurs des fonctions ordi- 
naires g(k). Une opération de symétrie en trans- 
formant les vecteurs k, transforme aussi les fonc- 
tions g(k). Désignons par la notation g%(k) les 
fonctions g(k) qui par suite de leur dépendance 
vis-à-vis de X, se transforment comme les va- 
riables x attachées à la représentation 1". Remar- 
quons que l’ensemble des fonctions g'#(k) forme 


No 


un système complet. Comme les fonctions g'°(k) 
sont des fonctions très générales, chacune de ces 
fonctions sera repérée par une indice particulier 
convenable. 

Pour comprendre exactement la signification des 
fonctions g*8(k) il suffit de remplacer le groupe des 
symétries du cristal par le groupe complet des 
rotations. Les fonctions g?(k) deviennent alors des 
fonctions s, p, d, etc. et les indices b et 8 doivent 
être identifiés avec les nombres quantiques habi- 
tuels Z'et m. 

Tout état d’exciton se transformant comme la 
variable x. attachée à la représentation I peut 
s’exprimer en fonction des états ®O$f* définis ci- 
dessus et de certaines fonctions gbP(k). 


Dé* = Lrerexro C(ec; a, & ; bf) gbP(R) DE (12) 
Ce qui revient à poser : Ê 
Eu = Dreraxroe C (ee; a, a; bB) fa dk). (13) 


Les seules quantités non calculables en théorie 
des groupes sont les fonctions g?f(k). Ce sont les 
«fonctions d’onde spatiales » de l’exciton. Pour les 
évaluer il faut résoudre l’équation de Schrodinger 
relative à l’exciton. 

Remarquons que les fonctions d’ondes spa- 
tiales g&Ë(k) doivent obéir à des conditions de nor- 
malisation. Nous pourrons poser : 


Z gb 9 (&)|? = Gap. 


Remarquons aussi que la norme Gw doit être 
indépendante de l’indice $B puisque les fonc- 
tions g?%(k) forment une base pour la représen- 
tation l',. 

Par ailleurs, l’état D devant être normé, nous 
aurons : 


(14) 


D |Féwl* = 1. (15) 
ve 


Cette condition implique la relation suivante : 
rer, xro Gao = 1. 


Les normales G;x mesurent la fraction de l’état 
de l’exciton associée au produit des représen- 
tations ([4 X 1%). La représentation l, est l’une 
des représentations irréductibles du produit PQ x L,3 
elle est liée à la symétrie intrinsèque d’une paire 
électron-trou. La représentation l, est associée aux 
fonctions d’onde spatiale g?8(k). 

Pour un exciton pur, nous aurions d’ailleurs 
sinplement : 


Gab = dasao db,be- (16) 


Dans ce cas, l’exciton ne possède qu’un seul 
groupe de fonctions spatiales g*®(k) et toutes ces 
fonctions se déduisent les unes des autres par 
symétrie. 

À partir des représentations 1, et [}, nous 


No1 


pouvons construire une série d’états d’excitons 
purs de représentation [, avec : 


ET = Tax Dh (17) 


En réalité, les fonctions propres des excitons de 
représentation L, sont formés par des combinai- 
sons linéaires des états d’excitons purs appartenant 
à la même représentation I. Le mélange des 
divers états ne change pas le nombre des niveaux 
et par conséquent nous aurons toujours : 


ZTe=ZE Ta X Lo = Ie X Te X ÊTé (18) 

Il est à remarquer d’ailleurs que si les masses 
effectives des diverses bandes de valence et des 
diverses bandes de conduction dégénérées pour 
k — 0, sont très différentes, les états des excitons 
sont presque purs. Ces remarques permettent de 
classer les niveaux d’excitons. 


IV. Photoproduction des exeitons. —— L’hamil- 
tonien extérieur associé à l’action de la lumière 
incidente sur les électrons, s’écrit : 


He = — (ele) A(E.v) 


& est le vecteur polarisation et À le module du 
potentiel vecteur. L'opérateur v est l’opérateur 
vitesse et peut s’exprimer facilement en deuxième 
quantification. 

L'expression de 7, est absolument générale. En 
particulier elle reste valable si l’on tient compte 
des phénomènes de corrélations et des interactions 
spin-orbite. Dans le cristal, nous pouvons expri- 
mer À, sous une forme un peu différente : 


(19) 


Hy=— ele À À E% pau, (20) 


Les opérateurs v%* sont des combinaisons liné- 
aires des opérateurs vz, Vy, et v,, ils se transforment 
comme les variables x, de la représentation [. 
L'opérateur Æ#, devenant invariant pour toute 
transformation du cristal, il s’ensuit que les coef- 


ficients 8% se transforment comme les variables x, 


de la représentation conjuguée [, (les représen- 
tations [, et I, peuvent être confondues) les opéra- 
teurs v%* peuvent s’exprimer facilement en théorie 
de la deuxième quantification. La probabilité de 
formation d’un exciton représenté par l’état 
s’écrit immédiatement : 


P = r-1< DE|H,|D, >|? (21) 


En exprimant //, de manière explicite, on 


obtient : 
Di 


2 
ñ2 ç2 


A? ADS D°*lve.e|D, = (2. 


(22) 


En effet, les fonctions propres de l’hamiltonien 4 
forment une base orthonormée pour les représen- 
tations irréductibles associées au groupe des symé- 
tries du cristal. Or la fonction v%*|®, > peut être 


EXCITONS DANS LES CRISTAUX DE TYPE WURTZITE 45 


développée par rapport aux fonctions propres qui 
appartiennent à la même représentation [, et qui 
se transforment comme la variable x, liée à cette 
représentation. Par conséquent, elle est orthogonale 
à toutes les fonctions appartenant à d’autres repré- 
sentations [;, ou même correspondant à d’autres 
variables x, de la même représentation 3. | 

L'égalité (4) entraîne la règle de sélection sui- 
vante. Pour qu’un exciton soit créé dans un 
état ®®%* appartenant à la représentation [°,, il faut 
que certaines combinaisons linéaires des compo- 
santes du vecteur & forment une base pour la 
représentation [: 

Nous pouvons appliquer cette règle au cas du 
sulfure de cadmium. Pour ce cristal, au centre 


de la zone de Brillouin, nous avons toujours l; = F',- 
Les combinaisons linéaires des composantes de & 
appartiennent à deux représentations l, et l.. 


Est = E— = (£r — ik) V2 
PAEse = 6e = (+ IVe 
1 EL Es 


Par conséquent, par irradiation optique, nous ne 
pouvons créer, par transition directe, que des états 
d’exciton appartenant aux représentations l', et I". 

En particulier : pour une polarisation perpen- 
diculaire à l’axe du cristal (£, — 0), on obtient des 
états l'; ; pour une polarisation parallèle à l’axe du 
cristal (£3 — Ëy — 0), on obtient des états l.. 

Cette remarque permet de classer les transitions 
permises. Toutefois, nous pouvons être encore plus 
restrictif. Exprimons en effet ®%® par rapport aux 
états D, 


(24) 


e? F2 , F ee . *bÉ 
nee er die Le? Tecra xroC*(es; ax ; bB) ga ‘(K) 


< Duel D, >]? (24) 


Si l’exciton n’est pas localisé, c’est-à-dire si la 
distance moyenne entre l’électron et le trou est 
grande par rapport aux dimensions de la maille 
cristalline, les modules des fonctions g?Ê(k) ne pos- 
sèdent de valeurs importantes que pour des vec- 
teurs k petits par rapport aux dimensions de la 
zone de Brillouin. Pour que cette condition soit 
réalisée, il faut que les masses effectives des bandes 
de valence et de conduction soient faibles et que la 
constante diélectrique soit importante, ce qui 
indique que les excitons sont faiblement liés. Dans 
ce cas, seuls contribuent à la probabilité de tran- 
sition, les éléments de matrice < Of{v®|D, > 
correspondant à des vecteurs k de petit module. 

Nous posons alors : 


< DE pvaeD, > << Dé lves|D, >. 


(25) 


Nous pouvons alors calculer une valeur appro- 
chée de la probabilité de transition; l’approximation 
faite est exactement la même que celle qu’a utilisée 
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Elliott [8] pour le calcul des probabilités de créa- 
tion d’excitons dans le modèle hydrogénoïde. 

Mais l’élément de matrice < ®D9%|v|®, > est 
indépendant de k ; il est donc nul si les représen- 
tations I, et l, sont différentes et plus spécia- 
lement si les indices (a, «) et (e, €) ne coïncident 
pas. 

Nous poserons : 


<< Due DS > = MW. (254) 


La probabilité de transition P s’écrit donc dans 
cette approximation : 
e? 


Le Ale Z rer xryC*lee;e, €; bB) es V8 (k) M. 


(26) 


La représentation [ qui figure dans (25) est 
forcément la représentation unité l, ; à cette repré- 
sentation correspond alors une fonction g{(k) qui 
pour toutes les symétries du cristal se comporte 
comme un scalaire. | 

Finalement, nous arrivons au résultat simple 
suivant : 


e2 
= A2 ce? 


P A2|Ee°|? [ME CILIE (27) 


NE Des $ No1 


Pour que la probabilité de transition ne soit pas 
négligeable, il faut non seulement que la repré- 
sentation JL}, soit compatible avec une transition 
optique directe mais encore que la représentation Le 
soit contenue dans le produit des représentations 
(Te X Lo): | 

Ces conditions suffisent à déterminer les transi- 
tions interdites, si les excitons sont composés. Si 
les excitons sont purs ou presque purs, il faut de 
plus s’assurer que les fonctions g{(k) ne sont pas 
nulles ou négligeables. Dans ce dernier cas, le 
nombre de transitions permises est encore restreint. 

Remarquons enfin que la représentation l, peut 


figurer plusieurs fois dans le produit (Te X nr 


S'il en est ainsi, et si elle y figure 7 fois, nous 


pourrons faire correspondre à ces n représentations 
équivalentes, z éléments de matrices #.. et n fonc- 


tions ge(K) la probabilité de transition s’écrit alors : 


re _ AJ PIE ME eb(R)?. (28) 
j k 


Comme application, nous considérons les niveaux 
d’excitons dans le sulfure de cadmium. L'égalité 
ET = 1, x Fe X ZT permet de les classer. Les 
résultats sont résumés dans le tableau I. (Les ni- 


TABLEAU I 


TRANSITION REPRÉSENTATIONS | REPRÉSENTATIONS DES FONCTIONS D’ONDES SPATIALES g(k) 
DONNANT LIEU ASSOCIÉES T 
A LA FORMATION AUX TRANSITIONS ÿ L 
D’EXCITONS BANDE À BANDE £ 3 à d À mn 
LR Te > T À y È F 
HSE Ge ie DARE ne T: Le EF 
NN Et r'; Ts r': 1 RER r'; FF To Fatbaes 
IDE À Pr EE APR A RROQS PE LES De 2 ES re NS RL LE 
nn. Ts Dm Tr TT) Ti Le Del Te a Le 
r'; r; l'; r; T'; r'; Te 
HE el Es T, Palr, r', r, Fe mi 


veaux d’excitons de la colonne s du tableau I ont 
déjà été donnés par Wheeler [9].) 

Seuls pourront être créés par des transitions 
optiques directes les états l, (pour & 0) et les 
états l'; (pour & +20) 

_Les niveaux soulignés dans le tableau I sont les 
niveaux permis dans l’approximation des excitons 
faiblement liés à condition de supposer que les 
excitons formés sont tous simples. Remarquons 
toutefois que les probabilités de création d’excitons 
doivent être plus fortes pour les niveaux l';(s) et 
l,(s) que pour les niveaux l';(d) et l,(d) et surtout 
les niveaux l';(p) et l'(p). En effet des transitions 
optiques ne sont permises dans ces cas qu’à cause 


de l’influence perturbatrice du champ cristallin et 
du couplage spin-orbite. Par ailleurs, nous voyons 
que l’état s étant le plus symétrique doit être le 
plus lié ; en conséquence, si une transition est per- 
mise au Centre de la zone de Brillouin, on doit 
toujours observer le niveau le plus bas d’exciton. 
C’est donc l’état s qui donnera lieu à l'absorption 
la plus intense et dans bien des cas le premier 
niveau s sera le seul observable par transmission. 
On peut s’en rendre compte en considérant la for- 
mule (27). En effet, si l’on définit la fonction 
enveloppe : 


Fr) D e‘kr g(k) 


N°1 


on voit que la quantité 
F(0) = 2 g4(k) 


k 


qui figure dans l’égalité (27) doit bien être maxi- 
mum pour le premier état ! s pour lequel l'orbite 
est la plus petite. 


V. Influence de la forme des bandes sur la struc- 
ture des niveaux d’excitons. — La classification 
des niveaux d’excitons donnée au paragraphe pré- 
cédent suppose implicitement que les fonctions 
Q%(k) introduites précédemment ne possèdent de 
valeurs appréciables qu’au voisinage de k — 0. En 
fait, il arrive souvent que les extrema des bandes 
de valence et de conduction ne coïncident pas avec 
le centre de la zone de Brillouin. Dans certains cas, 
il est encore possible de créer des excitons par 
transition directe, mais alors la classification des 
niveaux d’excitons doit être modifiée, car 1l faut 
tenir compte des contributions des différentes 
branches des bandes en dehors du centre de la zone 


courbes correspondent à des extrema de la fonc- 
tion e(k) ces extrema sont dans l’espace des k des 
minima ou des cols. Chacune des deux surfaces e(k) 
possède six minima ; les minima de la surface 
inférieure sont les points e, ou les points e;, les 
minima de la surface supérieure sont les points e, 
ou 65. 

En fait dans le sulfure de cadmium l’énergie 
cinétique associée à l’exciton, est égale à l’énergie 
cinétique associée à la transition d’un électron qui 
passe verticalement de la branche 7 de la bande 
de valence à la branche / de la bande de conduction 
Il y a deux configurations possibles pour chaque 
surface l’une correspondant au type I lorsque l’un 
des minima se trouve sur l’axe Ok,, l’autre corres- 
pondant.au type IT lorque l’un des minima se 
trouve sur l’axe Ok, (fig. 3). 

Les représentations irréductibles associées aux 
états excités (e) pour lesquels un électron de vecteur 
d’onde k est passé de la bande de valence à la 
bande de conduction s’obtiennent facilement si l’on 
remarque que lorsqu'on s’éloigne du centre de là 
zone de Brillouin, les représentations [, et 1, se 
décomposent de la même façon. 
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de Brillouin. Cette situation se présente dans le 


sulfure de cadmium. Nous avons montré dans un 


précédent article [4] que les extrema des bandes 
ne coïncident pas avec le centre de la zone (k = 0), 
mais se trouvent au voisinage de celui-ci. Lorsqu'un 
électron de vecteur d’onde & passe de la nappe 7 
de la bande de valence à la nappe ! de la bande de 
conduction, son énergie cinétique (énergie dans 
l’approximation des bandes) varie de la quantité 
e(k). 
e(k) = Ecu(k) — Ev,5(R). 


Référons-nous à la structure de bande donnée 
dans un précédent article [4] (fig. 4) et considérons 
les transitions qui ne font intervenir qu’une bande 
de valence et une bande de conduction, chacune 
des bandes étant formée de deux nappes qui se 
coupent au centre de la zone de Brillouin. Les 
variations des énergies e(k) suivant les directions 
des trois axes Ok,, Ok, et Ok, sont représentées 
sur la figure 2. : 

Les points e, e, e et e qui figurent sur les 


e(k) 


Suivant Ok, 


= 
7 


L x 
ns >, + > 


Suivant Ok, 


- 


7 
T7 
l'y 


7.) + TK). 


Les résultats seront identiques quelle que soit la 
bande de valence considérée, les représentations 
associées aux états (e) sont les suivantes [10] : 
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Dans la direction Ok, 
D: X DS 
D X 2 


DAT O2r 
2 
2 C2 He 


Dans la direction Ok, 
Talk) X Talk 
Te A en 
Talk) X Taf— k) 
Talk) X Ts(— À) 
Si l'énergie cinétique moyenne des excitons est 
petite par rapport aux séparations entre les di- 
verses courbes, il peut être intéressant de considérer 
la fonction d’onde globale ® d’un exciton comme 
la somme des 6 contributions associéés à l’un des 
groupes de six minima dégénérés. Nous suppose- 
rons que les fonctions d’onde spatiales associées à 
ces excitons sont des scalaires (représentations >; 
ou T,), c’est-à-dire localement des fonctions {s. 
Les fonctions d’excitons @ associées à chaque 
minima ont donc même représentation que les 
états excités (e) d’énergie ex, e2, &, e, qui leur 
correspondent. 


€1 


C2 


0 
1 


2 


No1 


Pour le groupe I, la fonction @ associée au site 7 
de représentation X, sera désignée par la nota- 
tion 4%, Pour le groupe IT, la fonction @ associée 
au site j de représentation 7, sera désignée par la 
notation ®;r,- 

On peut toujours choisir les coefficients de nor- 
malisation de ces fonctions de façon sûre, puisque C4 
étant l’opérateur de rotation autour de l’axe hexa- 
gonal, on aura toujours : 


Ce Pin — Pituzn tt C6 Piln = Pi+iTy- 


Les fonctions globales des excitons ® sont des 
combinaisons linéaires des fonctions d’excitons 
et forment des représentations associées au centre l 
de la zone de Brillouin (groupe de l’hexagone plan), 
nous poserons : 


DZ, = D dj Pj,Zn I 
dl 
ou 
DTn = Day Py.rr IT 
3 


Les coefficients « sont choisis de façon que les 
fonctions ® forment des représentations irréduc- 
tibles. 


TABLEAU II 


CC A CC) 
(1, 1, 1, 1,1, 1)/V6 > 
CRETE RE A Vr 
(L. j, k, l, j, k)/V6 

(LEE Le 


(2, Re ja it k, NET es 
EAN 
(l, — j, k, —1, j, — k)/V6 


j et k sont les racines cubiques de l’unité 
j = —1/2 + iV3/2 
k = —1/2— iV3/2. 


Des excitons sont créés : en lumière // s’ils ont 
une représentation l°, ; en lumière j s'ils ont une 
représentation [;. 

Le tableau IT indique que si les vallées de la 
fonction e(k) sont assez profondes, l’état fonda- 
mental des excitons associés aux transitions de 
l’une des trois bandes de valence à la bande de 
conduction est décomposé en deux groupes de 
& niveaux. 

En lumière polarisée parallèlement à l’axe du 


REPRÉSENTATIONS ASSOCIÉES 
AUX FONCTIONS 


€ OÙ € € OÙ € 
1(5) où HUT;) |I(3,) ou INT.) 
AU 2 | le 
D A 
PRE PET 
r'; Tr; Tr; r'; 


cristal, un seul niveau l°, est observable par absorp- 
tion, mais en lumière polarisée perpendiculaire- 
ment à cet axe deux niveaux ['; doivent être. 
observables. 

Les niveaux excités peuvent être traités de ma- 
nière analogue. Comparons maintenant les ta- 
bleaux Let IT, les représentations indiquées dans le 
tableau I sont obtenues en supposant implici- 
tement que les extrema des bandes de valence et 
de conduction sont extrêmement voisins du centre 
de la zone de Brillouin. Les résultats du tableau II 
tinnent compte du fait que les extrema peuvent 
être assez éloignés du centre de la zone de Brillouin. 
Dans ce dernier cas, la structure des états est plus 


_ Ne1 


complexe. Les états de ces excitons peuvent d’ail- 
leurs être décomposés en une somme d'états 
« d’excitons purs » En fait, on pourrait établir 
une correspondance entre les deux tableaux. 

Les classifications des excitons que nous avons 
données pourraient permettre de calculer les posi- 
tions et les intensités des pics d'absorption obser- 
vables. : 

Dans le Cds, il existe plusieurs séries d’excitons 
provenant de diverses bandes de valence, mais 
celles-ci peuvent empiéter les unes sur les autres. 
En effet, dans le cadre de la théorie simplifiée, on 
peut considérer trois séries d’excitons purs ou 
presque purs provenant de trois bandes de valence. 
Si, par contre, on tient compte de la complexité des 
bandes, qui se traduit par l’existence d’un ensemble 
de transitions privilégiées au voisinage du centre 
de la zone de Brillouin, on doit envisager la possi- 
bilité de création par absorption optique de trois 
séries d’excitons [', et de 6 séries d’excitons F;, 
observables respectivement en lumière polarisée 
parallèlement ou perpendiculairement à l’axe du 
cristal. 

Devant la complexité et le nombre considérable 
d'états possibles, lexamen des résultats expéri- 
mentaux doit être très critique. On peut supposer 
en premiere approximation que seulsles états de base 
de l’exciton sont créés par irradiation lumineuse et 


_observable en absorption puisque ce sont les tran- 


sitions les plus intenses. Dans ce cas, les bandes 
d'absorption observées peuvent avoir des origines 
différentes. Il est plausible d’attribuer des séries 
hydrogénoïdes de raies aux centres d’impuretés peu 
profonds et ne considérer comme dues aux excitons, 
que les transitions très intenses, au voisinage de la 
limite d'absorption. 


APPENDICE 


1) Toutes les substitutions linéaires associées 
aux opérations d’un groupe G peuvent être trans- 
formées de telle façon que la matrice associée à la 
substitution soit unitaire [11]. Les substitutions 
sont alors dites de forme standard. 
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2) Soit Gun groupe irréductible de substitutions 
linéaires de forme standard agissant sur les va- 
riables æ, ... «, et soit À un groupe isomorphe à G 
agissant sur les variables y,... yn ; le critère pour 
qu'il existe une forme bilinéaire en x et y, invariante 
lorsque les x et les y subissent des transformations 
correspondantes est le suivant : il faut et il suffit 
que les variables y obéissent aux transformations 


du groupe G conjugué de G. La forme bilinéaire 


invariante est alors : 
j=m 


2 y. 
j=1 
3) Soit G un groupe irréductible de substitutions 
linéaires agissant sur les variables x, ... 2, La 
norme ÏV du vecteur de composantes x; peut se défi- 
nir ainsi : | 


j=n 
N = X xŸ x;elle est invariante pour toutes les 


=1 
transformations du groupe. 

Un vecteur est dit normé si sa norme est égale à 
l’unité. 

4) Soit un vecteur Z de composante Z, appar- 
tenant à la représentation l, et obtenu en compo- 
sant les vecteurs x et y de composantes æ et y, 
associées aux représentations l, et Te nous le repré- 
senterons ainsi : 


Z = Y caal bflcy) Ze Yy 
BY 
El xl. 


Les coefficients c (ax|bB|cy) seront choisis de telle 
façon que si les vecteurs + et y sont normés, le 
vecteur Z soit lui aussi normé. 
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EXPOSÉ ET MISE AU POINT BIBLIOGRAPHIQUE 


DÉTECTION ET DOSAGE DE LA RADIOACTIVITÉ NATURELLE DE L’ATMOSPHÈRE 


Par DaniEL BLANC, Jacques FONTAN, RENÉ SOULHIER et Gizgerr VEDRENNE, 


Centre de Physique Nucléaire, Faculté des Sciences de Toulouse. 


Résumé. — Après avoir rappelé quelle est la nature de la radioactivité naturelle de l'atmosphère, 
on passe en revue les divers procédés de détection et de dosage de cette radioactivité : les méthodes 
directes (chambre d’ionisation, compteur proportionnel, scintillateur) et les méthodes indirectes 
(précipitation thermique, précipitation par centrifugation, précipitation électrostatique, filtration). 

Compte tenu de la diffusion de l’émanation et des conditions météorologiques qui modifient la 
radioactivité naturelle de l’air, on donne enfin l’ordre de grandeur de cette radioactivité. 


Abstract. — After recalling the nature of natural atmospheric radioactivity, we survey the 
various methods of detection and dosimetry of this radioactivity : direct methods (ionization chamber 
proportional counter, scintillator) and indirect methods (thermal, centrifugal and electrostatic 


precipitations, filtration). 


Finally, taking into account the diffusion of the emanation and the meteorological conditions 
that modify the natural radioactivity of air, we give the order of magnitude of this radioactivity. 


Origine et nature de la radioactivité naturelle de 
l’atmosphère. — En 1901, Elster et Geitel [El-1] ont 
mis en évidence l’existence d’une radioactivité de l’air 
atmosphérique. On reconnut assez rapidement qu’elle 
a pour origine deux gaz, le radon (?#Em) et le tho- 
ron (2?Em), qui sont deux isotopes de l’élément 86, 
appelé « Emanation ». Le premier appartient à la 
famille radioactive naturelle (uranium-radium) : il est 
produit par la désintégration alpha de *#$Ra. Le second 
appartient à la famille du thorium : il provient de la 
désintégration du thorium X (2%4Ra). 

Il existe un troisième isotope gazeux de l’élément 86, 
l’actinon : il appartient à la famille de l’actinium et 
résulte de la désintégration de l’actinium X (Ra); sa 
période est très brève (3,92 secondes)et l'atmosphère ne 
contient que de très faibles quantités de ce gaz ; la 
radioactivité correspondante est négligeable. 

Le radon et le thoron se forment à l’intérieur du sol, 
à partir de l’uranium et du thorium qui s’y trouvent en 
faibles concentrations ; ils diffusent ensuite vers la 
surface et se dégagent dans l’atmosphère. Les maté- 
riaux de construction habituels (briques, pierres, 
béton, etc.) contiennent des traces d'uranium ét de 
thorium : du radon et du thoron se dégagent à linté- 
rieur des bâtiments. Bien plus, selon Anderson et ses 
collaborateurs [An-1], [An-2], [An-3], une partie du 
radon présent dans l’atmosphère des régions indus- 
trielles proviendrait de la combustion du charbon. 

Dans l’atmosphère, le radon et le thoron se désin- 
tègrent à leur tour pour donner deux chaînes de des- 
cendants radioactifs (fig. 1 et 2). Tous ces descendants 
sont solides ; lors de leur formation, ils possèdent une 
grande mobilité et se fixent rapidement sur les aérosols 


et les poussières en suspension dans l’air, comme l'ont 
montré Labeyrie [La-2], ainsi que Carmichael et Tunni- 
cliffe [Ca-2]. Voir aussi : [Bu-1], [Po-1]. 

Ces aérosols, de très faibles dimensions, peuvent 
rester longtemps présents dans l’atmosphère, en l’ab- 
sence d’une cause extérieure entraînant leur préci- 
pitation au sol, et la radioactivité des descendants 
s’ajoute à la radioactivité du radon et du thoron ; le 
calcul de l’ensemble est assez complexe [Co-4], [Ha-1], 
[He-1}, [So-1]. 

On a mesuré l’activité de tous les descendants du 
thoron. Pour le radon, celles de RaA, RaB, RaC et 
RaC” ont été déterminées (fig. 1). La période du des- 
cendant suivant, le RaD, est beaucoup plus longue 
(22 ans): les aérosols sur lesquels il se fixe sont 
entraînés vers le sol, par leur poids ou par les préci- 
pitations, avant qu’un équilibre ait pu être atteint ; 
d’après Blifford et ses collaborateurs [B1-3], la période 
de séjour des aérosols dans l’atmosphère serait en effet 
voisine de 10 jours. Il en résulte que la radioactivité 
due aux descendants à vie longue (Ra), RaE, RaF) 
est très faible [Ki-1]. 


… Procédés de détection. — CHOIX DE LA MÉTHODE. — 
Selon la valeur de la concentration des aérosols ou des 


gaz radioactifs dans l’air, les méthodes utilisées sont 


très différentes. 

19 Si l’activité est supérieure à 10—? curie par litre 
(plus de 2 220 désintégrations par litre et par minute) 
on peut plonger directement dans l’air un détecteur, 
compteur de Geiger-Müller ou scintillateur, par 
exemple [Ak-1]. De tels niveaux d’activité sont supé- 
rieurs à 10 fois l’activité tolérable sans danger, d’après 


mi druide ht 


NOEL 


les normes de l'EURATOM. Ils ne peuvent être 
atteints qu’à l’intérieur des mines d'uranium. 

20 Si l’activité est de l’ordre de 10—1? curie par litre 
(environ deux désintégrations par minute et par litre), 


226. 


Ra 
1590 ans 


218 
RaA ,,Po 
8,05 min 


26,8 min 
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on peut enfermer une certaine quantité du gaz étudié 


dans un détecteur convenablement choisi. Il s’agit 
d’une méthode directe de mesure, analogue à celle 
employée dans létude des radioéléments gazeux 


82 
Stable 


F1G. 1. — Chaîne de désintégration des descendants du radium. 


ThX “Ra 
| 3,64 jours 


3,1 min 


F1. 2. — Chaîne de désintégration 
des descendants du thorium %, 


comme le tritium (%H){[B1-1], [Sh-1] ou #C sous la 
forme d’anhydride carbonique [St-3]. 

30 Pour des concentrations nettement plus faibles 
(40-25 eurie par litre, ou moins) la méthode précédente 
ne peut plus s'appliquer. On peut concentrer le gaz 
radioactif avant de l'envoyer dans le détecteur où opé- 
rer de façon indirecte en rassemblant sur un collecteur 
- les aérosols présents dans l’atmosphère ; l’activité 
recueillie est ensuite analysée. 


MÉTHODES DIRECTES. — l’air est envoyé à l’intérieur 
d’un détecteur (chambre d’ionisation, compteur pro- 
5 . . RS ? 
portionnel, scintillateur spécial). 


19 Chambres d'ionisation. — Les électroscopes uti- 


lisés par les premiers expérimentateurs tels que 


Ashman [As-1]et Eve [Ev-3] rentrent dans cette caté- 
gorie. Après concentration sur du charbon actif, très 
poreux, ou condensation à la température de l'air 
liquide, l’émanation était envoyée à l’intérieur de 
l’électroscope ; l’ionisation produite par les particules 
alpha avait pour effet de décharger l’élèctrode cen- 
trale. De la vitesse de ce phénomène, on déduisait la 
quantité d’émanation introduite ; l'appareil était préa- 
läblement étalonné à l’aide d’une source de radium, 
émettant par unité de temps une quantité connue de 
radon. 

La méthode fut progressivement perfectionnée : le 
gaz à analyser était envoyé dans une chambre d’ioni- 
sation dont l’électrode centrale était reliée à un élec- 
troscope ou à un électromètre [La-4], [Sa-1]. De nom- 
breux montages utilisant des chambres d’ionisation 
sont encore employés : ils diffèrent en général par la 
façon dont on mesure le courant d’ionisation : on utilise 
des tubes électromètres [Fi-1], des amplificateurs de 
courant continu | De-3], à condensateur vibrant [Ba-2], 
[Hu-3], ou même des électromètres classiques [An-2], 
[Ja-2],[He-51]. 

On peut aussi faire fonctionner la chambre en impul- 
sions et détecter le passage de chaque particule [Ca-21, 
[Cu-3] ; la hauteur d’impulsion étant proportionnelle 
à l’énergie perdue par la particule dans le gaz du 
compteur, on peut séparer les différents émetteurs 
alpha en associant la chambre d’ionisation à un sélec- 
teur d’amplitudes [Fa-2], [Fa-3], [Hu-1]. 
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Dans le cas d’activités très faibles, on peut utiliser un 
montage différentiel à deux chambres d’ionisation, dont 
la sensibilité est plus grande : l’une d’entre elles, 
fermée, enregistre l’ionisation due aux rayonnements 
parasites (rayons cosmiques, ete.) ; le courant qu’elle 
délivre est automatiquement soustrait de celui de 
l’autre chambre, dans laquelle on envoie le gaz à 
étudier [Ev-1], [Ja-1], [Wi-2]. Un tel montage peut 
fonctionner de façon continue. On trouvera sur la 
figure 3 l’ensemble réalisé par Labeyrie et Weill [La-3]. 


FiG. 3. — Chambre d’ionisation différentielle 
à fonctionnement continu. D’après [La-3]. 


Les! chambres d’ionisation et les électromètres con- 
timuent à être employés pour la mesure du radon et du 
thoron qui se dégagent des eaux thermales [Du-1], 
[Ku-2], [Tr-1]. 


20 Compteurs proportionnels. — Ils sont souvent uti- 
lisés pour mesurer la quantité de radon et de thoron 
qui se dégage des roches, des liquides, etc. [Be-21, 
[Br-2], [Si-2]. Le gaz est envoyé à l’intérieur d’un 
compteur proportionnel fonctionnant de manière clas- 
sique [BI-1]. Malgré sa grande sensibilité, la méthode 
ne s’applique pas à la radioactivité de Pair lui-même : 
l'oxygène et l’anhydride carbonique qu’il contient sont 
les sources d’un grand nombre d’ions négatifs qui pro- 
duisent des impulsions parasites. Nous citons cette 
méthode pour mémoire. 


30 Scintillateurs. — Cette méthode, plus récente, a 
trouvé de nombreuses applications. Le gaz est envoyé 
dans une chambre dont les paroïs internes sont recou- 
vertes d’un scintillateur pour rayons alpha, en général 
une poudre polycristalline de sulfure de zinc activé à 
l’argent [An-4]. Une fenêtre, ménagée dans la paroi, 
est occupée par la photocathode d’un photomulti- 
plicateur qui permet le comptage des scintillations. 
Ces scintillations sont produites par les particules alpha 
du radon, lorsqu'elles heurtent la paroi, et aux des- 
cendants du radon émetteurs alpha (RaA et RaC”), qui 
se fixent très rapidement sur les parois de la chambre, 


N°1 


comme l’a montré Labeyrie [La-2]. La détection du 
thoron, moins abondant, est plus délicate. 

Les chambres sont souvent des ballons de verre ou 
de pyrex. Il semble toutefois que le bruit de fond du 
détecteur soit très diminué lorsque l’on utilise des 
parois de Kovar ou d’acier inoxydable [Lu-1]. Le 
volume de l’appareil est en général de quelques cen- 
taines de em, [Ku-3], [Va-1]:0n ne peut alors mesurer, 
avec une bonne précision, que des concentrations de 
radon supérieures à 10—1% curie par litre, telles que 
celles employées dans certaines études de biologie 
[Da-21, [Ma-4], [Ma-5], [Va-1]. 

Pour mesurer la concentration du radon dans l’air, 
le service d’Électronique Physique du Commissariat à 
l'Énergie Atomique a mis au point un ballon de verre à 
fond plat, d’un volume de 1 litre, dont la paroi interne 
est recouverte de sulfure de zinc activé à l’argent : un 
tel volume permet de mesurer les concentrations de 
radon habituellement rencontrées dans l’atmosphère. 
Un ballon scintillateur analogue, d’un volume de 
500 cm3, a été utilisé par Pradel et Billard [Pr-1] pour 
le dosage du thoron au voisinage des mines de thorium ; 
on mesure ainsi les concentrations en thoron supé- 
rieures à .10—1% curie par litre. 


MÉTHODES INDIRECTES. — Le principe reste toujours 
le même : on rassemble sur un collecteur convenable les 
produits de désintégration du radon et du thoron en 
suspension dans l’air. Ces produits sont, le plus souvent, 
portés par des aérosols et les procédés industriels de 
dépoussiérage s’appliquent : précipitation thermique, 
précipitation par centrifugation, précipitation électro- 
statique, filtration. Nous examinerons successivement 
ces divers procédés. | 


19 Précipitation thermique. — Le principe de cette 
méthode est voisin de celui de la séparation isoto- 
pique : l’air chargé de poussières passe lentement entre 
deux parois, entre lesquelles existe un certain gradient 
de température. Les molécules situées au voisinage de 
la paroi chaude ont une énergie cinétique moyenne 
plus grande que l’énergie cinétique de celles situées au 


, 


Paroi chaude LE 


Fi1G. 4. — Principe de la précipitation thermique. La par- 
ticule P est déviée vers la paroi froide. D’après [De-2]. 


voisinage de la paroi froide. Finalement, les poussières 
et aérosols P transportés par l’air recevront davantage 
de chocs du côté chaud que du côté froid : ils auront 
tendance à se coller sur la paroi froide (fig. 4), [Br-1], 
[Wa-1], [Wi-6]. Ce procédé présente l’avantage de 
collecter les plus petites particules : un tel appareil 
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aurait une efficacité très voisine de 100 % pour des 
perticules de diamètres inférieurs à un micron [Pa-1]. 
Malheureusement, le débit ne peut guère dépasser 
10 cm° par minute, ce qui limite considérablement son 
emploi. 


20 Précipitation par centrifugation. — L'air est aspiré 
à travers une série d'ouvertures en face desquelles sont 
placées des plaques qui dévient brutalement les filets 
d'air. Cette déviation projette les poussières sur la 
plaque où elles sont recueillies par un enduit adhésif 
(colle, graisse de silicones, ete.) [Gr-1], [Ma-7], ou par 
un liquide mouillant [Da-3], [Fr-1]. A l'inverse de la 
précédente, cette méthode possède un rendement beau- 
coup plus grand pour les grosses particules que pour 
les petites. On peut l’associer à un autre procédé (fil- 
tration, précipitation électrostatique) pour pallier cet 
inconvénient. En faisant suivre d’une rondelle de 
papier filtre quatre jets successifs, Laskin [La-5], [St-2] 
est parvenu à obtenir un rendement variable mais 
parfois voisin de 100 %, pour des aérosols de 0,2 micron 
à 50 microns de diamètre. 


39 Précipitation électrostatique. — Nous avons groupé 
dans ce paragraphe les appareils qui utilisent la charge 
des particules pour les envoyer sur un collecteur grâce 
à un champ électrique. Selon l’origine de cette charge, 
il faut cependant distinguer deux catégories : à la pre- 
mière appartiennent les appareils dans lesquels les par- 
ticules se chargent en passant dans un champ élec- 
trique intense ; dans la seconde sont rangés les .appa- 
reils qui utilisent la propriété d’ionisation des atomes 
venant de subir une désintégration alpha. 


Première catégorie. — L'industrie emploie un certain 
nombre de procédés de dépoussiérage électrique dont 
le principe est le suivant : l’air à dépoussiérer est 
envoyé dans une région où règne un champ électrique 
intense, créé par une forte différence de potentiel 
(10 000 à 100 000 volts) établie entre un fil [Ha-5] ou 
une série de fils et un système d’électrodes. Au voisi- 
nage du fil, air est ionisé par effet couronne : les ions 
ainsi créés ionisent à leur tour les particules en sus- 
pension dans l’air ; ces particules sont entraînées sur 
l’électrode collectrice. Dans tous ces appareils, on place 
l’électrode collectrice à la haute tension positive : ainsi, 
le courant de décharge est légèrement plus élevé et la 
tension d’amorçage entre les électrodes sensiblement 
augmentée [De-11. 

Le mécanisme de l’ionisation des particules a fait 
l’objet de plusieurs théories. On pourra se reporter aux 
travaux de Pauthenier et de ses collaborateurs [Ch-3], 
[Pa-2], [Pa-3], de Loeb [Lo-2], de Lowe et Lucas 
[Lo-3], ete. 

L'application de la précipitation électrostatique à la 
mesure de la radioactivité de l’air fut suggérée en 1902 
par Sella, mais c’est Mlle Aliverti qui, en 1931, cons- 
truisit et mit au point le premier appareil fonctionnant 
sur ce principe [A1-2], [AI-3] ; l’effluve n’est plus créée 
par le gradient du champ au voisinage de l’électrode 
centrale, mais par une série de pointes dirigées vers 
l’électrode collectrice. Cette méthode fut reprise par de 
nombreux expérimentateurs [Ca-2], [Ma-1]. Wilkening 
[Wi-2] a utilisé ce procédé pour mesurer les dimensions 
des aérosols radioactifs de l'air. 

Deux électrodes courbes et coaxiales présentent l’in- 
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convénient de collecter les poussières sur une grande 
surface, et l’analyse radioactive devient plus difficile. 
On préfère souvent utiliser l’effluve produite entre une 
série de pointes et un plan pour que les poussières se 
rassemblent sur ce plan [Ch-21, [Pi-1]. Les différences 
de potentiel sont inférieures à 10 000 volts, la distance 
entre les pointes et la plaque étant de quelques milli- 
mètres. En remplaçant cette plaque par une bande 
d'aluminium en déroulement continu, on peut réaliser 
un enregistrement continu des variations de la radio- 
activité de l'air en un point donné[Be-4],[Be-5],[Ka-1] 
[Kn-1], [Wi-3] ; l'efficacité est en général excellente et 
voisine de 100 %. Bergstedt [Be-5] juge cependant que 
cette estimation est trop optimiste. 


Seconde catégorie. — Un atome qui vient d'émettre 
une particule alpha se transforme en ion positif ; le 
passage de la particule alpha dans le cortège électro- 
nique de l’atome modifie en effet profondément le 
champ électrique qui y règne et provoque le départ de 
plusieurs électrons. Les ions ainsi créés se déplacent 
dans atmosphère sous l'effet du champ électrique. On 
pourra les recueillir sur un fil ou une plaque chargés 
négativement. Cette méthode dite « d'activation » ou 
« de l’activité induite » a été souvent employée pour 
mesurer la radioactivité atmosphérique, dès les pre- 
mières recherches [As-1], [El-1], [Ja-2], et, plus récem- 
ment, par Evans et Godeman [Ev-2], Wilkening[Wi-4], 
Skorka[Sk-1].Vohra l’a utilisée pour mesurer la teneur 
de l'air en dérivés du thorium, et procéder à des repé- 
rages de gisements à distance [ Vo-1] ou à des mesures 
de contaminations atmosphériques [Vo-2}: autres 
auteurs ont employé ce procédé pour mesurer lés con- 
centrations de radon et de thoron au voisinage de 
sources thermales [Le-31, [Le-4], [Pr:3]. 

La méthode d’activation présente un inconvénient : 
les ions positifs créés dans l’atmosphère par désinté- 
grations alpha se fixent sur les aérosols de l’air, ce qui 
diminue leur mobilité ; ils peuvent aussi rencontrer des 
ions négatifs, également présents dans l’atmosphère, 
avec lesquels il se recombinent. Le nombre des parti- 
cules chargées collectées sera donc diminué du nombre 
des particules recombinées ; il faut tenir compte de cette 
recombinaison pour évaluer l’activité présente dans 
l’air. Certains des auteurs cités plus haut ont commis 
l’erreur de ne pas effectuer cette correction. 


49 Filtration. — On débarrasse l'air deses poussières 
en le faisant passer à travers une surface filtrante. Les 
méthodes de filtration sont appliquées industriel- 
lement, par exemple pour assurer la décontamination 
de l'air ayant séjourné dans les laboratoires de radio- 
chimie traitant des produits très actifs. Les caracté- 
ristiques de plusieurs modèles de filtres commerciaux 
ont fait l’objet d’études [Be-1], [Co-1], [Fi-2], [Ja-11. 

Le principe de la méthode est simple : l’air à dépous- 
siérer est aspiré à travers un filtre qui peut être cons- 
titué par un bourrage ou un feutrage de fibres de cel- 
lulose [St-4], de laine de verre, d'amiante, de matières 
végétales ou de matière plastique. Le diamètre moyen 
de chaque fibre va de quelques microns à quelques di- 
zaines de microns. 

La filtration des poussières provient : 1) d’un effet 
de tamisage en ce qui concerne les grosses particules ou 
les agglomérats de particules ; 2) d’un effet d'inertie, 
par déviation des filets d’air au passage sur les fibres ; 
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3) d’un effet de diffusion pour les particules les plus 
fines : les particules de diamètre inférieur à quelques 
dixièmes de micron, entraînées par le courant d’air, 
oscillent par agitation thermique autour d’une ligne 
moyenne ; l’amplitude des oscillations augmente quand 
les dimensions des aérosols diminuent. La probabilité 
de rencontre d’une particule avec une fibre, par effet 


F1G. 5. — EAR-600 vu de face. D'après [La-3]. 


À: Enregistreur MECI. — B : Panneau général de com- 
mande. — C : Panneau intégrateur. — D : Alimentation 
stabilisée des circuits électroniques (ALS). — E : Pan- 
neau libre. — F: Prise de masse, — G : Compteur 
horaire. 


de diffusion, augmente donc de la même manière ; 
4) pour certains filtres, s’ajoute un effet électrostatique 


1 


NS 


les fibres ayant gardé une charge électrique lors de leur 
compression. ; 

Dans le cas de la radioactivité atmosphérique, les 
filtres sont d’un type plus particulier : il s’agit, en 
général, de papiers filtres de quelques dixièmes de 
millimètre d'épaisseur. Il faut cependant citer une inté- 
ressante réalisation du Commissariat à l'Énergie Ato- 
mique : il s’agit d’un appareil portatif où l'air est filtré 
sur du charbon de bois qui retient également le radon 
[Ca-1]. Le comptage est assuré par un compteur de 
Geiger-Müller placé au centre de la cartouche de 
charbon. 


F1G. 6, — Vue arrière de l'EAR-600. D’après [La-31. 
À : Préampli GM (P,).—C: Préampli. — D: Détecteur GM. 


= E : Détecteur PM. — F : Dérouleur. G: Inter- ” 
rupteur de l'éclairage du coffre. — H : Manomètre à 
liquides. — l': Détendeur (caché en partie). J : régula- 
teur de tension. — K : Prise secteur. — LI, : Commuta- 
teur 110-220 V. — M : Trompe d'aspiration. — N : Bar- 
rette de distribution. — O : Aspirateur. — P. Coffre de 
prélèvement (ouvert), 
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Selon l'installation, le filtre peut être fixe [Ba-1], 
[Bo-1],[Da-1],[Ha-2],[Ha-3], [Ha-4], [Se-2], ou mobile, 
Dans ce dernier cas, on obtient un enregistrement con- 
tinu des variations de la radioactivité atmosphérique ; 
les appareils E. A. R. construits par le Commissariat à 
l'Énergie Atomique [Je-2], [La-1], [La-3], ainsi que de 
nombreux appareils étrangers [An-5], [Bo-1], [De-4], 
[Ku-1], [Ma-6], [Sa-41, [Sp-1], [Wa-2] rentrent dans 
* cette catégorie. On trouvera sur les figures 5 et 6 
V'EAR-600 du Commissariat à l'Énergie Atomique. 

La plupart des filtres ont une efficacité qui varie 
de 90 % à 100 % pour les poussières industrielles de 
diamètre supérieur à un micron, Cette efficacité n’est 
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Pair, préalablement filtré, est envoyé de façon continue 
dans une chambre de désintégration (1 000 litres). Les 
produits de désintégration du radon formés dans cette 
chambre sont recueillis sur un filtre et analysés. L’acti- 
vité minimale détectable est de 1,4.10—14 curie de 
radon par litre d’air. 


50 Mesure de l’activité des prélèvements. — Les aéro- 
sols radioactifs sont toujours recueillis sur un support 
et la moitié du rayonnement émis est dirigée vers l’in- 
térieur du collecteur ; elle ne pourra pas être comptée, 
sauf si le support est un papier filtre de faible épais- 
seur et si l’on opère par comptage bêta. 

La mesure de l’activité recueillie sur le collecteur 
peut se faire par comptage æ« ou par comptages f et y. 
Les aérosols radioactifs naturels sont émetteurs œ, Ê 
et y, alors que les aérosols radioactifs artificiels (pro- 
duits de fission, surtout) sont émetteurs $ et y ; un 
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plus que de 30 à 40 % pour des aérosols de dimensions 
voisines de 0,1 micron ; elle atteint à nouveau 90 % 
pour des gros ions de dimensions voisines de 0,001 mi- 
cron [Ja-1]. 

Nous pouvons classer dans cette catégorie notre 
appareil de dosage continu des gaz radioactifs [B1-2] : 
dans les méthodes qui précèdent, on suppose impli- 
citement que l’équilibre radioactif est atteint entre le 
radon et ses descendants, mais ce n’est pas toujours 
le cas, particulièrement en atmosphère agitée. Labeyrie 
et Weill [La-3] ont abordé ce problème en utilisant une 
chambre d’ionisation différentielle et en l’alimentant à 
travers un filtre (fig. 3). Dans notre appareil (fig. 7), 
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— Ensemble de dosage continu du radon dans l’atmosphère [BI-2]. 


comptage « permet de séparer les éléments naturels 
des éléments artificiels. Remarquons que, lorsque l’on 
mesure la radioactivité totale de l’atmosphère par 
comptage B, on peut séparer facilement les éléments 
naturels des éléments artificiels en étudiant la loi de 
variation des courbes de décroissance. 

Le premier appareil de mesure utilisé fut l’électro- 
scope. Il fut remplacé ensuite par la chambre d’ioni- 
sation, dans laquelle sont placés les fils ou les plaques 
chargés d’aérosols [Sa-1]. Cette méthode continue à 
être employée dans le cas de la précipitation électro- 
statique [AI-2], [AI-3] et dans celui de la filtration 
[Cr-1], [Fa-1], [Ma-3], [Sm-1]. 

Si la chambre fonctionne en impulsions, l'emploi de 
chambres spécialement étudiées permet de tracer le 
spectre des éléments recueillis par précipitation élec- 
trostatique ; c’est ainsi qu'ont opéré Facchini et Malvi- 
cini [Fa-1]}, [Fa-2]. 
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Quelques auteurs ont utilisé des compteurs propor- 
tionnels à l’intérieur desquels ils plaçaient l’électrode 
collectrice d’un précipitateur [Ca-2], [Wi-2]. Ce procédé 
est assez difficile, car le remplissage du compteur doit 
être renouvelé à chaque opération. 

Le compteur de Geiger-Müller est le détecteur le plus 
employé dans les mesures 8 ou y sur un collecteur. La 
mesure de l’activité recueillie sur un filtre peut se faire 
simultanément à l’aide de deux compteurs de Geiger- 
Müller placés de part et d’autre du filtre [Ph-1]. 
Quivy [Qu-1] a calculé quelle est, entre les deux comp- 
teurs, la différence du taux de comptage due à l’absorp- 
tion des rayons 6 à l’intérieur du filtre. 

Les appareils employant des scintillateurs sont aussi 
très nombreux. Certains auteurs utilisent simulta- 
nément un compteur à scintillations, pour la mesure 
des rayons &, et un compteur de Geiger-Müller, pour 
la mesure des rayons B, ces deux détecteurs étant 
placés de part et d’autre du filtre [Wa-2]. On peut 
employer, de même, deux compteurs à scintillations, 
l’un des scintillateurs étant une poudre de ZnS(Ag), 
l’autre un scintillateur plastique pour rayons £. 

Il est assez courant de placer les plaques ou le filtre 
à analyser au contact d’une émulsion nucléaire [Al-4], 
[AI-5], [Ch-1], [Fa-4], [La-2], [Le-1], [Ma-31, [Mo-1], 
[Mo-2], [Si-1]. On peut même recueillir directement 
sur la plaque photographique les aérosols, par préci- 
pitation [Er-1], [Je-1]. L’examen au microscope des 
traces obtenues dans les plaques développées donne 
des renseignements sur le nombre et sur la nature des 
particules recueillies sur le collecteur. 

On peut enfin faire appel à l’analyse chimique des 
poussières retueillies [Ki-1], [AI-1]. Avy et Raiïllère 
[Av-1], par exemple, emploient pour recueillir les pour- 
sières un filtre cristallin soluble. 


Causes diverses de variation de la radioactivité natu- 
relle de l’air. — DIFFUSION DE L’ÉMANATION. — La 
radioactivité naturelle de l'atmosphère varie beaucoup 
d’un lieu à un autre ; en un même lieu, on observe des 
variations fréquentes et importantes [Bl-4]. 

Cette radioactivité provient de la diffusion de radon 
et de thoron à partir du sol : elle est fortement in- 
fluencée par les conditions météorologiques des couches 
inférieures de l’atmosphère, et, en particulier, par 
l’existence de courants de convection au niveau du sol. 
Elle est liée à un phénomène d’échange dont il est 
possible de faire une étude quantitative. 

La fonction s(z), qui représente la distribution verti- 
cale du radon dans l’atmosphère, est la solution d’une 
équation différentielle de la forme [Co-3], [Is-1] : 


Ads 1dA ds _dE 
dei lodz d ‘°° d' 
s — densité de l’émanation (curie cm—ÿ). 
E — quantité d’émanation (curie em—? s—1), 
A — coefficient d'échange (g em—ts—t). 
À — constante radioactive du radon (s—1). 
p_— densité de l’air. 


Cette équation reste valable pour le thoron. Elle a été 
intégrée en supposant dA/dz — 0 [He-3], ou en 
donnant à À des valeurs arbitraires. Les résultats 
obtenus s'accordent grossièrement avec la valeur expé- 
rimentale de Æ : 


Pa=D0PI0E Cire ces (Is-1) 
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La concentration optimale a été obtenue par 
Priebsch [Pr-2] ; dans certains cas, les écarts sont im- 
portants : à New-York, par exemple, Cullen a trouvé 
E = 5.10—18 curie em—? sec—1 [Cu-1]. Ces écarts 
semblent dus à la nature des sols. 


sk carf midi RÉ 


INFLUENCE DES CONDITIONS MÉTÉOROLOGIQUES. — 


L’équation précédente n’est valable qu’en atmosphère 
calme (absence de vent, pas d’inversion de tempé- 
rature). 

L'influence du vent se manifeste de deux manières 
[Ch-2], [Fu-1], [He-2], [Hu-2]: la vitesse d'émission 
du radon est augmentée, ce qui entraîne une augmen- 
tation de la radioactivité ; d’autre part, le vent, emme- 
nant au loin les particules formées, provoquerait une 
diminution de leur concentration. Enfin, la direction 
du vent joue un rôle important [Ch-21, [Co-2], [Je-31, 
[Je-4], [Re-1], [Sk-1], [St-1], [Ta-1], et l’on a pu ainsi 
procéder au repérage à distance de gisements de tho- 
rium [Vo-1]. 

Sur la diffusion à partir du sol, l'influence de la 
température est très faible : il faut que la température 
s'élève de 1200 à 150 0C pour que la quantité de gaz 
émanant des roches soit doublée [Gi-1]. Par contre, son 
influence est très importante dans les variations quoti- 
diennes de la radioactivité de l’air. En général, en 
atmosphère calme, on trouve un minimum d’activité 
dans l’après-midi ou à la tombée de la nuit, au moment 
où les courants de convection de l’atmosphère sont 
favorisés par l’échauflement du sol, et un maximum 
durant la nuit ou aux premières heures de la matinée, 
alors que l’air froid demeure au voisinage du sol [Wi-41]. 
Se basant sur les données de Lettau [Le-5], Cotton 
[Co-3] a fait le calcul théorique du rapport des activités 
diurne et nocturne : ce rapport est d’autant plus grand 
que la radioactivité moyenne est plus faible. 

D'autre part, en un lieu donné, la radioactivité est 
maximale lors des inversions de température, c’est-à- 
dire quand les basses couches de l’atmosphère se main- 
tiennent constamment au voisinage du sol. 

La pluie provoque une décroissance de l’activité 
naturelle ; elle entraîne vers le sol les aérosols en sus- 
pension dans l’air, elle colmate les fissures et les capil- 
laires du sol par lesquels s’échappe l’émanation [Ko-1]. 
La neige et le gel produisent le même effet. Le brouil- 
lard retient au contraire les particules au voisinage du 
sol, ce qui augmente l’activité de l’air ; les aérosols 
formés ont aussi des dimensions plus grandes et leur 
probabilité de capture est augmentée. 

La radioactivité dépend de la pression atmo- 
sphérique [Sa-2] : le radon se dégage d’autant plus 
facilement que la pression atmosphérique est plus 
faible [He-2]. 

Enfin, on note une variation saisonnière marquée, 
présentant un maximum annuel, mais la période de 
l’année à laquelle se manifeste ce maximum est variable 
d’une région à une autre [Bl-4], [Wi-5]. 

Il est très difficile d'établir une théorie globale de la 
variation de la radioactivité naturelle de atmosphère, 
en raison de la complexité des facteurs météorologiques 
qui interviennent dans les observations. 


Ordre de grandeur de la radioactivité naturelle, — 
Compte tenu des divers facteurs que nous venons de 
citer, nous pouvons donner la valeur moyenne de la 


 No1 


radioactivité atmosphérique. Sur les continents, et au 
niveau du sol, elle est de l’ordre de 1 à 2.10—13 curie 
par litre d’air pour le radon, et de 2 à 3.10—15 curie 
par litre pour le thoron. 

Au-dessus des océans, les valeurs trouvées sont 
beaucoup plus faibles : 10—15 curie par litre pour le 
radon, 0,1.10—15 curie par litre pour le thoron [BI-4], 
[He-4], [He-5], [Is-17, [Sk-11. 

En montagne, on obtient des résultats plus forts, 
alors qu’à haute altitude, les mesures aériennesfaites par 
Wilkening [Wi-4] et par Garrigue [Ga-4], [Ga-2] mon- 
trent que la radioactivité décroît très vite quand l’alti- 
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tude augmente : au-dessus de 3 000 mètres, on ren- 
contre surtout des ions libres, non recombinés. 

La dose de tolérance est plus de dix fois supérieure 
aux plus fortes concentrations enregistrées. Il n’en est 
plus de même dans les mines d'uranium où des instal- 
lations de contrôle fonctionnant en permanence sont 
indispensables. 

Nous remercions, Mile Huguette Bourdès, Secrétaire 
du Centre de Physique Nucléaire, pour sa participation 
au dépouillement de la documentation. 

Cette étude a pu être menée à bien grâce à un contrat 
du Commissariat à l'Énergie Atomique. 


Manuscrit reçu le 27 juillet 1960. 
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MODULATION DE LA LUMINESCENCE 
EXCITÉE PAR DES RAYONS « 
DANS CERTAINS SULFURES DE ZINC 
SOUS L'ACTION DE CHAMPS ALTERNATIFS 


Par R. HENCKk et A. COCHE, 


Département de Chimie Nucléaire, 
Centre de Recherches Nucléaires, Strasbourg. 


Différents auteurs (Halsted [1], Daniel et coll. [21 
ont montré que la lumière émise par certains sulfures 
de zinc, excités par des rayonnements ultra-violets et 
soumis à des champs alternatifs, était modulée à la 
fréquence du champ au à la fréquence double, suivant 


la longueur d’onde de la lumière excitatrice. Jaffe a 


observé également ce phénomène lors d’une irra- 
diation par des rayons cathodiques [3]. 

Nous avons étudié les ondes de luminance qui appa- 
raissent en appliquant des champs alternatifs (50 Hz) 
à des sulfures mixtes de zinc et de cadmium (88 % SZn, 
12 %, SCd) activés au manganèse en proportion va- 
riable et excités par les rayons &« d’une source de ?21Po 
(0,5 mC). Les concentrations en manganèse des diffé- 
rents sulfures (1) sont comprises entre 10—? et 
1025 MnieS/nCd:A (107), B (5:40-%); C (10°), 
D (5.10—4), E (104). Les poudres luminescentes, dis- 
persées dans de l’araldite (épaisseur environ 100 p), 
sont placées entre les armatures d’un condensateur 
plan : l’une est constituée par une lame de verre con- 
ducteur transparent, l’autre par une couche métal- 
lique (argent ou aluminium) évaporée sous vide et 
suffisamment mince pour laisser passer les rayons «. 

Cette cellule est posée sur un photomultiplicateur 
(Dario 53 AVP) dont l’anode est reliée à l’un des 
amplificateurs (à courant continu) d’un oscilloscope à 
deux voies. La deuxième voie permet d’observer, en 
même temps que l’onde de luminance, la tension appli- 
quée à la cellule. On a examiné successivement les 
phénomènes dans les deux bandes d’émission (bleue 
et orange) des scintillateurs, en plaçant entre la photo- 
cathode et le condensateur luminescent, des filtres 


(1) Fournis par les Établissements Massiot-Fluor. 


interférentiels, dont les maximums de transmission se 
situaient respectivement à 4 790 et 5 840 À. 

La figure 1 montre les résultats obtenus dans la 
bande orange (a) et dans la bande bleue {b), avec le 


Rre1Na: F1G. 1 b. 
Scintillateur D. Scintillateur D. 
Bande orange. Bande bleue. 

1. —— Onde de luminance ; 2. — Excitation « ; 
3. — Tension appliquée. ; 4. — Niveau zéro. 


seintillateur D (5.10—4 g/g Mn). Dans ce cas, l’électro- 
luminescence produite par la tension appliquée (200 
volts, 50 Hz) est négligeable. Ces ondes de luminance 
présentent, dans les deux bandes d’émission, deux 
maximums par période du champ. 

Dans la bande de plus grande longueur d’onde, pour 
tous les scintillateurs étudiés, la valeur minimale Lin. 
de la luminance, atteinte quand la tension passe par 
zéro, est supérieure à la valeur Z, de la luminescence 
produite par les rayons « en l’absence de champ. Les 
maximums sont déphasés (légèrement en retard) par 
rapport à ceux de la tension, le maximum princi- 
pal Lirax1, Se produisant quand l’électrode métallique 
traversée par les rayons « est négative par rapport au 
verre conducteur. Jaffe a trouvé un effet identique (3), 
alors que Halsted observe le maximum de lumière 
lorsque l’électrode métallique est positive. Ces résultats 
sont confirmés en utilisant une tension obtenue par 
redressement d’une alternance (fig. 2a) : dans l’inter- 
valle t,, t3, la luminance reste supérieure à Z, et l’émis- 
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sion de lumière est maximale lorsque la valeur instan- 
tanée de la tension l’est également. 


Rrcmorar FG.22 0: 
Scintillateur C. Scintillateur E. 
Bande orange. Bande bleue. 
4. — Onde de luminance ; 2. — Excitation « ; 

3. — Tension appliquée.; 4. — Niveau zéro. 


Dans la bande bleue, les phénomènes sont différents : 


la valeur moyenne de l’onde de luminance est infé- 


rieure à Z,, les maximums sont atteints lorsque la 
tension est proche de zéro, tandis que les minimums 
sont en légère avance de phase par rapport aux maxi- 
mums de la tension, et le minimum principal Lin.1 Se 
produit lorsque l’électrode métallique de la cellule est 
négative. Avec la tension redressée utilisée précé- 
demment, on observe (fig. 2b) des minimums en phase 
avec la tension. Ceci montre que le mécanisme d’émis- 
sion de la lumière bleue est différent de celui de l’émis- 


sion orange. 
L'étude des variations relatives des grandeurs défi- 
; ; E —L 
nies sur les figures 1, a et b, soit ALymax. — y 
œ 
Lin. — L ; , 
MAUERRURE nr — en fonction de la tension 


œ 
appliquée au scintillateur, fait apparaître des diffé- 
rences caractéristiques, suivant la bande d’émission 
considérée. A titre d'exemple, on a représenté sur les 


AL 


200 (its) 200 


Prices 
Scintillateur D, 
Bande orange. 


Fic, 3 b: 
Scintillateur D. 
Bande bleue. 


figures 3a et 3b les résultats trouvés avec le scintil- 
lateur D mentionné précédemment. 
Dans la bande orange, les diverses grandeurs ALyax. 


N°10 


et ALyin. Sont positives (en ce qui concerne Ain, 
ceci n’est le cas qu’au-dessus d’une certaine tension 
— 70 volts dans l'exemple considéré — pour les scintil- 
lateurs ayant une concentration en manganèse infé- 
rieure ou égale à 5.10—4g/g) ; elles augmentent lorsque 


la tension croît, mais tendent vers une limite aux | 


tensions élevées, d’autant plus rapidement que la 
concentration en manganèse est plus grande. La valeur 
maximale de ALyax. est obtenue avec le scintillateur ” 
dont la teneur en manganèse est de 10—* g/g ; elle 
est de l’ordre de 5 à 300 volts. 

Dans la bande bleue, ALyax. est positif, croît avec 
la tension, mais tend vers une limite qui est atteinte 
d’autant plus rapidement que la concentration en man- 
ganèse est plus faible. AZyin. est négatif et diminue au 
fur et à mesure que la tension augmente. 

Si nous définissons un pourcentage de modulation M 


| j Max1. ___ Min.1 
par la relation : M — 100 Les Eee: 


vons des valeurs du même ordre de grandeur que celles 
signalées par Halsted (ces dernières étant obtenues 
pour une modulation à la fréquence du champ, alors 
que, dans nos expériences, elle a lieu à une fréquence 
double), soit 50 % dans les conditions les plus favo- 
rables. 

La figure 4 montre, pour les différents sulfures A, 
B, C, D, E, la variation de M, en fonction de la ten- 
sion, dans les deux bandes d'émission. On voit que M 


, Nous trou- 


augmente d’abord avec la tension, puis passe par un 


maximum plus ou moins accentué, ou tend vers une 


FIG 4 


limite. Les valeurs maximales s’observent dans la 
bande bleue des scintillateurs ayant des concentrations 
en manganèse faibles (10—4 et 5.10—1 g/g), où elles 
sont d'environ 50 %. Dans la bande orange, la modu- 
lation maximale ne dépasse pas 40 /. 
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._ SUR QUELQUES ASPECTS 
_ DE TRACES DE PARTICULES IONISANTES 
DANS DES ÉMULSIONS PHOTOGRAPHIQUES 
PEU SENSIBLES 


Par Marie ADER 
avec la collaboration technique de 
Marie Paule HouARD-CABANNES, 


Laboratoire de Physique Atomique et Moléculaire 
du Collège de France, 
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Nous utilisons des émulsions Ilford XO, X-1, X-2 et 
ces mêmes émulsions imprégnées de bore. On sait que 
ces dernières peuvent enregistrer des protons de plu- 
sieurs MeV (1). 

Nous montrons des traces de protons et de parti- 
cules de parcours anormalement long du Po, en inci- 
dence presque rasante dans des émulsions X0 et dans 
des émulsions X-1 et X-2 imprégnées de bore. 

Un examen attentif des émulsions ainsi irradiées (2) 
sugeèrerait quelques remarques : 
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Fi, 1.— Dans l’émulsion X,, des fins de traces de protons. 


Les traces de particules de faible énergie sont formées 
de segments denses, rectilignes, séparés par des lacunes 
souvent très petites et par suite peu apparentes ; les 
traces de particules d’énergie plus élevée sont des files 
de grains interrompues par des lacunes. Dans 70 % 
environ des cas, un changement de direction de la 
trace est précédé d’une lacune apparente ; souvent, et 


(1) 2e Colloque International de Photographie Corpus- 
culaire, Montréal, 1958. Les Presses Universitaires de 
Montréal, 1959, p. 264. 

. (2) Les irradiations ont été faites à l’I. N. S. T. N. à 
Saclay. 
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Les clichés 2 et 3 montrent des traces de protons dans une 
émulsion Æ-1 imprégnée de bore. Les déviations et 
variations de densité de grains après lacune sont appa- 
rentes. 


EG, & 
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ceci est remarquable surtout dans les traces de parti- 
cules de plus grande énergie, une lacune marque un 
changement de densité de grains, comme si le phéno- 
mène de capture suivi d'abandon, entraînait une 
variation sensible de l’énergie de la particule et, peut- 
être, l'émission d’une autre particule non enregistrée 
par l’émulsion ou neutre, emportant le supplément 
d'énergie. 


Fren5: 


FIG: 


Les Clichés 4, 5, 6 et 7 représentent des traces de parti- 
cules de parcours anormalement longs du polonium dans 
les émulsions Æ-1 et Æ-2 sensibilisées au bore. Les 
changements de densité de grains et'les déviations 
après lacune sont nettes. 


Lettre reçue le 8 novembre 1960. 


Not | 
POUVOIR D’ARRÊT DE L’OR 
POUR DES PARTICULES ALPHA DE 14,36 MeV 


Par C. GéRaRDIN, R. BiLWwESs et D. MAGNAG-VALETTE, “| 


Département de Physique Corpusculaire, C. R. N., 
Strasbourg-Cronenbourg. “| 


Pour détecter les particules chargées et mesurer leur … 
énergie, les physiciens utilisent souvent des feuilles 
métalliques dont le pouvoir d’arrêt est supposé connu. 
Or, nous avons constaté des différences appréciables . 
d’un auteur à l’autre. C’est ainsi que pour absorber 
des particules alpha de 14 MeV par l’aluminium, on 
trouve dans la littérature 29,4 mg/em? [1] et 
27,6 mg/cm? [2]. Les auteurs procèdent généralement 
de la manière suivante : à l’aide de considérations 
théoriques, ils déterminent la courbe d£ [dx et tracent 
la courbe intégrale dont ils ajustent les ordonnées à 
l’aide d’un point expérimental. . 

Pour nos expériences, nous avions à utiliser des 
feuilles d’or,.et ne disposions que d’une seule courbe [1] 
et d’un point expérimental à 4,5 MeV [3]. Ainsi, 
à 4,5 MeV, les épaisseurs de l’absorbant étaient de | 
14,6 mg/em? [1] et de 17 mg/em? [3]. Nous avons 
tout d’abord pensé que la valeur théorique de d£ [dx 
ne coincidait pas avec les valeurs expérimentales. , 
Notre article rend compte de cette vérification entre 9 # 
et 14 MeV. | 

A l’aide d’un spectrographe magnétique du type 
Buechner, associé au Van de Graff de 5,5 MeV de «… 
Strasbourg-Cronenbourg, nous avons mesuré l’énergie 
résiduelle des particules alpha de 14,36 MeV, après 
passage à travers des absorbants d’or d’épaisseurs M 
diverses. Les alpha provenaient de la réaction « 
10B(d, «)$Be, pour une énergie incidente Æa — 5 MeV 
et correspondaient à la formation de l’état fonda- 
mental de îBe. : 

La courbe d’étalonnage du spectrographe dans la 
partie qui nous a intéressés, est donnée figure 1. Nous 
avons maintenu les valeurs du rayon de courbure au 
voisinage d’une valeur donnée en modifiant le champ. 
La précision n’est pas la meilleure qu’on puisse obtenir 
avec le spectrographe, mais elle est largement sufi- 
sante, comparée au «straggling » important des parti- 
cules alpha. La fréquence de résonance nucléaire y 
utilisée pour la mesure du champ est une fonction 
bien connue de p?/2m : vie — K(p?[2mxey, p étant 
la quantité de mouvement des particules &« et m leur 
masse. La courbe 1 donne Æ en fonction de la dis- 
tance + des particules le long de la plaque qui les 
recueille. Terreur relative sur l'énergie due à l’étalon- 
nage est : AFJE = 3.10-$, 

Il faut ajouter à ceci l'erreur relative sur l’énergie 
des particules alpha de 14,36 MeV, erreur due à 
l’épaisseur de la cible et à l’angle solide d'émission. 
Nous n’avons pas à tenir compte de la largeur natu- 
relle du niveau fondamental de $Be qui est très faible 
(environ 5 eV). Nous avons estimé à 45 keV la demi- 
largeur à mi-hauteur du pic « sans absorbant. Nous 
nous sommes servis de la courbe d'étalonnage pour 
mesurer l’énergie de ce pic. En effet, la construction 
du spectrographe ne permet pas une mesure précise de 
l’angle d'émission et ainsi, l'énergie des alpha ne peut 
être déduite des données cinétiques de la réaction. 
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L'erreur relative sur l’énergie des & sans absorbant 
est alors de 4,5.10— et l’énergie de 14,36 MeV + 0,06. 
La courbe 2 représente la variation de l’épaisseur e 
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Fic. 1. —- Étalonnage du spectrographe. 


de l’absorbant en mg/em? en fonction de l’énergie 
résiduelle Æ; des particules «. Les barres horizontales 
indiquent la largeur du pic après absorption. Le tableau 


e épaisseur d'or traversce 
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HrC20 2: 


donne en même temps que les valeurs de Æ;, celles 

<E—E >? 
Es 

_ Comme on le voit, ces valeurs sont de l'ordre du 

« straggling » moyen à ces énergies (environ 210). 


de la fluctuation moyenne relative <l 
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TABLEAU 
n\ GG , CR 
RE A A CN RE pq 
EN mg/cm? EN MeV Et 
0 14,36 + 0,06 

6,72 + 0,05 13,45 + 0,14 0,9 

10,14 + 0,12 12,96 + 0,15 1 

12,96 + 0,17 12/47 220,22 4,7 

15,90 + 0,36 12,13 + 0,22 1,8 

17,56 + 0,14 11,824 0,25 | 

21,85 200,35 41,10 + 0,24 2,2 

DE ES 10:22 10,60 + 0,28 2,6 

25,76 + 0,30 10,56 + 0,29 2 

28,63 + 0,92 9,85 + 0,35 3,5 

30,61 + 0,57 9,61 + 0,38 3,9 


Soient À, le parcours des « de 14,36 MeV dans l'or, 
e l'épaisseur de l’absorbant et À, le parcours résiduel, 
nous avons À, — À + e. La courbe 2 donne e = g(Æ;). 
Nous avons donc 4% = R, — g(E;) = f(Ei), la 
courbe À — f(E) étant la courbe parcours-énergie 


(courbe 3). Nous voyons qu’il suffit de faire une 
R 


parcours en M9 /<m2 
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——— D'APRES WHALING 


F- NOS POINTS EXPERIMENTAUX 


symétrie de la courbe 2 par rapport à l’axe des énergies 
pour obtenir une courbe déduite de la courbe (3), par 
une translation dirigée suivant l’axe des ordonnées. 
Nous avons ainsi vérifié que la pente de la courbe de 
Whaling coïncide sensiblement avec la pente de notre 
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courbe, sauf en ce qui concerne les points à 9,61 et 
9,85 MeV. 

En effet, le pic des particules &« correspondant à 
l’état excité de 2,9 MeV du fBe s’élargit en même 
temps que le pic du fondamental, ce qui provoque un 
déplacement du maximum de ce dernier vers les basses 
“énergies. Nous n’avons donné ces points qu’à titre 
indicatif. 

Le point expérimental choisi par Whaling pour 
ajuster sa courbe est de détermination ancienne [4]. 
Il est visiblement en contradiction avec le résultat de 
Mne Faraggi [3]. Il nous reste donc à contrôler l’or- 
donnée absolue des deux points pour expliquer la 
divergence des résultats initialement constatée. 
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